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PART  I:  DETECTION  IN  GAUSSIAN  MIXTURE  NOISE 

1.0  INTRODUCTION 

1.1  BACKGROUND  SUMMARY 

In  many  applications,  the  passive  detection  of  signals  in  noise 
is  not  adequately  performed  if  the  detector  design  is  based  on  Gaussian 
noise  assumptions.  Particularly  in  the  LF  and  VLF  portions  of  the  electro¬ 
magnetic  spectrum  and  in  the  acoustic  environment,  the  ambient  noise  observed 
by  the  receiver  contains  impulsive  components  which  cause  the  noise  probability 
distribution  to  depart  from  the  Gaussian  model.  As  a  result,  the  sensitivity 
of  a  Gaussian-based  detector  is  degraded  for  given  false  alarm  probability. 
Therefore,  much  effort  has  been  focused  on  realizing  improved  detector  and 
receiver  designs  that  are  matched  to  realistic  noise  distributions,  using 
measured  values  of  noise  parameters. 

The  detector  design  problem  is  complicated  by  the  fact  that  the 
parameters  describing  the  noise  distribution  are  difficult  to  measure  and 
vary  with  time.  Since  it  is  not  practical  in  these  situations  to  employ 
theoretically  "optimum"  detectors  which  rely  on  a  priori  knowledge  of  noise 
parameters,  current  research  is  aimed  at  discovering  effective  detection  pro¬ 
cedures  which  are  "robust",  or  relatively  insensitive  to  uncertainties  in  the 
values  of  distributional  parameters,  if  not  "nonparametric" ,  or  "distribution- 
free".  Also  of  interest  are  realistic  noise  models  whose  parameters  are 
easily  measured  physical  quantities. 

Noise  due  to  man-made  sources  can  degrade  detection  and  communica¬ 
tion  performance  in  every  portion  of  the  spectrum.  Impulsive  interference 
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due  to  multiple-access  communications  emmissions,  for  example,  results  in 
a  non-Gaussian  received  noise  distribution.  The  threat  of  intentional  inter¬ 
ference  or  jamming  has  forced  system  designers  to  take  this  possibility  into 
account  by  modifying  their  detection  procedures.  Without  formulating  the 
problem  as  a  non-Gaussian  noise  problem  per  se,  many  heuristic  receivers  have 
been  developed  to  make  the  detection  and  communication  systems  operate 
satisfactorily  whether  jammed  or  not  jammed. 

The  success  of  certain  heuristic  detectors  in  adapting  to  non- 
Gaussian  noise  implies  that  their  designers  have  made  correct  assumptions 
concerning  the  noise  process.  For  the  most  part,  these  practical  detectors 
are  based  on  recognizing  the  "short  term"  or  time  varying  nature  of  the  noise, 
rather  than  its  "long  term"  or  marginal  (unconditional)  distribution.  There¬ 
fore,  it  seems  likely  that  modelling  non-Gaussian  noise  as  nonstationary  noise 
would  yield  better  results  than  the  present  theoretical  emphasis  on  marginal 
non-Gaussian  noise  models. 

In  order  to  show  the  motivation  for  our  efforts,  in  the  following 
paragraphs  we  give  further  background  on  non-Gaussian  modelling,  on  receiver 
design,  and  on  jamming  interference. 

1.1.1  Non-Gaussian  Noise  Models 

The  modelling  effort  evidenced  in  the  current  detection  literature 
has  been  approached  from  two  directions:  empirical  and  theoretical. 

A  good  example  of  the  empirical  noise  modelling  approach  is  the 
work  reported  by  Fennick  [1],  who  gathered  statistics  of  telecommunication 
channels  and  postulated  on  the  basis  of  measured  pdf's  that  the  underlying 
distribution  could  be  characterized  by  the  sum  of  a  Gaussian  pdf  and  an 
exponential  term: 
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pdf  =  Gaussian  +  exponential 

=  a  exp(-x2/2a2)  +  b  exp(-c|x|). 


(1.1-1) 


In  this  manner,  the  "larger  than  Gaussian"  tails  of  the  observed  statistics 
were  accounted  for. 

The  theoretical  approach,  one  working  from  the  physics  of  the 
situation,  can  be  exemplified  by  the  work  of  Hall  [2],  who  supposed  that 
the  received  noise  is  a  narrowband  Gaussian  process  multiplied  by  a  time- 
varying  weighting  factor.  Hall's  pdf  was  found  to  take  the  form 


-(m+l)/2 


(1.1-2) 


pdf  =  const,  (x2  +  a2) 


Good  fits  to  data  were  reported,  using  certain  values  of  the  parameters. 
However,  the  parameters  themselves  were  not  identified  with  the  physical 
processes  because  of  simplifications  chosen  to  make  the  mathematics  tractable. 


Middleton  [3,  4,  5]  asserts  that  he  has  found  tractable  pdf's 


that  fit  known  data  well  by  approaching  the  problem  from  the  physical/ 
theoretical  point  of  view,  and  that  the  parameters  of  the  distribution 
retain  physical  interpretations  (therefore  being  suitable  for  measurement). 
Moreover,  the  "Class  A,  B,  and  C"  noise  models  Middleton  has  developed  are 
claimed  to  be  "canonical",  or  capable  of  generating  the  wide  variety  of  ob¬ 
servable  statistics  while  keeping  the  same  functional  form.  For  example, 
the  pdf  for  Middleton's  Class  A  noise,  for  which  the  noise  bandwidth  is  said 
to  be  less  than  that  of  the  desired  signal,  takes  the  form  of  an  infinite 
sum  of  weighted  Gaussian  pdf's  with  increasing  variances: 


oo 


(1.1-3) 
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where  r  is  the  ratio  of  the  power  in  the  Gaussian  portion  of  the  interference 
to  that  in  the  impulsive  (Poisson)  component,  and  A  is  called  the  "impulsive 
index,"  a  kind  of  counting  function  for  the  impulsive  interference  and 
related  to  the  amount  of  overlap  in  individual  interference  waveforms 
(large  A  corresponds  to  a  trend  toward  Gaussian).  An  excellent  approximation 
for  small  values  of  A  is  given  by  Middleton  [5]  and  by  Vastola  [16]  as 


pdf  =  a  exp(-x2/2a2)  +  b  exp(-x2/2a|) , 


(1.1-4) 


or  the  weighted  sum  of  two  Gaussian  density  functions.  Although  we  have 
simplified  the  notation  somewhat  in  this  presentation,  each  of  the  parameters 
is  given  a  physical,  measurable  interpretation  by  Middleton. 


1.1.2  Non-Gaussian  Detector  Designs  (known  signal) 

If  N  independent  samples  jx..j  of  a  received  waveform  x(t)  are  to 
be  tested  as  to  whether  x(t)  contains  noise  only  or  known  signal  s(t)  plus 
noise,  the  log-likelihood  ratio  takes  the  form 


logA(xls  x2. 


N 

xN)  =  Y.  109  Pdf(xilsi)/Pdf(xilsi=°) 

i=l  *■ 


s+n 

^  threshold  (1.1-5) 
n 


where  x.  =  x(t. )  and  s.  =  s(t. ).  For  stationary  Gaussian  noise  the  resulting 
statistical  test  is  the  linear  detector 


N 


I 

■5=1 


Vi 


s+n 

> 

< 

n 


threshold; 


(1.1-6) 
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the  receiver  needs  only  to  perform  a  weighted  sum  (filtering)  of  the  input 
data  and  to  compare  the  value  of  that  sum  to  a  threshold  value. 

If  the  noise  is  not  Gaussian,  the  likelihood  function  may  be 
difficult  to  interpret  in  terms  of  a  discrete  component  or  analog  implementa¬ 
tion,  depending  on  the  form  of  the  noise  pdf.  For  Hall's  pdf,  (1.1.2),  the 
closed  form  permits  solving  directly  for  the  optimum  receiver  structure  as 


threshold 


(1.1.7) 


where  a  is  the  parameter  shown  in  (1.1-2).  For  pdf's  which  are  not  closed 
forms,  implementation  may  be  performed  using  digital  processing;  however,  analysis 
may  require  making  some  approximations.  For  example,  even  for  the  two-term 
approximation  to  Middleton's  pdf  given  by  (1.1-4),  the  likelihood  ratio  is 


logA 


Z  109 


1-1 


a  exp[-  (xi  -  si)2/2af]  +  b  exp[-  (x^  -  si)2/2a|j 

a  exp(-x?/2af)  +  b  exp(-x?/2a§) 


(1.1-8) 


A  way  out  of  the  analytical  difficulty  which  has  been  used  extensively  is  to 
treat  the  special  case  of  weak,  or  "threshold"  signals  and  therefore  to  obtain 
what  are  termed  locally  optimal  or  threshold  receivers.  How  this  approach 
works  may  be  explained  as  follows:  since  the  signal  is  "small",  to  a  good 
degree  of  approximation  the  pdf  may  be  written  as  the  pdf  for  no  signal 
plus  a  first  order  term  in  a  Taylor  series  expansion  ^5,  6]  ,  giving 

A.  .  p(x>  liil  .  i  -  it.  p'(^ 

1  p(xi  -  0)  p(x.) 


=  1  -  s. 


—  hog  pU,) 

3X .  L 


(1.1-9) 
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This  results  in  the  statistical  test 


N  N 


Z  Si  yi  =  E  Si  9(Xi) 


<  threshold  =  n» 


(1.1-10) 


i=l  i=l 


which  resembles  the  "linear"  test  (1.1-6)  after  the  data  has  been  transformed 
by  the  nonlinearity 


(1.1-11) 


The  form  of  this  nonlinearity  is  highly  sensitive  to  the  parameters  and 
functional  form  of  the  assumed  pdf,  as  illustrated  in  Figure  1.1-1.  In  that 
figure,  we  observe  that  for  Gaussian  noise,  g(x)  is  simply  a  linear  depend¬ 
ence,  whereas  for  the  other  assumed  distributions  shown  the  transformation 
can  be  almost  any  form,  depending  on  the  shape  and  parameters  of  the  noise 
pdf. 


The  form  of  the  nonlinearity  requires  either  knowledge  of 


adaptation  to  the  noise  conditions  which  exist.  Therefore,  Middleton  [5,  8] 
stresses  the  correspondence  between  the  parameters  of  his  pdf  model  and 
measurable  quantities.  The  ability  of  the  receiver  to  perform  satisfactory 
adaptation  to  time-varying  noise  parameters  can  make  the  difference  between 
success  of  failure  for  practical  non-Gaussian  detectors  [9],  Martinez  and 
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x(t) 

o — 


NOISE  DENSITY 

Gaussian 

p(x.j)  =  K  expj-x?/2a?j 


TRANSFORMATION 


x./ct? 

v  l 


Xi 


DOUBLE  EXPONENTIAL 
p(Xi)  =  k  exp|-3i  |xi  |J 


3isgn(xi) 


Bi 


~Bi 


xi 


BISTATIC  WAVEFORM  IN  GAUSSIAN 
p(x.)  =  K^x.)  +  e2(xi)| 

e2  =  expJ-(xi-A)2/2a?J 

e2  =  expj-(xi+A)2/2a?J 


FIGURE  1.1-1.  LOCALLY  OPTIMUM  NON-GAUSSIAN  DETECTOR  STRUCTURE 
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Thomas  [lOj  ,  for  example,  report  that  sampled  Artie  under-ice  ambient 
acoustic  noise  data  reveal  a  nonstationary  process  with  impulsive  components; 
this  behavior  is  well  illustrated  by  the  time  history  of  sample  variance 
shown  in  Figure  1.1-2.  Because  the  impulsive  components  occur  relatively 
infrequently,  large  sample  sizes  are  required  to  determine  accurate  estimates 
of  impulsive  parameters.  This  requirement  conflicts  with  the  need  to  take 
time-varying  noise  properties  into  account  over  a  smaller  sample  size. 

Since  the  form  of  the  nonlinearity  for  locally  optimum  detection 
is  often  complicated  as  well  as  sensitive  to  the  accuracy  of  measured  noise 
distribution  parameters,  investigations  have  been  made  to  determine  whether 
simplified  or  approximate  versions  of  the  transformations,  which  use  fewer 
parameters,  can  be  used  with  success.  For  example,  Miller  and  Thomas  £llj 
report  that  relatively  simple  piecewise  linear  approximations  can  give 
nearly  as  good  asymptotic  performance  (relative  to  a  Gaussian  detector)  as 
the  locally  optimal  nonlinearity.  Among  the  simplified  nonlinearities 
they  studied  were  the  "amplifier-limiter"  (also  known  as  a  "clipper1  or 
"soft  limiter"),  the  "hard  limiter",  and  the  "noise  blanker"  approximations 
to  the  optimum  nonlinearity  for  detection  in  a  Gaussian-Laplace  mixture 
noise  distribution.  These  transformations  are  shown  in  Figure  1.1-3. 
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SAMPLE  VARIANCE 


0.6 


FIGURE  1.1-2. 


AMPLE  VARIANCE  VS  RECORD  NUMBER  , 
inTTP  i!MnCD_  T  PF  AMRTFNT  ArnUSTIC  NOISE 


(FROM  [10]) 
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(a)  Amplifier-limiter 


FIGURE  1.1-3.  SIMPLIFIED  NONLINEARITIES  FOR  LOCALLY  OPTIMUM  DETECTORS 
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We  note  that  in  equation  (1.1-9)  above,  if  the  true  signal 

values  are  the  same  (s^  =  s,  Vi),  the  receiver  design  becomes 
N  N 

s  Y.  9  (xi}  <  n  or  Y.  9(xi}  <  ?  =  n'  *  (1-1-12) 

i=l  i=l 

The  receiver  does  not  use  any  information  on  the  value  of  the  (constant) 
signal,  since  the  thresholdn"  is  to  be  determined  by  false  alarm  probability 
requirements.  Therefore,  an  advantage  of  threshold  or  weak  signal  locally 
optimum  detector  designs  is  that  they  do  not  require  prior  knowledge  of  the 
signal  parameters  if  they  are  "slowly  varying"  over  the  N  samples  of  the 
input.. 

In  general,  the  log-likelihood  ratio  for  known  signal  values 
js^j  depends  strongly  on  these  values.  Therefore,  a  penalty  in  loss  of 
performance  can  be  expected  if  the  actual  received  signal  values  differ 
(for  example,  due  to  loss  of  synchronization).  This  statement  is  true 
even  for  a  constant  signal  (s^=s)  except  for  the  special  case  of  Gaussian 
noise,  as  seen  from  equation  (1.1-6),  so  we  expect  to  see,  for  example, 
close  agreement  of  detection  performance  between  locally  optimum  and  optimum 
detectors  when  the  signal  is  actually  weak,  and  a  loss  in  the  locally 
optimum  detector's  performance  relative  to  the  optimum  when  the  signal  is 
actually  strong. 

For  a  single  sample  (N=l),  all  distributions  of  transformed 
variables  giving  likelihood  ratios  which  are  monotonic  yield  the  same 
detection  performance,  since  then 

Pr  |a ( x ! ) >  nj  =  Pr  |xx  >  nT  (1.1-13) 
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1.1.3  Detector  Design  for  Unknown  Parameters. 

When  either  signal  or  noise  parameters  (or  both)  are  unknown, 
we  can  choose  to  use  estimated  values  of  these  parameters,  risking  a  loss 
in  detection  performance  if  the  actual  values  differ  significantly  from 
those  assumed.  However,  if  the  distributions  of  the  unknown  parameters 
are  known,  decision  theory  for  composite  hypothesis  testing  [l9]  requires 
that  we  use  the  unconditional  likelihood  ration  (LR)  test 

fde1  p(X[01,  Hi)  pfejHj) 

A(X)  =  J7 - -  (1.1-14) 

J de0  p ( X 1 0 0 ,  H0 )  p(e0  |H0) 

where  the  data  X  and  the  parameters  e0  and  ex  can  be  vectors  or  sets. 

If  the  distrubutions  for  the  unknown  parameters  are  not  avail¬ 
able,  or  if  they  are  considered  "unknown  nonrandom",  the  best  test  procedure 
is  not  clearly  specified  by  decision  theory.  However,  since  the  optimum 
performance  would  be  achieved  if  somehow  a  perfect  measurement  were  made 
of  the  unknown  parameters,  it  is  reasonable  to  use  the  generalized 
likelihood  ration  (GLR) 

max  p  (X|e1  ,  Hx) 

A  g(X)  =  - -  (1.1-15) 

max  p  ( X j e 0  ,  H0) 
eo 

For  example,  testing  the  hypothesis  H0 :  x1-=G(0,  a2)  against 
Hj :  x.=G(m,  a2),  where  neither  the  mean  m  nor  the  variance  a2  is  known, 
results  in  the  test  [20] 
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N 


I 


X. 


1 


i=l 


N 


z 


2 


X  . 


1 


n(7)2 


N 


> 

<  n  > 


Z  (xrx)2 

i=l 


(1.1-16) 


which  is  equivalent  to  comparing  the  ratio  of  estimates  (m)2/a2  to  a 
threshold.  (If  m  is  the  value  of  a  constant  "signal",  then  the  ratio  is 
a  form  of  estimated  signal -to-noise  ratio.) 


1.1.4  Intermittent  Gaussian  Interference  as  Non-Gaussian  Noise 


An  understanding  of  non-Gaussian  noise  as  arising  from  inter¬ 
mittent  or  time-varying  noise  is  evident  in  the  present  communications 
literature,  which  reflects  much  concern  over  the  disruption  of  communica¬ 
tions  and/or  detections  due  to  intentional  interference  or  jamming.  In 
frequency  hopping  communications  systems,  for  example,  it  has  been  shown 
that  for  limited  jammer  power  an  effective  jamming  strategy  is  one  for 
which  the  jamming  is  present  in  the  dehopped  bandwidth  some  fraction  of 
the  time  (y),  rather  than  continuously.  Thus  the  marginal  pdf  of  the 
received  noise  at  a  given  time,  assuming  Gaussian  noise  jamming,  is 
given  by 

/2iF  p(x)  =  (1-y)  —  exp(-x2/2a2) 

aN 


V 


+  i 


background  noise  power, 


;  expj-x2/2(a2  +  a2)] 

=  jamming  noise  power) 


(1.1-17) 
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which  has  the  same  form  as  (1.1-4).  However,  this  Gaussian-Gaussian  mixture 
type  of  non-Gaussian  noise  in  these  jamming  situations  is  the  result  of  the 
nonstationary  or  time-varying  properties  of  the  total  noise,  which  is 


Gaussian  at  a  given  time.  That  is. 


p(xt)  =• 


/2tT 


exp 


-x^/2 


(1.1 -18a ) 


where 


with  probability  1-y 
with  probability  y. 


(1-1. 18b) 


Thus  we  observe  that  in  this  case  a  correspondence  or  analogy 
exists  between  "non-Gaussian"  noise  and  "nonstationary  Gaussian"  noise.  For 
example,  the  joint  pdf  of  L  independent  samples  of  noise  from  the  distribution 
given  by  (1.1-4)  is,  using  a=(l-y)/a1v/^nr  and  b=y/a2'/2ri:  » 


PL(x) 


exp< 


exp<  * 


PG(xk;  cti  )  +  y  PGUk;  oi) 


(1-y) 


L-2. 


(2«i 


L-  £ 

I 

kx-l 


* 


(2ira§  ) 


1/2 


(1.1-19) 
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This  may  be  written 

L 

P|_(x)  =  p2,  (1.1-20) 

£=0 

where  may  be  interpreted  as  the  probability  that  "a  of  the  noise  samples 
have  variance  a|  and  l-i  have  variance  cr^",  and  p^(xj£)  is  the  joint  density 
of  the  samples  conditioned  on  this  event. 

With  this  viewpoint  we  can  interpret  Middleton's  Class  A  pdf 

(1.1-3)  as 

p(xt)  =  PG(xt ;  ot)  where  o£  =  with  probability 

p  =  e"^  Am/m!,  m  =  1,  2,  . ..;  (1.1-21) 

m 

that  is,  conditionally  the  noise  is  Gaussian,  with  the  variance  selected 
randomly  from  a  discrete  set  of  values. 
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1.2  DETECTION  THEORY  FOR  MIXTURE  NOISE 

1.2.1  PDF  Formulations  and  Examples 

We  consider  the  situation  in  which  the  noise  is  modelled  as 
having  the  mixture  probability  density  function  (pdf) 


pn(x>  -Y.  =1 


(1.2-1) 


m 


m 


where  each  function  f  (x)  is  a  pdf  weighted  by  a  positive  constant 


m 


m 


0  <  <  1 .  The  argument  x  may  be  considered  a  scalar  for  lowpass  noise  or 
complex  (two-dimensional)  for  bandpass  noise.  For  example,  the  model  pro¬ 
posed  by  Middleton  [3]  has  the  form  in  which  each  pdf  fm(x)  is  Gaussian: 

-x2/2 a2 


»n(x)  -  Y. 


-X.m 
e _ X_ 

m! 


fG(x;  »„).  fG(x;c)  ■ 


m=0 


r/2iT 


(1.2-2) 


Often  two  terms  are  used  in  the  mixture  model,  for  example,  the 
Gaussian-Gaussian  mixture 

pn(x)  =  (1-e)  f g ( x ; a x )  +  e  fG(x;a2),  a2>al,  (1.2-3) 

or  the  Gaussian-Laplace  mixture 

pn(x)  =  (l-e)fG(x;a1)  +  e  •  |  e'alXL  (1.2-4) 

In  that  the  two-term  models  express  the  condition  that  the  noise  is  nearly 
Gaussian  (the  first  term)  but  with  distributional  tails  higher  than  the 
Gaussian  (contributed  by  the  second  term),  they  are  sometimes  called 
"contamination"  models. 
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1.2.2  Interpretation  of  the  Mixture  pdf 

Consider  the  class  of  mixture  densities  in  which  the  component 
pdf's  f  (x)  take  the  same  functional  form,  but  with  different  parameter 
values,  that  is, 

f  (x)  =  f ,(x;c_) , 
m  '  A  m 

where  f^(  ;  )  is  the  common  functional  form  for  the  family  "A"  of  the  pdf's, 
and  5  is  a  particular  value  of  a  parameter  (or  set  of  parameters)  determining 
the  scaling,  location,  etc.,  of  the  pdf.  Then  the  mixture  pdf  has  the 
interpretation 


pn(x)  ‘m  fA(X;  «■> 


=  J  < i{  Pc<5).  (1.2-5) 

Under  this  view,  Pn( x ) is  the  result  of  averaging  the  parametric 
pdf  f^(x;c)  over  a  discrete  pdf  for  the  values  of  5: 


p{<«>  =  Zv1'  ( 

m  m 

For  example,  in  Middleton's  non-Gaussian  noise  model,  the  parameter  5  is 
the  variance  of  a  zero-mean  Gaussian  pdf,  and  the  Km  describe  a  Poisson 
(discrete)  pdf  for  occurrence  of  the  variance  values. 

1.2.3-  Form  of  Likelihood  Ratio 

For  additive  and  independent  signals  and  noise,  the  pdf  under 
the  alternative  hypothesis  under  the  mixture  model  of  noise  becomes 

Wx)  =  pn(x-s)  =  ZKm  Vx-S)- 

m 


17 


J.  S.  LEE  ASSOCIATES,  INC. 


The  likelihood  ratio  (LR)  then  takes  the  form 

V  K  f,  (x-s) 
m  nr  ' 


,  ,  &  ^s+n^  _  m 

A(x>  ‘  TTxT 


f  (x) 

m  nr  ' 

V"  K  f  (x)  [f(x-s)/f m(x)] 
/  .  m  mv  '  1  m  '  m 


m 


y  k  f_(x) 

m  nr 


m 


ZWm(x)  Am(x)-  I“m(x)  =  *• 


(1.2-8) 


m 


m 


a  mixture  of  individual  LR’s  Am(x)  weighted  by  the  (nonconstant)  functions 


K  f  (X) 

W  (x)  =  m  m 
m 


y ■<  f  (x) 

/ .  m  nr  ' 


(1.2-9) 


m 


For  example,  for  the  Gaussian-Gaussian  mixture  (1.2-3), 

/,  X  -!  -x2/2ct? 
(l-e)axie  1 


Wl(x)  = 


n  (  -i  -x2/2a2  ,  -i  -x2/2a? 
(l-elaj/e  1  +  e  a2  e  2 


1  -  W,(x).  (1-2-10) 


In  the  case  that  the  signal  has  unknown  parameters  (e } ,  the  LR 
formulation  (1.2-8)  becomes  the  generalized  likelihood  ratio 


A(x) 


- ■  w 


=  I  Vx>  Ee{Am^x;0)}  • 


(1.2-11) 


m 
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1.2.4  Superposition  of  Detection  Measures 

It  has  been  shown  that  the  pdf  for  mixture  noise  is  an  additive 
combination  of  individual  pdf's: 


p  (*) =  f  (x), 

Knv  '  L  m  nr  '  m 


(1.2-12) 


m 


For  this  reason  the  probability  of  false  alarm  (PpA)  or  Type 
I  detection  error  is,  given  the  threshold  n> 

pFA(n)  =  Pr{A(x)  >  n | H0 } 


=  Pr{x  eR  |H0> 
n 


--J  dx  pn(x) 


dx 


f  (x) 
m 


”  PFA,m(n)‘ 

m 


(1.2-13) 


We  see  that  PpA(n)  is  the  superposition  of  PFA's  arising  from  the  individual 
terms  in  the  pdf.  Similarly,  the  detection  probability  (Pp)  is  a  super¬ 
position  of  individual  probabilities: 


19 


J.  S.  LEE  ASSOCIATES,  INC. 


PD(n;s)  =  Pr{A(x)  >  n | Hi > 


=  Pr{X  £  R  I  Hi } 
n 1 


Z<  f  dx  f(x-s) 

ml  nr 

m  R 

n 


P 


D,m 


( n ;  s ) 


m 


(1.2-14) 


Thus  in  a  significant  manner  the  use  of  mixture  densities  to  characterize 
non-Gaussian  noise  facilitates  calculation  of  detection  measures. 
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1.3 


TREATMENT  OF  BANDPASS  NON-GAUSSIAN  NOISE 


For  zero-mean,  stationary  bandpass  Gaussian  noise,  it  has  long 
been  understood  that  the  sample  function  of  the  noise  random  process  can  be 
represented  in  quadrature  form  by  the  decomposition 

n(t)  =  nc(t)  cosa)0t  -  ns(t)  sinujyt, 

where  f0  =  cj0/2ir  is  the  center  frequency  of  the  band,  nc(t)  is  an  "in-phase" 
random  process,  and  ns(t)  is  a  "quadrature"  random  process.  Both  nc(t)  and 
ns ( t )  are  lowpass  Gaussian  random  processes,  with 


(1.3-2) 


The  correlation  functions  pertaining  to  the  quadrature 
components  nc(t)  and  ns(t)  are 


|nc(t)nc(t+ t)J  =  Ejn$(t)ns(t  + x)J 


(1.3-3) 


and 


B 


-B 


(1.3-4) 
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From  these  correlation  functions  we  observe  that  at  the  same  time  instant 
(t=0),  Rcs=0  and  therefore  nc  and  ng  are  uncorrelated.  Further,  if  the 
noise  power  spectral  density  S  (f)  is  even  about  the  center  frequency,  f0, 
then  Rcs(t)  =  0  for  all  x,  implying  that  nc  and  n$  are  uncorrelated  for  all 
pairs  of  time  instants. 

For  nc(t)  and  ns(t)  Gaussian,  zero  correlation  is  equivalent 
to  statistical  independence,  and  we  can  write  their  joint  pdf  as 


^nc,ns(a’^  ^nc^^ns^ 


2^2-  exp 


n2  +  n2 
c  s 

- 2^" 


(1.3-5) 


Now  the  question  we  wish  to  address  is  how  to  model  the  joint 
pdf  of  the  quadrature  components  of  a  bandpass  non-Gaussian  process.  Under 
the  assumption  that  the  noise  spectrum  is  even  about  the  center  frequency, 
we  can  state  that  nc  and  ns  are  uncorrelated.  If  they  are  not  Gaussian,  the 
lack  of  correlation  no  longer  implies  independence,  although  of  course 
independence  implies  they  are  uncorrelated. 


1.3.1  Independent  Quadrature  Component  Assumption 

One  assumption  which  can  be  made  concerning  the  quadrature 
components  of  bandpass  non-Gaussian  noise  is  that  they  are  independent.  Thus, 
for  example,  for  Gaussian-Gaussian  mixture  noise  we  could  write 


22 


J.  S.  LEE  ASSOCIATES,  INC. 


W'0'6’  =  pnc(a)pns(e) 


(1"e)  c\  Pg(“1)  +  E  o2  ps(°2) 
(1-£)  ^  pg(^t)  +  c^pg(^ 

a2  +  6Z 


11 -e): 

2™f 


"2oT 


+  2^0 x ©2  exp 


2^  " 


+  exp 

2iTa^cr2 


a2  e2 

A  e2 

,  +  32 1 

rol 

Q 

n>n>] 

.  2ttoJ  exp 

‘  2a2  ] 

(1.3-6) 


Such  a  model  was  used  by  Trunk  [20]  to  study  the  performance  of  radar 
detectors  in  sea  clutter. 

For  noncoherent  detection,  we  are  interested  in  the  distribution 
of  the  envelope  of  the  noise.  Using  the  independent  quadrature  version  of 
Gaussian-Gaussian  noise,  the  pdf  of  the  envelope  is 


penv(u) 


2ir 

J  0 


d<$>  p  (u  cos4> ,  u  sini>) 
nc,ns 


=  (1-E>2  •  TT  pr(tt)  ^(37) 

+  2e( 1 — £ )  u 


a  i°2 


exp 


r a 2  +  a2\ 

1  2 


"°1°2 


0 1°2 


(1.3-7a) 


where  pR(*)  is  the  Rayleigh  pdf: 


pR(x)  =  x  e~x  //2,  x  5  0. 


(1.3-7b) 
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It  can  be  observed  from  the  above  expression  that  the  assumption 
of  independent  quadrature  components  yields  an  envelope  pdf  with  some 
complexity.  With  a  signal  present,  the  envelope  pdf  is  very  difficult  to 
analyze,  not  having  a  closed  form  [20]. 

Another  aspect  of  this  noise  modelling  assumption  is  revealed  by 
the  pdf  of  the  phase  of  the  noise,  found  to  be 

oo 

Pphase'*'  *  J iu  u  pnc,ns(ucos*’  us1n*) 


(1-e)2  +  e2 

2tt 


e( 1-e)  +  m  o  o\ 

+  2ir  CTiCT2 +  ^CTi  ”  a2 ^ Z  n2(j)C0S2tf)  (1.3  8) 

Obviously,  the  phase  of  the  bandpass  Gaussian-Gaussian  mixture  noise  is  not 
uniformly  distributed  when  it  is  assumed  that  the  noise  quadrature  components 
are  independent,  nor  is  the  phase  independent  of  the  envelope. 

1.3.2  Circularly  Symmetric  Quadrature  Assumption 

If  it  is  understood  that  the  quadrature  components  nc  and  ng  of 
the  bandpass  noise  are  in  general  statistically  dependent  when  the  noise  is 
nonGaussian,  then  finding  the  form  of  the  joint  density  of  nQ  and  n$  becomes 
the  problem.  Two  methods  [21]  are  available  for  constructing  a  joint  pdf 
when  it  is  assumed  that  the  joint  pdf  has  the  form 


pnc.ns(“-6)  =  9  (^2+S2). 


(1.3-9) 


some  function  of  /a2  +  32  that  is,  the  pdf  possesses  "circular  symmetry". 

9 
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The  first  method  is  to  select  the  function  g(R),  recognizing  that 
R  is  the  envelope  of  the  noise,  and  assuming  that  the  phase  of  noise  relative 
to  the  center  frequency  is  independent  of  the  envelope.  Then  the  joint  pdf 
of  envelope  and  phase  becomes 


R,e 


(u,e)  =  K^ 


(1.3-10) 


where  K  is  a  normalization  constant.  This  approach  has  been  used  to  analyze 
the  performance  of  communications  in  VLF  atmospheric  noise  [22],  with  the 
envelope  assumed  to  have  a  log  normal  distribution.  The  marginal  distributions 
of  the  individual  quadrature  components,  Pnc(<*)  and  Pns(s)  follow  from  the 
joint  pdf  and  we  do  not  have  control  over  their  form  using  this  method. 

The  second  method  for  constructing  a  joint  pdf  is  based  on 
generalizing  the  marginal  pdf's  of  the  quadrature  components.  Let  the 
characteristic  function  of  the  individual  quadrature  components  be 

jvn 


(  JvrV 

C(v)  =  E  e 


=  E 


(1.3-11) 


A  circularly  symmetric  joint  distribution  can  be  assigned  to  the  quadrature 
components  by  defining  the  joint  characteristic  function  to  be 


cnCinS(v,u}  =  C  (Ai2Tii2  , 


(1.3-12) 


For  example,  if  the  quadrature  components  are  Gaussian-Gaussian  mixtures, 
then 


c(v)  =  (l-e)e  +  £  e 


(1.3-13) 


and  the  circularly  symmetric  joint  characteristic  function  is 


Cnc,ns(v>u)  =  (1-e)  e 


(v2  +  vi2)  (1.3-14) 


+  e  e 
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From  this  assumption  it  follows  that 


"  (1_^-  CXD 

1 

i  +  B2- 

,  £ 

\ 

a 2  +  | 

2no2  eXP 

"  2a2 

2ira2  exp 

2a  ^  l 

1 

1 

z 

2 

(1.3-15) 


Throughout  our  subsequent  analysis,  we  shall  use  the  circularly  symmetric, 
rather  than  independent,  quadrature  assumption. 
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2  0  NONCOHERENT  DETECTION  IN  BANDPASS  GAUSSIAN-GAUSSIAN 

MIXTURE  NOISE  USING  CONVENTIONAL  DETECTORS 

In  this  section  we  examine  the  effects  on  detection  performance 
that  occur  as  the  noise  departs  from  a  Gaussian  distribution.  The  detectors 
will  be  based  on  the  assumption  that  the  noise  is  Gaussian. 

As  we  develop  these  results,  several  purposes  are  in  mind.  First, 
the  performance  of  Gaussian  detectors  in  the  non-Gaussian  noise  environment 
will  serve  as  a  useful  reference  or  yardstick  for  evaluation  of  the  per¬ 
formance  of  improved  detectors.  Second,  the  results  will  reveal  whether 
the  use  of  multiple  samples  will  overcome  the  loss  in  single-sample  detector 
performance,  so  that  acceptable  performance  is  achieved  in  spite  of  the 
fact  that  the  noise  is  not  Gaussian.  Third,  by  comparing  the  performance 
of  single-channel  Gaussian  detectors  in  non-Gaussian  noise  with  that  of 
detectors  utilizing  two  channels  of  data  (specifically,  a  correlator- 
detector),  we  will  learn  whether  one  type  of  these  common  detectors  is  less 
vulnerable  to  the  degradation  from  non-Gaussian  noise  than  the  other. 

2.1  DETECTION  FORMULATIONS 

As  a  particular  case  of  detection  of  signals  in  non-Gaussian 
noise,  we  turn  to  the  problem  defined  as  follows:  on  the  basis  of  the 
received  waveform  r(t),  0<t$T,  we  wish  to  accept  or  reject  the  null 
hypothesi s 

H0:  r(t)  =  n(t)  (noise  only) 

when  the  noise  is  bandpass  Gaussian-Gaussian  mixture  noise, 

n(t)  =  nc(t)cosiid0t  -  ns(t)sinw0t 


(2.1-1) 


(2.1-2) 
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where  f0  =  aj0/2ir  is  the  center  frequency  of  the  band,  and  the  joint  pdf  of 
the  quadrature  components  nf(t),  ng(t)  at  a  given  instant  is  the  bivariate 
Gaussian-Gaussian  mixture  pdf 


nc,ns 


1-e 

a2  +  B2 

-1-  — — — -  oYn 

a2  +  32 

p 

-  2a? 

2iro|  XP 

Zal 

(2.1-3) 


Detection  occurs  when  Ho  is  rejected  in  favor  of  the  alternative 

hypothesis 


Hx :  r(t)  =  n(t)  +  Acos[cj0t  +  e(t)],  (2.1- 

in  which  the  signal  amplitude  A  is  constant  during  the  observation  interval* 
and  the  signal  phase  e(t)  is  random.  Two  different  assumptions  will  be  made 
about  the  phase:  (a)  the  random  phase  is  constant  ("slowly  varying")  during 
the  observation  interval  (type  1  signal);  or  (b)  the  random  phases  of  samples 
taken  during  the  observation  interval  are  independent  (type  2  signal). 

We  assume  that  K  samples  of  rc  and  r$,  the  quadrature  components 
of  r(t) ,  are  taken  on  the  interval  (0,T).  Under  the  two  hypotheses,  the 
joint  pdf's  of  these  samples  are 


H0 : 

Prc,rs(2-’  =  pnc,ns(s-’^ 

(2. l-5a) 

Hi: 

Prc.rs'i-  il«l’  =  pnc,ns(i- V  £'*s> 

(2. l-5b) 

where  the  vector  notation  signifies 

*  In  certain  cases  shown  below,  the  extension  of  the  analysis  to  amplitude 
variations  is  possible  in  a  simple  manner.  The  assumption  of  constant 
amplitude  signals  is  commonly  made  as  a  means  to  simplifying  the  analysis, 
realizing  that  in  practice  variation  is  almost  always  observed. 
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Lc  =  [rc(ti),  rc(t2),...,rc(tK)] 
Is  =  ['"s(t1),  rs(t2),...  ,r$(tK)] 
^  =  [sc(t1 ) ,  sc(t2),... ,sc(tK)] 


=  Acose  [1,  (Type  1) 

=  Afcose!,  cos02 , . . .  ,cos0^,]  (Type  2) 

=  [ss(t1),  ss(t2),...,ss(tK)] 

=  Asin0  [1,  1, . . .1]  (Type  1) 

=  A[sin8i,  sin02,...,  si ne^]  (Type  2). 


(2.1-5c) 


The  test  for  rejecting  H0  in  favor  of  Hi 
generalized  likelihood  ratio  (GLR) 


is  to  be  based  on  the 


Ar(“>i) 


Mprc,rs(2-’£lHi>4 

prc,rs(-’^lHo) 


(2.1-6) 


We  also  shall  consider  the  extension  of  this  formula  to  the 
situation  in  which  two  channels  of  data  (from  perhaps  two  sensors)  are  to 
be  tested  for  the  presence  of  the  same  signal.  In  this  case,  the  GLR  becomes 

Ar  r  («i»  3 1»  62)  =  il)  ^(“2*  £2)*  (2.1-7) 

» 

assuming  that  the  noises  in  the  two  channels  are  independent,  and  expanding 
the  notation  of  (2.1-5)  and  (2.1-6)  in  an  obvious  way. 
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2.2  PERFORMANCE  OF  SINGLE  CHANNEL  GAUSSIAN  DETECTORS  IN  GAUSSIAN- 

GAUSSIAN  MIXTURE  NOISE 


When  the  mixture  parameter  e  in  the  noise  pdf  (2.1-3)  is  zero, 
the  likelihood  ratio  (2.1-6)  becomes 


K 

Ar  (a,  g)  =  E  n  exp 
L  ^  k=l 


(ak  -  A  cos9k)  +  (3k  -  A  sineK) ^  ak2+gk2 


2ai2 


2ai5 


exp 


KA 

2^ 


|expj-j-  ^  (ak  cos  ek  +  sk  sin  ek^| 


(2.2-1) 


assuming  the  samples  (ak,  $k)  are  independent. 

2.2.1  Forms  of  the  Detectors 

For  the  Type  I  signal,  the  phases  {ek}  are  all  equal,  and  the 
result  of  the  expectation  taken  in  (2.2-1)  is 

Ar(a,  §)  =  exp  j-  f^2  j  Iq  j  i)  }  (2.2-2) 


where  I0(-)  is  the  modified  Bessel  function  of  the  first  kind  and  order 
zero,  and 

f(».  i)  -  ^  +  «k)  (2-2-3) 

Since  the  likelihood  ratio  is  monotonic  or  directly  proportional  to 
f(a,  _§_) ,  testing  the  likelihood  ratio  is  equivalent  to  testing  f(a,  B. ) , 
that  is, 
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Ar  (“>  A)  <  Hi 


Hi 

> 


f(a,  _§_)  <  n2 

H  n 


(2.2-4) 


The  detector  based  on  (2.2-4)  is  diagrammed  in  Figure  2.2-l(a). 
For  the  Type  II  signal,  the  phases  { e^}  are  all  independent,  and 
the  result  of  the  expectation  taken  in  (2.2-1)  is 

KA2  i  K 


A  (a,  b)  =  exp  ]- 
r  1  2  a? 

KA2 


IJ:  A2  +  ) 

kij  £ni«(4“V“k2 + b3|  .  <2'2'5> 


=  exp  —  + 

I  2a2  k= 


When  the  signal  is  weak,  we  can  simplify  (2.2-5)  greatly  by  using  the 
approximation 


In  I, 


( 

\  of 


“k +  sk 


A  H + 


(2.2-6) 


4of 


to  arrive  at  the  equivalent  detection  test 

K 

I  (“I  +  6k  >  <  "• 

k=l  H0 


(2.2-7) 


This  practical  implementation  of  the  test  of  (2.2-7)  is  diagrammed  in 
Figure  2.2-l(b). 

2.2.2  Single-Sample  Detector  Performance  in  Gaussian  Mixture  Noise 
For  one  sample  (K=l),  the  two  detectors  implement  the  test 


r2  +  r2 


Hi 

R2  <  n. 
Ho 


(2.2-8) 
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COSOjQt 


Hi 

H0 


(a)  Type  1  signal  (constant,  random  phase) 


R  k  =  r  ck  +  r  sk 

(b)  Type  2  signal  (Independent,  phase  samples) 


Figure  2.2-1  Detectors  for  a  bandpass  signal  in  Gaussian  noise. 
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where  R  is  the  envelope  of  the  received  waveform  at  the  sample  time.  Thus 
instead  of  quadrature  sampling,  we  may  employ  an  envelope  or  square-law 
envelope  detector.  It  is  well-known  that  in  Gaussian  noise ,  R2  is  a  constant 
times  a  noncentral  chi -squared  random  variable  with  two  degrees  of  freedom 
and  noncentrality  parameter 


x  = 


/  2p ,  Hi  true  ,  p=  A/202; 


0  ,  H0  true. 


(2.2-9) 


The  performance  of  the  quadrature  detector  for  Gaussian-Gaussian 
noise  in  terms  of  false  alarm  and  detection  probabilities  is  found  to  be 


Pp^(n)  =  Pr |r2  >  n |H0j 


-n/2o2  -n/2o| 

=  (1-e)  e  +  e  e 


(2.2-10) 


and 


=  (1-e)  Q(A/ai,  Ai/cr2  )  +  £  q(a/o2,  frf /o|  j  , 

(2.2-11) 

where  n  is  the  detection  threshold  and  Q(a,b)  is  Marcum's  Q-function. 

The  effect  of  the  mixture  parameters  e  and  V2  -  o2/o2  on  the  false 
alarm  threshold  for  the  square-law  envelope  detector  is  shown  in  Table  2.2-1. 
It  is  evident  from  this  table  that  for  Pp^  small  (<0.1),  the  tendency  is  for 
the  first  term  in  (2.2-10)  to  be  negligible,  resulting  in 


-2v2*n  (PFA/‘). 


(2.2-12) 
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TABLE  2.2-1 

FALSE  ALARM  THRESHOLDS 

FOR  SQUARE-LAW  ENVELOPE  DETECTOR  (SINGLE  SAMPLE) 


PFA 

Gaussian 

Noise 

n  h\ 

GAUSSIAN- 

■GAUSSIAN 

NOISE 

V2  =  10 
n/of 

£  =  0.1 

V2  =  100 
n/of 

V2=  10 
r\/o\ 

£  =  .01 

V2  =  100 

n/crf 

0.1 

4.60517 

6.8662 

10.3657 

4.7494 

4.7905 

0.01 

9.21034 

46.0517 

460.5170 

10.9215 

14.5094 

0.001 

13.8151 

92.1034 

921.0340 

46.0517 

460.5170 
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and  the  threshold  is  raised  to  a  value  much  higher  than  that  for  Gaussian 
noise.  Therefore,  for  the  same  false  alarm  probability,  a  higher  SNR  will 
be  required  to  achieve  the  same  detection  probability.  Note  that  if  either 
e=0  or  V2  =  1,  then  the  probabilities  of  false  alarm  and  detection  become 
those  for  the  Gaussian  noise  case. 

In  Figures  2.2-2  to  2.2-5,  the  false  alarm  probability,  as  a 
function  of  the  normalized  threshold  n h\,  is  plotted  for  different  values 
of  the  mixture  parameter,  e,  and  the  variance  ratio,  V2: 


e  =  (.01,  .1,  .2,  .5) 

V2=  (1,  2,  10,  100,  1000). 


(2.2-13) 


It  is  evident  from  these  figures  that  for  PFA  1  ess  than  or  equal  to  the 
mixture  parameter  e,  the  threshold  is  determined  by  the  second  term  in 
(2.2-10),  as  noted  already  in  (2.2-12). 

In  Figures  2.2-6  to  2.2-9,  the  detection  probability  is  plotted  as 
a  function  of  SNR  for  the  same  parametric  conditions  as  described  by 
(2.2-13),  and  for  fixed  false  alarm  probabilities  of  10  *,  10  2,  and  10  3. 

As  anticipated;  for  each  value  of  e,  the  detection  probability  decreases 
as  V2,  the  variance  ratio,  increases,  except  in  some  cases  for  high  SNR. 

We  observe  also  that  the  degradation  in  performance  for  Pp>. 5  is  proportional 
to  e,  and  that  in  general  the  amount  of  degradation  is  greater  for  smaller 
values  of  Pp^,  the  false  alarm  probability.  This  result  is  consistent 
with  the  fact  that  the  tails  of  the  distribution  are  extended  by  the 
contaminating  Gaussian  noise  with  variance  o§.  Even  with  e  as  small  as 
0.01  (Figure  2.2-6),  a  loss  in  detectability  of  over  4  db  is  experienced 
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Gaussian  mixture  noise,  mixture  parameter  e  =  0.01. 
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o 


Gaussian  mixture  noise,  mixture  parameter  e  =  0.1. 
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o 


Gaussian  mixture  noise,  mixture  parameter  e=0.2. 
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o 


Gaussian  mixture  noise,  mixture  parameter  e=0.5. 
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Figure  2.2-7  Receiver  operating  characteristics  for  Gaussian  detector 
in  Gaussian-Gaussian  mixture  noise,  mixture  parameter 
e  =  0.1 


41 


DETECTION  PROBABILITY 


42 


DETECTION  PROBABILITY 
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SIGNAL  TO  NOISE  RATIO  CdB) 

Figure  2.2-9  Receiver  operating  characteristics  for  Gaussian  detector 
in  Guassian-Gaussian  mixture  noise,  mixture  parameter 
e  =  0.5. 
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for  V2=10;  this  is  the  increase  in  SNR  required  to  maintain,  say,  90% 
detection  probability  as  V2  goes  from  the  value  1  to  10. 

For  smaller  values  of  PQ  and  V2  >  10,  we  observe  in  Figures 
2.2-6  to  2.2-9  that  the  Gaussian  detector  performance  actually  improves 
slightly  as  e  increases. 

As  we  consider  next  the  performance  of  the  Gaussian  detectors 
of  Figure  2.2-1  in  Gaussian-Gaussian  mixture  noise  for  multiple  samples, 
we  shall  be  interested  to  learn  whether  detection  losses  can  be  compensated 
by  using  multiple  samples. 


2.2.3  Multiple-Sample  Detector  Performance  In  Gaussian  Mixture  Noise 


The  joint  pdf  of  K  multiple,  independent  samples  of  the 
quadrature  components  of  the  input  waveform  is  simply  a  K-fold  product 
of  the  single  sample  pdf  (2.1-3).  Because  this  pdf  has  a  two-term  or 
binomial  form,  the  K-sample  joint  pdf  for  the  noise  only  case  can  be 
wri tten 


K  I  _ 

pr  r  £.)  =  n  — ^exp 

- C ,  — S  k=l  2ircr^ 


ok2  + 


6k2} 


2a? 


r 


exp  « 


V+A]  1 

2a2  >  J 


(a?rK+m 


-m 


(2.2-14) 

This  expression  conveys  the  information  that  the  joint  pdf  consists  of 
K+l  terms,  each  of  which  is  a  (weighted)  joint  pdf  for  K  independent 
Gaussian  random  variables,  m  of  which  have  variance  a2  and  K-m  of  which 
have  variance  o2,  provided  that  the  indexing  or  time-ordering  of  the 
samples  ( rc|< »  r^)  is  arbitrary. 
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Another  method  for  writing  the  multiple  sample  pdf  is  to 
consider  the  samples  as  conditionally  Gaussian: 


Pr  r 
— c  5 — s 


s) 


(2.2-15) 


where  the  variance  multipliers  |v£j  can  take  the  values  1  or  Q\l°\  -  V2, 
that  is, 


Pv2  (y)  =  (1-e)  <5  (t-  1 )  +  e6(y-V2).  (2.2-16) 

vk 

2.2.3. 1  Sum  and  Square  Detector 

First  we  consider  the  "sum  and  square"  quadrature  detector  of 
Figure  2.2-l(a),  the  GLR  for  a  Type  1  (constant  phase)  signal  in  bandpass 
Gaussian  noise.  The  test  statistic  for  this  detector  may  be  written  as 


z 


(2.2-17) 


Since,  given  the  variance  multipliers  jv£j,  The  quadrature  samples  are 
independent  Gaussian  random  variables,  so  are  their  sums;  that  is, 
conditionally, 

K  K 

Z  rck  =  Z  G(A  cos6-  VH) 

k=l  k=l 


=  G  KA  cose 


•Z 


vk°i 


k=l 


(2.2-18a) 


sk 


K  /  K  \ 

G (a  sine,  V2a2)  =  g|kA  sine,  v^a2j 
k=l  V  k=l  / 


(2.2- 18b ) 
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Therefore  conditionally  z  is  a  factor  times  a  chi-squared  random  variable 
with  two  degrees  of  freedom: 


with  the  noncentrality  parameter 


X  = 


* 


(2.2- 19b ) 


using  p  =  A2/2af.  Although  the  detector  form  is  based  on  constant  phase 
and  "slowly-varying"  signal  amplitude  A,  the  performance  of  the  detector 
can  be  evaluated  for  time  varying  amplitude  and  phase.  The  modification 
necessary  for  this  evaluation  is  to  interpret  the  noncentrality  parameter 
X  in  (2.2-19b)  as 


(2.2-19c) 


(2.2-19d) 


Conceivably  the  physical  process  giving  rise  to  the  non-Gaussian 
noise  can  introduce  dependence  among  the  variance  multipliers  jv^j.  For 
example,  at  one  extreme,  slow  variation  or  low  bandwith  in  factors 
affecting  the  power  of  a  conditionally  Gaussian  process  may  allow  us  to 
consider  that  vj^v2,  k=2,  3,  ...,  K;  that  is,  the  value  of  v^  is  a  random 
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constant.  In  this  instance  it  can  be  shown  that  the  detection  and  false 
alarm  probabilities  achieved  by  the  sum  and  square  detector  are 


PD(n;  p,  K)  =  Ey2 [Pr  {z>n|H1}] 


( 1  - e )  Q 2 Kp  ,  Vn/1<cri  J 

+  eQA/2Kp/V2,  Vn/KafV2  \ 


where  Q(a,  b)  is  Marcum's  Q-function,  and 


PpA  (n;  K)  =  Pn  (n;  0,  K) 


D 


-  (1-c) 


v^=vf ;  (2.2-20) 


exp  |  -  n/2Kof|  +e  exp  j-  n/2Kaf  V‘ 


(2.2-21) 


Another  extreme  case  of  the  statistical  relationship  among  the 
jv2j  is  that  in  which  high  bandwidth  in  the  variation  of  the  power  of  a 
conditionally  Gaussian  process  permits  us  to  assume  that  the  jv|j  are 
independent.  Then  the  detection  and  false  alarm  probabilities  are 


v*  -  *  ■  i  (:) 

m=0  x  '  x 


n  ja\ 


m+mV2  ’  y  K-m+mV2 

independent  v^;  (2.2-22) 


and 


PFA  (n;  K) 


K 

Y  (l-e)K"m  £m  exp  {-  n/2af(K-m  +  mV2)} 
m=0  ' 


(2.2-23) 
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Figure  2.2-10.  Receiver  operating  characteristics  for  the  sum-and-square 
Gaussian  detector  in  Gaussian-Gaussian  mixture  noise  <e=0.1,  V2  = 100)  when 
multiple  samples  are  used  and  the  noise  power  is  slowly  varying. 
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The  difference  in  sum-and-square  detector  performance  due  to  the 
two  different  assumptions  about  the  dependence  of  the  variance  multipliers 
can  be  observed  by  comparing  Figures  2.2-10  and  2.2-11,  in  which  e=0.1  and 
V2=100. 

In  Figure  2.2-10,  the  equal  (v  £}  case  is  evaluated,  showing  that 
an  improvement  in  detection  holds  for  this  case.  A  close  look  at  the  equal 
{v2}  Pp  and  Pp^  expressions,  (2.2-20)  and  (2.2-21),  reveals  that 

PD  (n;  p,  K)  =  PD  (Pp^;  KP,  1),  equal  V|<2;  (2.2-24) 

that  is,  the  detectability  of  the  signal  is  increased  by  the  factor  K.  Thus 
the  single-sample  detection  losses  tabulated  previously  can  be  made  up  by 
using  an  appropriate  number  of  samples,  assuming  that  the  variance  is  indeed 
constant  for  the  K  samples.  However,  relative  to  the  performance  achievable 
in  Gaussian  noise  with  the  same  number  of  samples,  the  detectability  loss 
remains  the  same,  regardless  of  the  value  of  K. 

The  sharp  rise  in  the  curves  for  PpA  =  io“2  in  Figure  2.2-10  can 
be  explained  as  follows:  the  false  alarm  probability  (2.2-21)  is  dominated 
by  the  second  term,  yielding 

n/o2  *  -2KV2  £n(PpA/e)  =  460K.  (2.2-25) 

The  detection  probability  (2.2-20)  is  due  to  the  second  term  for  small 
SNR;  when  the  SNR  approaches  the  value  at  which  the  first  term's  Q-function 
becomes  0.5,  or 

2Kp  -  n/Ka2  =  460,  (2.2-26) 
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°  Gaussian  detector  in  Gaussian-Gaussi an  mixture  noise  (e=0.1,  =100) 

Q 

a 

when  independent  multiple  samples  are  used. 
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this  term  rapidly  increases  from  zero  to  (1-e)  -  0.9.  For  example,  if 
K=10  this  occurs  at  p=  23  =  13 . 6dB . 

The  sum-and-square  detector  performance  for  independent  {v|},  ar>d 
hence  independent  samples,  is  plotted  in  Figure  2.2-11.  A  very  different 


behavior  from  that  of  Figure  2.2-10  is  observed:  for  the  case  of  PpA  =  10  1 , 


the  independent  sample  performance  actually  gets  worse  initially  as  K 
increases;  however,  for  the  case  of  PpA=  ,the  detection  probability 
increases  uniformly  with  the  number  of  samples  K,  and  is  better  than  th*-* 
results  shown  in  Figure  2.2-10.  We  conclude  that  for  low  PpA>  the  per¬ 
formance  of  this  detector  is  better  for  independent  noise  samples,  and 
that  effects  of  the  non-Gaussian  noise  can  be  countered  by  using  multiple 
samples. 

2. 2. 3. 2  Square  and  Sum  Detector 


Next,  we  consider  the  "square  and  sum"  quadrature  detector  of 


Figure  2.2-l(b),  the  optimum  detector  for  a  weak  Type  2  signal 
(independent  phase  samples)  in  bandpass  Gaussian  noise.  The  test 
statistic  for  this  detector  may  be  written  as 


(2.2-27) 
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Conditionally,  z  is  the  sum  of  K  weighted  chi-squared  random  variables  with 
two  degrees  of  freedom, 

K 

z  =  Z  Vk  *2(2lXk)’  Xk  =  2p/vk-  (2.2-28a) 

k=l 

The  evaluation  of  this  detector  for  non-constant  signal  amplitude  involves 
modifying  (2.2-28a)  by  considering 

X£  =  2p|</v|,  =  A^/2a^  .  (2.2-28b) 

Since  equally-weighted  chi-squared  random  variables  combine,  we 
see  that  unconditionally  the  probability  of  detection  is 


independent  v£;  (2.2-29b) 

where  CL  (a,  b)  is  the  generalized  Marcum's  Q-function  [24],  and 

P  (n;  p,  K)  =  P r  (afx2[2(K-m),  2(K-m)Pr  +  aTV2x2[2m,  2mP/V2]>n} 

Dm 

=  QK_m  (V2(K-m)P ,  VnM  ) 

/•n/ai  - - \ 

+  /  du  Pl(u)  Qm(V2mp/V2,  V(n-aiu)/a?V2)  ) 

o  (2.2-29c) 
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with 


Pi (u )  =iexpj-|-  (K-m)p 


{ 

r  u/2 

1 

(K-m)p 

(K-m-l)/2 


I( 


y2(K-m)pu 


(2 . 2-29d) 

In  these  expressions,  we  may  treat  non-constant  signal  amplitude 
by  replacing  Kp  by  Kp',  where 


A? 


=  i  y  _l_  =  jl  y  _ i  =  i 

K  y  2of  KAt  y  of  K 


(2.2-29e) 


in  which  is  the  signal  energy  and  N0  is  the  equivalent  noise  spectral 
density  in  the  bandwidth  B  =  1/At. 


The  performance  of  the  square-and-sum  detector  for  equal  (v^l  is 
shown  in  Figure  2.2-12,  for  c=0.1  and  V2  =100,  with  several  values  of 
K  and  Pp^.  Unlike  that  of  the  sum-and-square  detector  (Figure  2.2-10)  for 
the  same  dependent  noise  sample  case,  as  the  number  of  samples  K  increases, 
the  detection  probability  for  this  detector  tends  to  jump  from  a  small 
value  (-Pp^)  to  a  high  value  at  a  fixed  value  of  SNR.  That  is,  increasing 
the  number  of  samples  does  not  improve  the  detectability  in  terms  of 
reducing  the  required  SNR  for  a  given  value  of  P^.  This  can  be  shown 
analytically  as  follows:  for  PpA  <  e  and  large  K,  the  detection 
threshold  approaches 

rT  =  n/af  =  inverse  of  =  PFA/,£| 


2V2/K[xq(PFA/e)  +  A  ], 


(2 . 2-30a) 
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Figure  2.2-12  Receiver  operating  characteristics  for  the  square-and-sum 
Gaussian  detector  in  Gaussian-Gaussian  mixture  noise 
(e  =  0.1,  V2  =  100)  when  multiple  samples  are  used  and 
the  noise  power  is  slowly  varying. 
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where  Q  (•)  is  the  Gaussian  complementary  probability  integral  and  (•) 
is  defined  as  its  inverse: 

Xq  (p):  Q  (Xq)  =  P-  (2.2-30b) 

Given  this  threshold,  the  detection  probability  (2.2-29a)  for 
large  K  is  approximately 


PD  =■  (1-C)  Ql 


JV-2K(l+p)l 

+  e0 

*  n  VV2-2K(1+p/V2) 

[  2/K  /l+2p j 

2/K/1+2p/V2  J 

(2.2-31) 


in  which  the  Q-function  in  the  first  term  is  equal  to  0.5  for  tT-2K(1+p)=0, 
or 

P  =  nV2K  -1 


=  V2 


1  + 


Xn(PFA/e) 

v'K 


-1. 


(2.2-32) 


For'example,  if  PpA=10  and  e=  0.1,  then  Xq=  1.28  and  the  will  jump  at 
approximately  p  =  128  =  21  dB  for  K=20  and  V2=100.  Although  the  approximations 
no  longer  apply,  for  K=1  the  value  is  p  =  227  =  23.6  dB,  which  can  be  verified 
by  examining  Figure  2.2-10;  in  that  figure  the  switching  occurs  at 
p  =  17dB  for  K=5,  which  implies  that  it  occurs  at  p  =  17+7  =  24  dB  for  K=l. 

For  Pp A=e  as  is  the  case  for  the  PpA=10-1curves  in  Figure  2.2-12, 
the  foregoing  analysis  is  inadequate  since  Xq(l)  =  -  °°>  both  terms  in  the 
probability  expression  then  contribute  to  the  PpA-  Nevertheless  from 
Figure  2.2-12  we  observe  that  the  switching  phenomenon  exists,  and  for  this 
case  as  K  increases  the  PQ  is  approximately  0.1=PpA  until  p=6  dB. 
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The  performance  of  the  square-and-sum  detector  for  independent 
Gaussian-Gaussian  mixture  samples  (independent  )  is  shown  in  Figures 

2.2-13(a)  and  (b)  for  e  =  0.1  and  V2  =  100.  PpA  =  10_1  in  Figure  2.2-13(a) 
and  PF.  =  10"2  in  Figure  2.2-13(b),  and  the  numbers  of  samples  for  which 
(2.2-29b)  was  computed  are  K=l,  2,  5,  and  10. 

For  P,-fl  =  10”1,  we  observe  in  Figure  2.2-13(a)  that  a  loss  in 

1  A 

performance  occurs  when  two  samples  are  used,  but  for  smaller  detection 
probability  values  (Pq  =  .5),  this  loss  becomes  less  as  K  is  increased 
further. 


For  Pcfl  =  10  ,  in  Figure  2.2-13(b)  we  see  that  the  detection 

r  A 

performance  for  the  square-and-sum  detector  improves  uniformly  as  the 
number  of  samples  increases. 

In  both  part  (a)  and  (b)  of  Figure  2.2-13,  a  (temporary) 
saturation  or  leveling  of  the  PQ  curve  is  evident.  For  the  case  of 
p  =  10"1,  this  effect  can  be  explained  as  follows.  For  large  K,  the 

rn 

central  limit  theorem  suggests  that 
Pr  {  z  >  =  Q 


n  -  2Ka2  (l-£  +  £ V2  +  p) 


2o2v/K  /1-e  +  eV4  +  2  (1-e  +  eVz)p_ 

1  1 


(2.2-33) 


where  Q(  )  is  the  Gaussian  Q-function.  For  PpA  -  10  1 ,  p-0  in  (2.2-33) 
gives  for  the  threshold 

n /a2  =  (80. 99)/K  +  (21.8)K,  (2.2-34) 

2  t 
when  e  =  0.1  and  V  =  100.  Substituting  this  threshold  in  (2.2-33)  gives 

the  approximate  detection  probability 
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Gaussian  detector  in  Gaussian-Gaussian  mixture  noise  (e=0.1,  V2  =  100) 

Q 

when  independent  multiple  samples  are  used,  and  Pp^  =  10  1. 
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Figure  2 . 2-13 ( b )  Receiver  operating  characteristics  for  the  square  and 

sums  Gaussian  detector  in  Gaussian-Gaussian  mixture  noise  (e=0.1,  V2  =  100) 

_2 

when  independent  multiple  samples  are  used,  and  Pp^  = 10  . 
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p 


D 


Q 


(80.99)/K  -  2Kp 

2/1 K  /1000.9  +  21. 8p 


(2.2-35) 


Thus  we  anticipate  that  Pg  =  0.5  when  p  =  81/2/K  and  observe  that  the  rate 
of  increase  in  the  detection  probability  as  p  increases  will  start  slowing 
down  at  the  "breakpoint"  p  =  1001/22  =  16.6dB.  For  K=10,  both  these  effects 
are  observed  at  or  about  these  values  of  SNR,  even  though  ten  samples  are 
too  few  to  invoke  the  central  limit  theorem. 
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2.3  PERFORMANCE  OF  CORRELATION  DETECTOR  IN  BANDPASS  GAUSSIAN-GAUSSIAN 

MIXTURE  NOISE 

Representative  of  the  class  of  multi-sensor  detectors  is  the 
bandpass  correlator  diagrammed  in  Figure  2.3-1.  Two  bandpass  waveforms 
u1(t)  and  u2(t),  where 

u.(t)  =  A.  cos  (to  t  +  e. )  +  n.(t),  i  =1 ,  2;  (2.3-la) 

1  I  V  .  I  I 

=  u  .(t)  cos  to- 1  +  u  .(t)  sin  to-t,  (2.3-lb) 

Cl  c  si  c 

are  multiplied,  then  lowpass  filtered  to  produce  the  output  y(t).  For 

the  ideal  assumption  of  zonal  lowpass  filtering,  the  output  can  be  expressed 

as 


y(‘)  -  ?Sucl(t>  uc2(t)  +  uSl(t)  us2(t)| ' 


(2.3-2) 


The  distribution  of  y(t)  for  Gaussian  noises  nx(t)  and  n2(t) 
was  shown  in  [28].  Briefly,  for  that  case  the  vector  of  quadrature  components 


u  =  (u  ,  u  ,  u  ,  u  ) 
-  cx  Si  c2  s2 


(2.3-3) 


is  multivariate  Gaussian  with  mean 


m^  =  (A1cos  e1,  Axsin  els  A2cos  e2,  A2sin  e2) 


(2.3-4) 
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and  covariance  matrix 


0  £<70.  ra  CT. 

a  b  a  b 


■re  a.  e,oo 


a  b  a  b 


5Vb 


ra  a. 
a  b 


‘rVb  % 


E,a  a. 

a  b 


(2.3-5) 


In  (2.3-5)  a2  is  the  variance  of  ni(t),  al  is  the  variance  of  n2(t),  e,  is 
the  correlation  coefficient  between  quadrature  components  which  agree  in 
phase  and  ±r  is  the  correlation  coefficient  between  components  which  are  in 
phase  quadrature. 

Given  these  definitions,  in  [28]  is  shown  that  the  correlator 
output  y(t)  is  equivalent  to  the  difference  between  two  independent  scaled, 
noncentral  chi-squared  random  variables  with  two  degrees  of  freedom: 


y  -  klX2(2;  xx)  -  k2x2(2;  X2). 


(2.3-6) 


In  this  equivalence,  we  have 


ki  =  aaab^  /ir^2+  S)/4 


(2.3-7a) 


k2  =  oa<7b(vT=r2  "  5)/4 


(2.3-7b) 
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and  the  noncentrality  parameters  are  given  by 


Pi  +  p2  +  2/L-rz  /p !P2 cos (6 1  -  02)  +  2r ✓p1p2  sinfej.  -  e2  ) 

/l-rz  ( vl-rz  +  ?) 


(2.3-8a) 


Pj  +  p2  -  2/l-rz  /pxp2 cos ( 0 i  -  e2)  +2r/Plp2  sin(e1  -  e2) 

X2  =  -  (2.3-8b) 

✓l-rz  (/l-rz  -  5) 

using  the  channel  SNR's  defined  as 

Pi  “  Af  /  2o| ,  p2=  A|/2az  (2.3-9) 

When  the  bandpass  spectrum  of  n^t)  and  n2(t)  is  symmetric,  the 
cross-quadrature  correlation  coefficient  r  is  zero.  If  we  assume  that  this 
condition  holds,  and  also  that  the  correlator  delay  shown  in  Figure  2.3-1 
is  such  that  ex  =  e2  ,  then  the  parameters  in  (2.3-6)  to  (2.3-8)  become 

kjL  =  aaab(l  +  5)/4,  k2  =  aaab(l  -  $)/4  (2.3-10a) 

Xi  =  (^pT  +  /pl)2/(l  +  z)  (2.3-10b) 

2 

X2  =  (l/pY  -  /p^)  /(I  -  5) 

2.3.1  Analysis  for  Gaussian  mixture  noise 

The  assumption  of  Gaussian-Gaussian  mixture  noise  for  n-^t)  and 
n2(t)  can  be  treated  by  considering  the  covariance  matrix  (2.3-5)  and  the  SNR's 
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p!  and  p2  (whose  values  follow  from  the  covariance  matrix)  to  be  random 
variables.  That  is, 

K  =  K  with  probability  n  .  (2.3-11) 

m 

For  example,  we  may  consider  the  case  that  the  parameters,  with 
probability  nx  =  1-e,  take  the  values 

(<V  V  r’  ^  =  (ai>  °i»  °s  °)»  prob’  =  1-e»  (2-3-12a) 

and  with  probability  n2  =  e, 

( Og »  »  £ )  —  (o2  ,  o2»  0,  5j(t)),  prob.  -  e.  (2-3-12b) 

The  parameter  values  in  (2.3-12a)  reflect  the  assumption  of 
a  "background"  Gaussian  noise  situation  with  equal  noise  power  at  the  two 
inputs  and  no  correlation.  The  parameter  values  in  (2.3-12b)  are  suggestive 
of  a  combined  background  plus  "impulsive"  Gaussian  noise  with  equal  power 
cr§  =  V2a|  at  each  input,  but  with  a  finite  correlation  between  the  noise 
inputs.  This  correlation  is  further  suggested  in  (2.3-12b)  to  be  a  function 
of  the  relative  delay  of  the  inputs,  or  the  direction  of  arrival  of  the 
impulsive  noise,  due  perhaps  to  discrete  events  located  in  a  specific 
direction  relative  to  the  two  sensors  whose  waveforms  are  being  correlated. 

2.3.2  Probability  integral 

For  cases  such  as  given  in  (2.3-12),  where  the  noise  powers  are 
equal  at  the  two  inputs  and  it  is  also  assumed  that  the  signal  amplitudes 
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are  equal  (A1=A2),  then  the  probability  that  the  correlator  output  exceeds 
a  threshold  is  [29] 


Pr{y  >  nl  =  (1-e)  Pr{y  >  n;  oa  =  ab  =  als  5  =  r  =  0,  =  A2) 

+  e  Pr{y  >  n;  aa  =  ab  =  axV,  r=0,  tfQ,  Ax  =  A2} 

=  (1-e)  f(n>  Pj  erf,  0) 

+  e  f  (n/V2;  p/V2,  afV2,  5) 


(2.3-13) 


where 


f(n;  pi  a2.  5)  9  Pr 


X2(2;  x2(2;  0) 


>  n 


=  Pr|x2(2;  ^  x2(2)  + 


(2.3-14a) 


The  probability  above  may  be  solved  in  terms  of  Marcum's  Q-function,  yielding 


f(  n  5  p»  a2,  c) 


•V 


4t \£  q2 


i+e  r 


-  p 


2iVa2 

1-5 


(2.3-14b) 
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2. 3. 2.1  False  Alarm  probability 


For  the  SNR,  p  =  0,  (2.3-14)  becomes 


f(n;  0,  a2,  ?)  =  exp 


2n/o2 
"  1+5 


Hz il 


exp 


2ri /a2  4n/o2) 

1-5  " 


=  Ulil 


exp 


2n/o2 

"  1+5 


(2.3-15) 


Thus  the  false  alarm  probability  for  the  case  considered  is 


P„  -  (l-e) 


FA 


exp 


2n MV2 

1+5 


(2.3-16) 


where  5  is  the  correlation  coefficient  of  the  sum  of  background  and 
"impulsive"  Gaussian  noise.  We  observe  that  (2.3-16)  is  very  similar  to 
(2.2-10),  the  P fA  for  the  quadrature  detector,  with  the  correlation  coefficient 
5  acting  as  V2,  that  is,  extending  the  threshold.  This  influence  of  5  can  be 
seen  in  the  thresholds  given  in  Table  2.3-1. 

2. 3. 2-2  Detection  performance 

By  substituting  (2.3-14b)  into  (2.3-13)  for  nonzero  SNR  p,  we 
obtain  the  correlation  detector's  probability  of  .detection : 


PD(n;  p)  =  (l-e) 


QUv'p’,  2/n h\) 
j  exp  j-p  +  2r\/ of |  Q(/2p",  2/2n/a^)j 
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GAUSSIAN 

GAUSSIAN-GAUSSIAN  NOISE 

NOISE 

e= 

:0.1 

e=0, 

,001 

Correl ation 
Coeff.,  5 

PFA 

(e;°  l 

n/crf 

V2=10 

n Jo\ 

V2=100 

n/cif 

V2=10 

n/ erf 

V2=100 

n/af 

0.1 

0.8047 

1.0143 

1.0880 

0.8214 

0.8249 

0.0 

0.01 

1.9560 

8.0472 

80.4719 

2.1478 

2.2756 

0.001 

3.1073 

19.5601 

195.6012 

8.0474 

80.4719 

0.1 

0.8047 

1.0546 

1.1389 

0.8239 

0.8275 

0.1 

0.01 

1.9560 

9.3761 

93.7611 

2.1819 

2.3256 

0.001 

3.1073 

22.0403 

220.4033 

9.3761 

93.7611 

0.1 

0.8047 

1.2552 

1.4179 

0.8344 

0.8382 

0.5 

0.01 

1.9560 

15.1118 

151.1177 

2.3484 

2.5956 

0.001 

3.1073 

32.3812 

323.8116 

15.1118 

151.1177 

Table  2.3-1  False  Alarm  Thresholds  for  Correlation  Detector 
(Single  Sample) 
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(2.3-15) 


For  e=0,  this  expression  reduces  to  the  Gaussian  noise  case. 

The  performance  of  the  correlation  detector  in  uncorrelated 
Gaussian  noise  is  given  in  Figure  2.3-2  as  a  reference.  When  compared  to  similar 
cases  for  the  (single  channel)  envelope  detector  as  presented,  for  example, 
in  the  V2  =  1  curves  of  Figure  2.2-7,  we  observe  that  the  correlation  detector 
achieves  a  50%  detection  probability  for  2.5  to  3.0  dB  less  SNR  than  required 
by  the  square-law  envelope  detector.  However,  the  correlation  detector  uses 
two  channels  and  thus  has  twice  the  signal  energy  to  use  for  detection.  Thus, 
in  Gaussian  noise  the  correlation  detector  may  be  said  to  be  -0.5  to  0.0  dB 
worse  than  the  envelope  detector  in  Gaussian  noise. 

Figures  2.3-3  and  2.3-4  show  the  detection  performance  of  the 
correlation  detector  in  bandpass  Gaussian-Gaussian  mixture  noise  (e  =  0.1) 
f0r  v2  =  10  and  V2  =  100,  respectively.  We  observe  from  these  figures  that 
positive  values  of  the  correlation  coefficient  g  in  (2.3-15)  degrade  the 
detector's  performance,  except  for  high  values  of  SNR  and  relatively  high 
false  alarm  probability  such  as  PpA  =  0.1.  We  also  note  that  the  performance 
is  degraded  in  proportion  to  values  of  the  variance  ratio,  V2. 

How  the  correlation  detector  performs  relative  to  the  square-law 
envelope  detector  in  the  same  Gaussian-Gaussian  mixture  noise  is  learned 
by  comparing  Figures  2.3-3  and  2.3-4  with  the  previous  figures,  such  as  2.2-7. 
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Figure  2.3-2.  Performance  of  bandpass  correlation  detector  in 
uncorrelated  Gaussian  noise. 
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Figure  2.3-3.  Performance  of  bandpass  correlation  detector  in  Gaussian- 
Gaussian  mixture  noise  (£=0.1,  V2=10)  for  different 
degrees  of  correlation  in  "impulsive"  component. 


70 


DETECTION  PROBABILITY 


Q.  9999- 
0.  9997- 

0.  9990' 

0.  9970- 

0.  9900- 
0.  9800' 

a.  9500' 

0.  9000 

0.  8000' 

0.  7000' 
Q.  6000 
0.  5000 
0.  4000 
0. 3000 

Q.  2000 

o.  laao 

0.  0500 

0.0200 

0.0100 

0.  0030 

0.  0010 

.0.  0003 
0.  0001 


SIGNAL  TO  NOISE  RATIO  CdB) 

Figure  2.3-4.  Performance  of  bandpass  correlation  detector  in  Gaussian- 
Gaussian  mixture  noise  (e=0.1,  V2=l 00 )  for  different 
degrees  of  correlation  in  "impulsive"  component. 
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For  50%  detection  and  V2  =  10,  we  find  that  the  correlation  detector  (after 
subtracting  3dB  because  of  the  two  channels  used)  effectively  requires 
0. 5dB  less  SNR  for  PpA=10”1  and  1.4dB  less  for  PpA=10-2.  When  V2  =  100, 
the  correlation  detector  in  effect  requires  2.4dB  less  for  Pp^lO-1  and 
1 . 6dB  less  for  PpA=l(T2.  Thus  in  this  type  of  noise  the  correlation 
detector  seems  to  have  a  slight  advantage  over  the  conventional  noncoherent 
Gaussian  detector  when  the  noise  is  uncorrelated  U  =  0).  However  when 
the  correlation  is  positive  (5  >  0),  this  advantage  is  decreased.  The 
interpretation  of  this  fact  seems  to  be  that  the  correlator,  in  using 
two  sensors,  acts  as  a  two-element  array.  When  this  array  is  "steered" 
toward  a  signal  source  such  that  01=02  and  such  that  the  (directional) 
impulsive  Gaussian  component  is  either  not  correlated  at  the  sensors 
(effectively  at  90°  or  270°  bearing)  or  negatively  correlated,  the  effect 
of  this  interference  is  lessened. 
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3.0  OPTIMUM  NONCOHERENT  DETECTION  IN  BANDPASS  GAUSSIAN-GAUSSIAN 

MIXTURE  NOISE 

In  the  previous  section  we  considered  the  performance  of  detectors 
which  are  designed  for  Gaussian  noise,  when  the  actual  noise  environment 
is  a  Gaussian-Gaussian  mixture.  Departure  from  Gaussian  conditions  was 
shown  to  result  in  a  loss  of  detector  performance  whose  severity  depends 
on  the  mixture  parameter  e  and  the  variance  ratio  o^/o\  -  V2.  It  was  shown 
also  that  using  multiple  samples  to  decide  whether  a  signal  is  present  can, 
but  does  not  necessarily,  improve  the  detector's  performance,  depending 
on  the  time  variability  of  both  signal  and  noise. 

Now  we  consider  the  performance  attainable  by  optimum  detectors 
for  signals  in  bandpass  Gaussian-Gaussian  mixture  noise,  that  is,  detectors 
based  on  the  generalized  likelihood  ratio,  and  describe  how  that  performance 
depends  on  knowledge  of  signal  and  noise  parameters. 

3.1  NON-GAUSSIAN  DETECTOR  FORMULATION 

As  a  particular  case  of  detection  of  signals  in  non-Gaussian 
noise,  we  turn  to  the  problem  defined  as  follows:  on  the  basis  of  the 
received  waveform  r(t),  '0<t$T,  we  wish  to  accept  or  reject  the  null 
hypothesis 

Ho^  Kt)  =  n(t)  (noise  only)  (3.1-1) 

when  the  noise  is  bandpass  Gaussian-Gaussian  mixture  noise, 

n(t)  =  nc(t)cosoj0t  -  ns(t)sinuj0t  (3.1-2) 
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where  f0  =  o)0/2n  is  the  center  frequency  of  the  band,  and  the  joint  pdf  of 
the  quadrature  components  n  (t),  n  (t)  at  a  given  instant  is  the  bivariate 
Gaussian-Gaussian  mixture  pdf 


1-e 


pnc,ns<“’s>  =  2^Texp 


o?~  +  )  l  £  l  ci^  +  3^ 

+  2^2 exp \—^r 


(3.1-3) 


Detection  occurs  when  Hq  is  rejected  in  favor  of  the  alternative 


hypothesis 


Hj :  r(t)  =  n(t)  +  Acos[u0t  +  e(t)],  (3.1- 

in  which  the  signal  amplitude  A  is  constant  during  the  observation  interval* 
and  the  signal  phase  e(t)  is  random.  Two  different  assumptions  will  be  made 
about  the  phase:  (a)  the  random  phase  is  constant  ("slowly  varying")  during 
the  observation  interval  (type  1  signal);  or  (b)  the  random  phases  of  samples 
taken  during  the  observation  interval  are  independent  (type  2  signal). 

We  assume  that  K  samples  of  rc  and  rg,  the  quadrature  components 
of  r(t),  are  taken  on  the  interval  (0,T).  Under  the  two  hypotheses,  the 
joint  pdf's  of  these  samples  are 


H0: 

prc.rs(a-  ^Ho)  =  pnc.ns(^ 

(3 . l-5a) 

Hi: 

prc,rs<i'  6|H„  I)  *  P„Cf„s<2.-V  !-is> 

(3 . l-5b) 

where  the  vector  notation  signifies 


*  As  noted  in  Section  2,  detector  performances  can  be  evaluated  for  non¬ 
constant  signal  amplitudes  in  a  straightforward  manner.  However,  we 
continue  to  assume  constant  amplitude  in  detector  formulations,  in 
order  to  simplify  the  analysis  somewhat. 
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Lc  =  [rc(ti),  rc(t2),...,rc(tK)] 
£s  =  [rs(t1)>  rs(t2),...,rs(tK)] 
If.  =  [sc(t1),  sc(t2),...,sc(tK)] 


(3 . l-5c) 


=  Acose  [1,  (Type  1) 

=  A [cose x ,  cose2,. . . ,cos9k]  (Type  2)  (3 . l-5d) 

Is  =  [ss(t1),  ss(t2),...,ss(tK)] 

=  Asine  [1,  1....1]  (Type  1) 

=  A[sinels  sine2,...,  sine^]  (Type  2).  (3.1-5c) 

The  test  for  rejecting  H0  in  favor  of  \ix  is  to  be  based  on  the 
generalized  likelihood  ratio  (GLR) 

Ee{Prc,rs(i’SlHi'i>| 

=  PrCiI.s(2.,6lH0)  •  (3.1-6) 

3.1.1  Conditional  Gaussian  Approach 


The  noise  pdf  (3.1-3)  for  a  single  sample  may  be  viewed  as  the 
average  over  the  variance  of  a  conditionally  Gaussian  pdf.  Let  v2=  oz/o\ 
be  a  variance  multiplication  factor;  then  (3.1-3)  can  be  understood  as 


^nc,ns 


(3. l-7a) 
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with 


PJO  =  (l-e)fi(Y-l)  +  £  6(Y-V2);  V  =  — j  . 

1 


(3. l-7b) 


Extension  of  this  concept  to  multiple  samples  takes  the  form 


(2’°i)’K  .J,  £  “k  +  Bk 


Pnc,n_s^-’-^  Ev2]  vivi  •  *  -vk  8XP 


k=l  2vH 


(3 . l-8a) 


with 


V'L>  -  I  c„  I5™'1 


m=l 


m=l 


(3. l-8b) 


Using  this  form,  the  expectation  of  the  Hi  joint  pdf  over  the  signal 
phase  becomes 


Eilprc,rs(2-'^Hi-5.)| 


=  E„ 


(2ira|)_ 

exp 


r*  (“k-Acos9k)2  +  K-Asin6k'2 

L  2vpf 

k=l  K 


=  E. 


1  1 

(2™! ) 

c 

5 

( 

1 

Y'  ak +  ek +  A2 

]v?vl  •• 

’  H  2vk3^ 

k=l  K  J 

x  E0{exp 


I 

k=l 


akcosek  +  ek  sinek 


v2 


(3.1-8c) 
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3.1.2  GLR  For  Constant  Signal  Phase 

For  type  1  signals,  the  above  operations  yield 


(3.1-9) 


With  this  result  the  GLR  becomes 


Ar(a,6.)  - 


Z 

vm=l 


C 


1 


m  vlmv2m' 


Tv2"- 

Km 


exp 


z 

k=l 


+  Sf 


-l\-l 


2 v2  a2 
km  1 


Z 

,  m=l 


■'m  ~\T1 

m  vlmv2m‘ 


Km 


exp 


-z 

k=l 


“Zk  + 


8j[  +  A2 


2v2  a2 

km  l 


M 

=  ^  Wm  ^  ’ 

m=l 


(3.1-10a) 


(3.1-10b) 


where 
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C  fv2  v2 


W  (a  ,6)  = 
m - 


mV  lm  2m ’ "  km 


exp 


K 

E 

k=l 


i2 

‘k 


(3.1- 10c ) 


y  (numerator) 
m=l 


and 


Am(— 


=  exp 


-z 

k=l 


A2 


2  a2V? 


1  km 


km 


=  exp 


-p  £v 

k=l 


km 


using  p  =  A2  /ic\  . 

3.1.3  GLR  For  Independent  Phase  Samples 

For  the  Type  2  signal,  each  0^  e(0,  2-rr )  is  assumed  to  be  independent, 
yielding  the  Hi  joint  pdf 


Kr  r 
— c s 


1  K 

(27ra^  V1 

f  A\yr\ 

-  c.2  +  62  +  A2‘ 

- 

TT 

k=l 

* - 1 

O 

A  L,  2  +  p2 

exp 

l  2aivk  J 

l^TTk  kJ 

(3.1-11) 


and  the  GLR 


A  (a,e)  = 


E  2{pr  r  («,6.;A|v2)J 

C  >  S 


r  —  —  E 


:i2|prCir$l-'i;0|-")| 


(3.1-12a) 
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with  W  (a 

m  — 


3-1.4 

vectors 


(3 . 1- 12b ) 


M 

*  Y.  Vl'i)  (3-1-  12c ) 

m=l 


given  in  (3.1-lOc),  and 


exp 


4^; 


(3. l-12d) 


=  exp 


-p 


(3.1-1 2e ) 


Example 

]/ 

For  example,  let  K=2.  Then  there  are  M=2  =4  possible  variance 
(vlm>v2m . vMm): 


vf  =  (1,  1),  v|  =  (V?  1),  v23=  (1,  V),  vj  =  (V2  t)  (3.1-13a) 


79 


J.  S.  LEE  ASSOCIATES,  INC. 


with  probabilities 

Cm  =  p>-jv2=4|}  ;  ■-!.  2,...,H=2K. 

In  this  case,  the  numerator  and  denominator  of  the  GLR  become 


(3. l-13b) 


Ci 

-4  exp 


al  +  ^1  +  a2  +  ^2  +  ^A2 


2af 


I, 


°AVv 


ai  +  a2\z  /B i +  &2\2 


G2 

+  °FT exp 


of  +  gf  +  A2  of  +  ef  +  A2 


2af 


M  ^X2  /Bi  e2\2 


+  — t*  "'■s  exp 

al°2 


af  +  gf  +  A2  af  +  8f  +  A2 


2af 


2a| 


’I 


a2  / 


,ai 


+  ^  exp 


a2  +  82  +  af  +  +  2A2 


al +  a2  \  /Si +  32 


(3. l-14a) 


and 


T^r  exp 


al  +  S2  +  Otf  +  $2 


■^r 


+  T^z  exp 


a20l 


af  +  Sf 


a*  + 


"T 


C3 

a2  +  g2 

+  S| 

C  4 

+  exp 

°2 

of +  ef  +  of  +  ef 

— 7”“^ —  6XD 

°ia2 

2a2 

2cf  ‘ 

2a2 

(3. l-14b) 
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3.2  EFFECT  OF  VARIANCE  DISTRIBUTION  ASSUMPTIONS 

So  far  we  have  not  specified  the  joint  pdf  p  2(y)  i  except  to 
describe  it  in  (3.1-8b)  as  taking  discrete  vector  values  {v^},  possibly  2  =M 
values  since  each  noise  sample  has  a  discrete  variance  multiplier  pdf  with 
two  possible  values  (1  or  V2).  The  joint  pdf  should  reflect  the  behavior 
over  time  of  the  non-Gaussian  noise  process,  which  we  are  modelling  by 
a  bandpass  Gaussian-Gaussian  mixture  process. 

For  example,  at  one  extreme  we  may  say  that  the  probabilities 

Cm  of  (3.1-8b)  are 

Ci  =  1  a  Cm=0,  m  >  1.  (3.2-1) 

This  corresponds  to  a  single  variance  °2v£=a2  f°r  the  samples,  or, 
in  effect,  Gaussian  noise.  The  GLR  (3.1-10)  then  becomes 


(3.2-2) 

As  illustrated  previously  in  Figure  2.2-l(a),  implementation  of  this  Gaussian 
GLR  for  the  type  1  signal  requires  separate  summing  of  samples  (integration)  in 
the  two  quadrature  channels,  then  forming  the  combination  of  the  squares  of 
these  sums  and  comparing  this  quantity  to  a  threshold.  For  the  type  2  signal, 
the  combining  of  the  squares  takes  place  before  summing,  as  shown  in  Figure 
2 . 2-1 ( b) . 

3.2.1  Slowly-varying  Noise  Power. 

Another  limiting  case  is  described  by  Cx  =  1-e;  CM  =e  ;  Cm  =  0, 
m  f  1,  M.  That  is. 
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Pv2(y)  =  (l-e)6[X  “  (1»  1 . D1 

+  e  6[y  -  (V2  V^. . . ,V2)  ] ; 


(3.2-3) 


the  a2v2  are  all  equal  to  a2  with  probability  1-e  and  all  equal  to  a2  with 
probability  e.  This  case  corresponds  to  a  "slowly-varying"  non-stationary 
Gaussian  noise  model  with  two  possible  variances.  For  this  case,  the  GLR 
is  (Type  1  signal) 


Ar(“>JL) 


e,  =  l- 


[1  -  W(a_,j3)]e 


KA2/2a? 

Io 


_A_ 

a2 

1 


1.  h 

k=l  / 


where 


+  W(a,Bj 


(3.2-4a) 


W(a_,3_) 


(l-e)ai 


2  K 


e 


-2K 

a2 


exp 


- r 

exp 


K 

I 


k=l 


exp 


(3.2-4b) 

We  note  that  the  noise  samples  for  this  slowly-varying  nonstationary  model 
are  not  independent,  although  they  are  uncorrelated.  As  shown  in  Figure  3.2-1, 
the  appropriate  detector  is  a  modification  of  the  detector  for  Gaussian  noise. 
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W ( Z2  i  £»  ^1 >  ^2 )  — 


e  ojj^  exp  {-z2/2cr|} 


(l-e)o^^  exp  {-z2/2of}  +  ea2^  exp  {-z2/2a|} 


-2K 


Figure  3.2-1  Optimum  detector  for  sinusoidal  signal  in  slowly-varying  Gaussian-Gaussian 
mixture  noise. 
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3.2.2  Independent  Noise  Samples 

Independence  of  (a^,  3^)  quadrature  sample  pairs  implies  the 
following  set  of  probabilities  in  the  pdf  for  the  variance  vector: 

Ci  =  Pr(v>  (1,...  ,1)}  =  (l-e)K 

C2  =  pr{v2=  (V2l,...l)}  =  e(l-£)K_1  (3.2-5) 

CM_!  =  Prlv2=  (V2. . .  ,V2  1)}  =  e^-1(l-e) 

CM  =  Pr{v2=  (V2V2...,V^>  =  £K 


In  general , 


(1-e) 


K-w 


m 


(3 . 2-6a ) 


where  w  is  the  number  of  one's  in  the  binary  representation  of  (m-1): 
m 

wm  =  weight! (m-l)2l •  (3.2-6b) 

We  observe  that  direct  implementation  of  the  GLR  (3.1-10)  for  the 
case  of  Type  1  signal  and  independent,  multiple  samples  appears  to  be 
undesirable  for  numbers  of  samples  of  appreciable  size,  since  the  number  of 
separate  LR's  A  (  )  in  (3.1-10b)  grows  exponentially  with  the  sample  size 
(M  =  2K). 

For  Type  2  signals  and  independent  noise  samples,  the  weighting 
concept  expressed  by  (3.1-12)  becomes  cumbersome,  since  it  is  simpler  to 
express  the  GLR  as  the  product  of  single-sample  GLR's: 

A  (a,  3)  =  n  A(a. ,  3.),  (3.2-7a) 

1  k=l  K  K 
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where 


A(ct  ,  6  ) 
k  k 


(l-£)ai2  exp 


.2  o2 


k  +  k  + 


-2 

+  e  a2  exp 


2al 

■i  *  ♦  *2 


L°i 


(l-£)ai2  exp 


ak  ±  Sk 


2a| 


-2 

+  £  a2  exp 


“k  +  ek 


2o2 


(3.2-7b) 


The  implementation  of  (3.2-7)  can  be  accomplished  with  accumulators 
as  illustrated  in  Figure  3.2-2,  since  the  likelihood  ratio  test  is 
equivalent  to 

K 


H, 


1°9  Ar  .  -  ^  log  A(ak,  Pk)  <  n'. 


(3.2-8) 
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BPF 


In  A(z;  A,  e,  of,  a\)  =  In 


In 


(l-e)ol2  exp  {-  I„  [Ar  /z]  +eoiexp  j-  §±£j  I„[^  /f] 

(l-e)oj2  exp  j-z/2o^|  +  e  o|  exp  |-z/2o^|| 


Figure  3.2-2.  Optimum  detector  for  noncoherent  signal  and  independent  Gaussian-Gaussian 
mixture  noise  samples. 
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3.3  PERFORMANCE  OF  SINGLE-SAMPLE  DETECTORS 

Before  evaluating  the  performance  of  detectors  based  on  the  GLR 
for  multiple  samples  (K>1),  in  this  section  we  evaluate  the  detection 
performances  of  the  GLR  and  other  detectors  for  a  single  sample  pair  of 
quadrature  components  in  order  to  develop  an  understanding  of  the  non- 
Gaussian  detection  problem.  For  a  single  sample,  the  Type  1  and  Type  2 
signal  phase  models  are  the  same. 

3.3.1  Form  of  GLR  Detector 


For  one  sample  (K=l),  the  GLR  in  (3.1-10)  becomes 


r 

-A2/2a2 

"  A  - - 

[l-w(Vs)J 

e  I0 

—  Jr2  i-  r2 
_ai  V  c  s 

-A2/2 a\ 


+w(>vrs)e  2l°£V 


.4.  y.2 

a  2  Y  c  s 


(3.3-la) 


where 


e  a 22  exp 

i+1 

(1-e)  a“2exp  * 

I  rc  +  rs  ( 

+  £  aj2  exp 

rz  +  r2 

c  s 

“  2o?  , 

2a| 

(3.3-lb) 


We  observe  that  the  GLR  is  a  function  of  the  envelope  only,  so  that  its 
implementation  involves  formation  of  x=  R2  =  r2  +  r2,  followed  by  a  nonlinear 
function: 
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Ar  =  f(x;  e,  V2,  a2,  p) 


( 1-e) V2  exp 

X 

- -  p 

I  o  ( ^2p  x/ <7 ! )  +  e  exp 

1  f^T  +  p)| 

2a2 

V2  \ 2a?  /' 

I0  (  /2o x/c^V2) 


(l-e)V2  exp  | - — 

+  e  exp 

x  ) 

1  2a2 

2a2V2  i 

(3.3-2a) 

V2  =  —  and  P  =  —  . 

(3.3-2b) 

using 


2a2 


This  function  is  parametric  in  signal  as  well  as  in  noise  parameters.  To 
illustrate  the  influences  of  the  various  parameters,  plots  of  (3.3-2)  are  shown 
in  Figures  3.3-1  to  3.3-6.  In  each  figure,  £  and  V2  assume  given  values 
while  p  is  the  parameter  indexing  the  family  of  curves,  which  are  plotted 
as  functions  of  x/crf,  the  squared  envelope  sample  value  normalized  by  the  smaller 
of  the  two  mixture  noise  powers. 

The  curves  in  Figures  3.3-1  to  3.3-6  show  that  the  form  of  the 
detector  for  this  type  of  non-Gaussian  noise  is  very  dependent  on  the  a_ 
priori  value  of  p,  the  SNR  which  pertains  when  the  signal  is  present.  This 
behavior  is  significantly  different  from  the  Gaussian  case  (single  sample) 
in  only  one  respect:  the  characteristic  is  not  monotonic.  What  is  meant 
by  this  statement  can  be  explained  as  follows: 
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-Figure  3.3-1.  Likelihood  ratio  for  Gaussian-Gaussian  mixture  noise 

( e  =  0 , 0 1 ,  V2  =  10),  parameterized  by  SNR. 
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LIKELIHOOD  RATIO 


to 


Figure  3.3-2.  Likelihood  ratio  for  Gaussian-Gaussian  mixture  noise 

(e  =  0.01 ,  V2  =  100) ,  parameterized  by  SNR. 
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Figure  3.3-3.  Likelihood  ratio  for  Gaussian-Gaussian  mixture  noise 
(e=0.01,  V2  =  1000),  parameterized  by  SNR. 
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NORMALIZED  ENVELOPE  SQUARED 

Figure  3.3-4.  Likelihood  ratio  for  Gaussian-Gaussian  mixture  noise 


(£=0.1,  V2=10),  parameterized  by  SNR. 
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NORMALIZED  ENVELOPE  SQUARED 

Figure  3.3-5.  Likelihood  ratio  for  Gaussian-Gaussian  mixture  noise 
(e  =0.01,  V2  =  10),  parameterized  by  SNR  (detail  of  Figure  3.3-4). 
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LIKELIHOOD  RATIO 


NORMALIZED  ENVELOPE  SQUARED 


Figure  3.3-6.  Likelihood  ratio  for  Gaussian-Gaussian  mixture  noise 

(e=0.1»  V2  =  100),  parameterized  by  SNR. 
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Even  in  the  Gaussian  case,  straightforward  formulation  of  the 
GLR  gives  a  function  which  is  dependent  on  SNR,  namely 


A(x;p)  =  e~pI0  (j  /2px).  (Gaussian  GLR)  (3.3-3) 


The  false  alarm  probability  for  this  case  is 


P 


FA 


Pr 


Pr 


x >  n  (a,p)  |H0  , 


(3.3-4) 


since  a(x;p)  is  monotonic.  That  is,  the  probability  that  the  GLR  is  greater 
then  some  threshold  X  is  entirely  equivalent  to  the  probability  that  R2=x  is 


greater  than  another  threshold, 

n(a,p) 


I'0‘  (efx(p)) 


n  >  0, 


(3.3-5) 


using  the  notation  IoM*)  to  indicate  the  inverse  of  the  function  I0(  )•  In 
Neyman-Pearson  detection,  the  Pp^  is  fixed  at  some  value  a.  This  in  turn  fixes 
n(a,p)  =  n  (o)  tor  all  values  of  p.  Therefore,  in  the  Gaussian  case  no  matter 
what  the  SNR  is,  the  receiver  decision  is  based  on  comparing  x  to  the  false 

alarm  threshold  n  .  The  same  statement  is  true  for  all  GLR's  which  are  monotonic 

a 

functions  of  x.  It  follows  that  the  single-sample  detection  probabilities  for 
all  decision  variables  which  are  monotonic  functions  of  x  are  also  equal,  since 


PD  =  Pr 


f(x;p,a)  >  Aa(p,a)|Hi 


=  Pr 


x>  na(cr)  |HX 


(3.3-6) 


for  all  of  them.  Now,  as  illustrated,  in  Figures  3.3-1  to  3.3-6,  the  GLR  for 
Gaussian-Gaussian  noise  is  not  monotonic;  consequently,  the  inverse  mapping 
of  the  GLR  to  x  yields  false  alarm  thresholds  which  depend  on  the  SNR. 

A  further  description  of  the  dependence  of  the  detector  structure 
on  the  a  priori  SNR,  p,  is  provided  by  the  following  considerations.  Suppose 
the  desired  false  alarm  probability  is  Pp^=a;  for  monotonic  GLR's  this  require- 
. ment  is  equivalent  to  x>  na  with  probability  a.  But  for  non-monotonic  GLR's 
the  requirement  in  general  is 
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a  =  Pr  | A(x;p)  >  X(p)j 

=  Pr  |m(p)  <  x  <  n2  (p)|  +  Pr  |  x  >  n3 (p )  |  ,  (3.3-7) 

as  illustrated  in  Figure  3.3-7a;  the  region  of  x  corresponding  to  A>X  involves 
three  thresholds.  However,  for  moderate  p  values,  the  non-monotonic  behavior 
of  a(x;p)  is  confined  to  a  certain  region  of  (A,x)  as  shown  in  Figure  3.3-7b, 
and  if  the  quadrature  detector  false  alarm  threshold  na  (Table  2.2-1  and  Figures 

2.2- 2  to  2.2-5)  falls  outside  this  region  (na  <  n  or  na  >  na^  in  Figure 

3 .3 - 6),  the  performance  of  the  GLR  receiver  is  the  same  as  for  square-law 
envelope  detector  For  large  a  priori  p,  then  we  expect  the  non-Gaussian  GLR 
to  yield  the  same  performance  as  the  square-law  envelope  detector. 

Implementation  of  the  GLR  for  a  single  sample  can  be  performed  in 
concept  using  the  configuration  diagrammed  in  Figure  3.3-8.  This  implementation 
approach  makes  use  of  the  mixture  form  of  the  GLR  discussed  in  Section  1.2.3. 

We  observe  that  two  function  generators  are  required,  one  for  the  parametric 
likelihood  ratio  and  one  for  the  weighting  function.  It  is  also  clear  from 
the  diagram  that  the  a  priori  information  required  consists  of  the  parameters 
of,  e,  V2,  and  p.  In  practice  these  parameters  may  be  imperfectly  estimated; 
however,  for  the  present  we  shall  assume  that  they  are  known  or  estimated 
precisely. 

3.3.2  False  Alarm  Probability 

The  generalized  likelihood  ratio  (3.3-2)  depends  on  the  SNR  p  in 
such  a  way  that  an  equivalent  statistic  not  parametric  in  the  SNR  does  not 
exist  for  testing  the  hypotheses  Hi  (signal  present,  with  given  SNR)  against 
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Fig  3.3-7  Dependence  of  detector  false  alarm  thresholds  on  shape  of 
likelihood  ratio. 
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A(x)  =  [l-W(x) ]  Ax(x)  +  W(x)  A2(x) 

where  A l{x)  =  f  (X,  p),  A2(x)  =  f(-^-  ,  ^  ) 

f  (a ,  b)  =  e“b  I0(/2ib) 

and  W(x)  =  e  exp{-x/2V2}  /  [(l-e)V2  exp  {-  |}  +  e  exp  {-  ~pr)] 


FIGURE  3.3-8  CONCEPTUAL  IMPLEMENTATION  OF  GLR  DETECTOR 
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the  hypothesis  H0  (noise  only).  That  is,  a  uniformly  most  powerful  (UMP) 
test  statistic  does  not  exist  for  this  problem  ([19],  p.  91).  Therefore  it 
is  necessary  to  distinguish  between  an  a  priori  or  assumed  SNR,  which  we 
shall  denote  by  p0,  and  the  actual  SNR,  which  we  continue  to  denote  by  p. 

With  this  understanding,  the  false  alarm  probability  for  bandpass 
Gaussian-Gaussian  mixture  noise  can  be  expressed  by 

Pra  =  Pr  {Ar(x'>  Po)  >  *|p-0> 


=  Pr  (x  e  Rx  (A,  p0)  |  p=0} 


(3.3-8) 


where  x  is  the  sampled  envelope-squared,  X  is  a  threshold  for  the  likelihood 
function;  and  R  (x,pQ)  is  region  for  x  defined  by  X  and  p0.  As  explained 
previously  with  the  help  of  Figure  3.3-7,  this  region  involves  either  one 
threshold  or  three  thresholds: 


R*  (x,p0) 


x  >  m,  n  =  A-1  (X;  p0)  single-valued; 

< 

Hi  <  X  <  Tl2»  X.  >  H3f 


(3.3-9) 


n  =  A-1  ( X ;  pQ  )  multiple-valued 


Thus  the  single-sample  detector  may  be  implemented  simply  by  comparing  the 
squared-envelope  sample  value  to  the  threshold(s) ;  it  is  not  necessary  to 
implement  the  GLR  directly,  as  diagrammed  in  Figure  3.3-8. 
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Explictly,  the  false  alarm  probability  is,  using  (3.3-9)  in 

(3.3-8), 

PFA  =  po  (ni),  n  =  A-1  U;  p0)  single-valued;  (3.3-10a) 


where 


=  po  (ni )  -  Po  (n2)  +  po  (^3) > 

n  =  A-1  (A;  p0)  multiple-valued;  (3.3- 10b ) 


po  (n)  *  (1-e)  e' 


■n/2 0? 


+  e  e 


■n/2a^V2 


(3.3-11) 


Since  P0(n)  is  identical  to  the  false  alarm  probability  for  the  Gaussian 
detector  (2.2-10),  we  can  observe  from  (3.3-10)  that  the  PpA  for  the  optimum 
detector  is,  for  the  same  threshold  m»  l-ess  than  or  equal  to  that  obtained 
using  the  Gaussian  detector.  It  is  less  when  m  falls  in  the  non-monotonic 
region  of  the  likelihood  ratio,  and  equal,  otherwise.  This  interpretation 
of  (3.3-10)  is  confirmed  by  the  curves  shown  in  Figures  3.3-9  to  3.3-12. 

Figures  3.3-9  and  3.3-10  give  the  false  alarm  probability  for 
the  SNR-dependent  GLR  when  e  =  0.01  and  p0  =  0  dB  and  15  dB,  respectively. 

In  order  to  compare  with  previous  results,  we  have  plotted  PpA  vs  m» 
the  smallest  of  the  thresholds  when  there  are  more  than  one.  These  two 
figures  correspond  to  the  case  shown  previously  in  Figure  2.2-2  for  the 
Gaussian  receiver.  The  most  obvious  effect  seen  in  the  new  figures,  as 
compared  with  the  old,  is  that  the  false  alarm  thresholds  are  reduced  by 
orders  of  magnitude.  The  effect  is  especially  pronounced  for  high  SNR; 
the  sudden  decreases  in  PpA  as  the  threshold  increases  occur  when  the 
threshold  X  =  A“1(ni;p0)  begins  to  fall  in  the  non-monotonic  portion  of  the 
GLR  characteristic  (compare  Figures  3.3-3  and  3.3-10). 
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Gaussian-Gaussian  mixture  noise  (e  =  0.01),  with  assumed  SNR  of 

0.0  dB  and  parametric  in  variance  ratio. 
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a.  Gaussian-Gaussian  mixture  noise  (e=0.01),  with  assumed  SNR  of 

a 

a 

15  dB  and  parametric  in  variance  ratio. 
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Similar  results  for  e  =  0.1  are  given  in  Figures  3.3-11  and 
3.3-12  for  p0  =  0  and  15  dB,  respectively.  These  may  be  compared  to  the 
Gaussian  detector's  false  alarm  probability  as  shown  previously  in 
Figure  2.2-3,  and  the  threshold  values  for  which  the  sudden  decreases  in 
PpA  occur  can  be  identified  with  the  non-monotonic  GLR  behavior  for  e  =  0.1 
given  in  Figure  3.3-7  for  p0  =  15  dB. 

3.3.3  Detection  Performance. 

For  fixed  PpA  and  assumed  or  a  priori  SNR  value  p0,  the  probability 
of  detection  for  the  optimum  detector  in  Gaussian-Gaussian  mixture  noise 
can  be  written 

PD  =  PD^p;  PFA’  p°^ 

=  Pr{A  (x;  p0)>A  | p^O}  (3.3-12) 

i  ot 

where  A  is  defined  by  the  constraint 

a 

PFA  =  Pr{Ar(x;  Po)>Aa|p=0}=a  .  (3.3-13) 

As  discussed  in  the  previous  section,  the  event  of  the  likelihood  ratio 

A^(x;  Po)  exceeding  some  threshold  A^  is  equivalent  to  the  squared  envelope 

x  being  in  some  region  R  (A  ,  p0),  which  can  be  defined  by  a  single  threshold 

n  or  by  three  thresholds  (m>  n2>  n3),  depending  on  whether  the  GLR  is 

monotoni'c  at  A  =  A  .  Therefore  the  Pn  can  be  calculated  by  the  expression 
r  a  u 

PQ  =  P1(p;  m)  ,  n  =  A-1  (Aa;  p0)  single-valued;  (3.3-14a) 

=  Pi(p;  ni)  -  Pi (p ;  n2)  +  pi(p;  n3), 

n  =  A-1 ( Aq;  Po)  multiple-valued;  ( 3 . 3- 14b ) 
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Figure  3.3-11.  False  alarm  probability  for  optimum  detector  in 
Gaussian-Gaussian  mixture  noise  (e=0.1),  with  assumed  SNR  of 
0.0  dB  and  parametric  in  variance  ratio. 


104 


FALSE  ALARM  PROBABILITY 


Q 

O-, 


PJ 


or 

CL 


EPSILON  =  0.  10 
RHO  =  15.  0  dB 

V2  =  a22/a12 


I 

o- 


i  mm — 

1C1 


Mrrr 


'  ’  "I  . 

10; 


-r~ it  r »'  1 1 


lCr 


10‘ 


104 


NORMALIZED  THRESHOLD 

Figure  3.3-12.  False  alarm  probability  for  optimum  detector  in 
Gaussian-Gaussian  mixture  noise  (e  =0.1),  with  assumed  SNR  of 
15  dB  and  parametric  in  variance  ratio. 
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where 

P 1  ( P ;  n)  £  (1-e)  Q(/2p,  +  eQ(S2^7W,  A, /af  V*),  (3.3-15) 

and  Q  (*,  •)  is  Marcum's  Q-function. 

3.3.3. 1  Results  for  known  SNR. 

The  best  detector  performance  can  be  expected  when  the  a_ 
priori  SNR,  p0,  is  correct,  that  is,  equal  to  the  exact  SNR,  p.  Calculations 
of  this  best  performance  were  made  using  (3.3-14);  since  we  assume  that 
p  =  p0,  different  false  alarm  thresholds  were  obtained  for  each  value  of  SNR 
(The  computer  program  used  is  listed  in  Appendix  3A.)  For  e  =  0.01,  the 
detection  probability  for  these  assumptions  varies  with  p  as  shown  in  Figure 
3.3-13.  In  this  figure  the  is  shown  for  Pp^  =  10"1,  10  2  and  10  3,  and 
for  V2  =  1,  10,  and  100.  For  V2  =  1,  the  noise  becomes  Gaussian  and  GLR  reduces 
to  the  optimum  Gaussian  detector,  so  that  we  can  observe  from  the  V2  f  1  curves 
the  effects  of  the  non-Gaussian  parameters  e  and  V2. 

The  comparable  performance  results  for  the  single-sample  Gaussian 
detector  were  shown  previously  in  Figure  2.2-6.  In  comparison  with  those 
results,  we  observe  the  optimum  detector  in  Gaussian-Gaussian  mixture  noise 
for  e  =  0.01  performs  about  as  well  for  (a)  Pp^  =  10  1 ,  and  (b)  V2  =  10  and 
Pq  >  .5.  For  Pq  <  .5  a  great  improvement  is  accomplished,  particularly 
as  V2  increases.  For  example,  in  Figure  2.2-6,  we  find  that  an  SNR  in  excess 
of  20  dB  is  required  for  V2  =  100  to  achieve  Pp^  =  10  3  and  P^  >  3x10  3, 
while  in  Figure  3.3-13,  a  PQ  of  0,5  can  be  achieved  for  the  same  Pp^  when 
p  =  9  dB. 
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SIGNAL  TO  NOISE  RATIO  (dB) 


Figure  3.3-13.  Receiver  operating  characteristics  for  optimum  detector 
in  bandpass  Gaussian-Gaussian  mixture  noise  (e  =  0.01), 
for  different  false  alarm  probabilities  and  variance 
ratios. 
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Perhaps  the  most  interesting  feature  in  Figure  3.3-13  is  the 
behavior  of  the  as  a  function  the  SNR  for  V2  =  100  and  P^_  =  10  3. 

For  p  less  than  about  12  dB,  the  Pp  increases  with  p;  it  increases 
again  for  p  >  17  dB,  but  for  12  dB  <  p  <  17  dB,  the  PQ  decreases  with  p. 

From  the  mathematical  expression  for  the  Pp,  equation  (3.3-14),  we  observe 
that  such  a  decrease  is  possible  since  the  second  term  is  negative.  Because 
we  are  assuming  p0  =  p,  that  is,  the  a  priori  SNR  equals  the  actual  SNR, 
the  thresholds  (m,  n2,  n3)  change  as  p  changes  in  order  to  maintain  a 
constant  false  alarm  probability.  Evidently,  these  thresholds  change 
within  such  a  way  as  p increases  as  to  produce  the  "dip"  in  PQ  we  observe 
in  Figure  3.3-13. 

For  fixed  thresholds  (m,  n2,  n3),  which  corresponds  to  having 
both  p0  and  x  fixed,  we  expect  that  PQ  will  "dip"  because  the  action  of  the 
likelihood  ratio  characteristic  is  to  "pass"  or  accentuate  certain  values 
of  the  detected  envelope,  and  to  "suppress"  other,  high  values.  That  is, 
the  detector  discriminates  against  a  range  of  high  values  of  x,  in  effect 
considering  them  to  be  due  to  noise  impulses. 

When  the  mixture  parameter  is  increased  to  e  =  0.1  or  to  0.5, 
representing  a  greater  departure  from  Gaussian  noise,  the  optimum  detector 
performs  as  shown  in  Figure  3.3-14  and  3.3-15,  in  contrast  to  that  of  the 
Gaussian  detector,  shown  previously  in  Figure  2.2-7  and  2.2-9.  We  observe 
that  for  V2=  10,  the  detection  probability  is  practically  the  same  for 
both  detectors,  except  for  smaller  values  of  SNR,  at  which  the  optimum 
detector's  characteristic  is  not  a  mono tonic  function  of  the  detected 
envelope  for  the  false  alarm  thresholds  considered  (see  Figure  3.3-4). 
However,  the  optimum  detector  performance  is  much  improved  for 
V  2  =  100  and  PpA  <  10  "  l.  It  is  notable  that  for  much  of  the  range  of  SNR 
values  shown,  the  performance  of  the  detector  is  better  for  V2=  100  than 
for  V  2  =  10. 
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SIGNAL  TO  NOISE  RATIO  (dB) 


Figure  3.3-14.  Receiver  operating  characteristics  for  optimum  detector 

in  bandpass  Gaussian-Gaussian  mixture  noise  (e  =  0.1)  for 
different  false  alarm  probabilities  and  variance  ratios. 
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Figure  3.3-15. 


Receiver  operating  characteristics  for  optimum  detector 
in  bandpass  Gaussian-Gaussian  mixture  noise  (e  =  0.5) 
for  different  false  alarm  probabilities  and  variance 
ratios . 
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3. 3. 3. 2  Results  for  fixed  a  priori  SNR. 

Since  it  is  usually  the  case  that  the  SNR  of  the  signal  to  be 
detected  is  not  known,  the  detector  performances  shown  in  Figures  3.3-13 
to  3.3-15  must  be  considered  as  upper  bounds  to  what  may  be  realized  in 
practice.  Now  we  consider  the  effect  on  detection  performance  as  a  function 
of  actual  SNR  when  the  detector  GLR  characteristic  A^(x;  p0)  is  a  fixed 
design,  due  to  assuming  that  the  SNR  takes  a  certain  fixed  value,  p0. 

For  given  Pp^  and  p0,  this  results  in  a  comparison  of  the  squared  envelope 
x  to  a  threshold  n.  or  possibly  to  three  thresholds  (nls  n2>  03),  using 
the  test  (3.3-9)  discussed  previously. 

A  typical  plot  of  detection  probability  vs  SNR  for  fixed  p0 
is  shown  in  Figure  3.3-16,  for  the  case  of  e  =  0.1  and  V2  =  100.  For 
both  Pp£  =  10  1  and  10  2,  the  Pp  curves  for  p0  =  -10  dB  and  p0  =  +10  dB 
are  given,  and  are  compared  to  the  comparable  optimum  (p0  =  p)  and  Gaussian 
detector  results.  As  expected,  the  fixed-p0  P^  achieves  the  best  performance 
at  the  points  for  which  p0  =  p,  the  actual  SNR.  When  the  actual  SNR  is  not 
equal  to  p0,  we  may  distinguish  two  different  consequences,  depending 
on  p0. 

For  p0  =  -10  dB,  the  detector  is  predicated  on  the  assumption  of  a 
weak  signal.  From  Figure  3.3-16  we  observe  that  the  resulting  P^  values 
are  almost  as  high  as  the  optimum  values,  for  actual  SNR  as  much  as  5  dB, 
or  15  dB  different  than  the  assumed  value!  However,  the  P^  then  falls  to 
very  low  values,  less  than  the  false  alarm  probability,  before  rising  again 
for  extremely  high  SNR. 


Ill 


DETECTION  PROBABILITY 


Figure  3.3-16  Receiver  operating  characteristics  for  optimum 
detector  in  bandpass  Gaussian-Gaussian  mixture 
noise:  effect  of  using  fixed  value  of  a_  pri ori 
SNR. 


112 


J.  S.  LEE  ASSOCIATES,  INC. 


For  p0  =  +10  dB,  the  performance  falls  to  low  values  more  quickly 
when  p  >  p0,  that  is,  it  is  more  sensitive  to  p0  when  p  >  p0.  But  the 
degradation  in  PQ  for  p  <  p0  is  relatively  minor:  the  degraded  performance 
is  still  better  than  that  achieved  by  the  Gaussian  detector. 

These  results  suggest  that,  when  a  fixed  value  of  a  priori  SNR 
is  used,  it  should  be  somewhat  higher  than  the  average  value  of  SNR 
anticipated. 

3.4  PERFORMANCE  FOR  MULTIPLE  SAMPLES. 

We  have  observed  for  a  single  sample  that  the  generalized  likelihood 
ratio  (GLR)  detector  for  signals  in  bandpass  Gaussian-Gaussian. mixture 
noise  achieves  an  improvement  in  detection  performance  for  low  SNR  over  that 
of  the  conventional  (quadrature  or  square-law  envelope)  detector. 

Now  we  consider  the  detector  performances  achieved  when  the 
number  of  samples  (K)  is  greater  than  one.  Several  issues  become  significant 
for  K  >  1  which  do  not  exist  for  K  =  1.  First,  it  becomes  necessary  to 
distinguish  between  signal  models  that  we  have  called  Type  1  (slowly  varying 
phase)  and  Type  2  (independent  phase  samples).  In  both  cases  the  detection 
is  incoherent  in  the  sense  that  the  instantaneous  signal  phase  is  unknown 
(assumed  to  be  uniformly  distributed  on  (0,  2tt ) ;  the  difference  between  the 
two  models  is  in  the  assumed  rate  of  phase  fluctuation  (bandwidth)  over 
the  observation. time  during  which  the  K  samples  are  taken. 

The  second  issue  which  becomes  significant  for  K  >  1  is  the 
assumption  concerning  the  joint  distribution  of  the  K  noise  samples.  Because 
of  the  Gaussian-Gaussian  mixture  probability  density  function  (pdf)  assumed, 
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(3.4-1) 


the  quadrature  samples  jn^,  j  are  uncorrelated  or  linearly  independent 
regardless  of  the  assumed  pdf  for  the  variance  mixture,  of  the  general  form 


M 


M 


(3.4-2) 


As  discussed  in  Section  2.1.2,  the  number  of  terms  in  the  pdf  p  2(<i),  and 
their  weights  |Cm[,  depend  on  the  assumptions  about  the  independence  or 
dependence  of  the  non-Gaussian  quadrature  sample  pairs. 


3.4.1  Forms  of  the  GLR  for  Multiple  Samples 


For  K  >  1  samples  of  the  quadrature  components  of  the  received 


waveform,  the  generalized  likelihood  ratio  (GLR)  for  the  various  assumptions 
about  the  signal  phase  and  the  variance  multiplier  take  the  various  forms 
discussed  previously  in  Section  3.1.  We  observe  that,  in  general,  the 
independent  signal  phases  (Type  2)  plus  independent  noise  variance  case 
has  the  most  convenient  form,  sinch  the  detection  probability  admits  the 
development 


Pr 


{Ak  (— )  >  "} 


(CASE  II) 


(3.4-3) 
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Thus  an  equivalent  detector  for  this  case  when  there  are  multiple  samples 
is  a  simple  extension  of  the  detector  Ai(r^)  for  a  single  sample. 

The  second  most  convenient  form  is  that  for  independent  phase 
and  equal  v£,  since  it  can  be  constructed  using  the  statistics 

zi(L)  =  Y.  K  +  Bk)  =  IRI’  (3.4-4) 

k  k 

2m  2. 

the  sum  of  samples  of  the  squared  envelope,  and,  for  v  =  1  or  V  , 


=  -Kp/v2  +  pZl(r)/2v4af,  p  small  (3.4-5b) 

~  -Kp/v2  +  Rk;  P  lar9e-  (3.4-5c) 

V 

The  GLR  for  this  case  is  formed  from  zx  and  z2  as 

A(r)  =  {l-W[Zl(r)]}e"Z2(-’  1}  +  W[zr(r)]  e"Z2(r’  V2)  (3.4-6a) 
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eV~^  exp { -z i/2af  V2} 

where  W(z  )  =  - 

( 1-e)  exp{-z i/2a2}  +  eV-^  exp{-z \/2a\  V2} 


(3.4-6b) 


Similarly,  the  detector  for  constant  signal  phase  (Type  1 
signal)  and  slowly-varying  noise  variance  (equal  (v^l)  can  be  reasonably 
implemented  using  the  statistic  Zi(r_)  and  the  statistic 

^>=(l«k)2  ♦(l^*  <3. 

since 


(3.4-8) 


The  GLR  for  this  case  is  formed  from  z 3  and  z3  as 


A(r)  ={1-W[z1(r)]}  e”Kp  I0(  /?pz3(r)/of ) 

+  W[z3 (r) ]  e“Kp/,V  I0  (/2pz3(r.)/a^V4) , 


(3.4-9) 


where  W  [z1(r_)]  is  given  by  (3.4-6b). 

The  most  complicated  detector  is  that  for  equal  phase  and 

K 

independent  v£,  which  requires  generating  Zi(_r)  and  M=2  statistics  of  the 
form 
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z 


4m 


;  m=  1,  2 


...M; 


(3.4-10) 


as  noted  previously  in  Section  3.1. 


3.4.2  Numerical  Results  for  Independent  Samples 

For  independent  samples  of  the  envelope  of  the  signal  plus 
bandpass  Gaussian-Gaussian  mixture  noise,  the  performance  of  the  optimum 
detector  (3.4-3)  for  Type  2  signal  can  be  computed  using  the  numerical 
convolution  technique  shown  in  Appendix  3B.  Since  the  anticipated  effect 
of  using  multiple  samples  is  to  increase  the  detection  probability,  it  is 
sufficient  to  consider  but  a  few  cases  to  verify  this  effect. 

Figure  3.4-1  gives  the  detection  probability  for  the  optimum 
detector  for  Pp^=10-1,  e=0.1,  V2=10,  and  the  number  of  independent  samples, 
K,  equal  to  1,  2,  5,  10,  and  20.  The  accumulated  signal  energy  allows 
90  percent  detection  with  about  7.3  dB  less  SNR  for  K=10  than  for  a  single 
sampl  e . 

Figure  3.4-2  is  similar  to  Figure  3.4-1,  except  that  the 
variance  multiplier  is  increased  to  V2=100.  For  this  case,  about  7.5  dB  in 
detectability  is  gained  by  observing  K=10  samples. 
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PROBABILITY  OF  DETECTION 


Figure  3.4-1.  Performance  of  optimum  detector  for  signals 
in  bandpass  Gaussian-Gaussian  mixture  noise 
(e=0.1 ,  V2=10)  for  P^IO"1  and  multiple, 
independent  samples. 
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PROUAGIUTY  OF  DETECTION 


Figure  3.4-2.  Performance  of  optimum  detector  for  signals 
in  bandpass  Gaussian-Gaussian  mixture  noise 
(e=0.1,  V2=100)  for  PF-=10- 1  and  multiple, 
independent  samples. 
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4.0  SUBOPTIMUM  DETECTORS 

4.1  APPROXIMATIONS  TO  THE  OPTIMUM  DETECTOR 

Exact  implementation  of  the  generalized  likelihood  ratio  (6LR) 
for  bandpass  Gaussian  mixture  noise,  given  by  equation  3.3-1,  requires 
generating  the  function 


f  (a ,  b)  =  e“b  I o ( /2aF) ,  (4.1-1) 

where  the  argument  "a"  is  proportional  to  the  squared  envelope  of  the  received 
waveform  and  "b"  is  proportional  to  the  a  priori  SNR,  p0.  (See  Figure  3.2-1). 
We  now  consider  simplifications  which  result  from  approximating  f(a,  b)  and 
the  impact  on  detector  performance 


4.1.1 


Large  SNR  Approximations 

For  large  values  of  SNR  (b»l),  we  have 

-b+/2aF 


f (a,  b) 


,  b>>l. 


(4.1-2) 


'  2tt  /2ab 

with  the  resulting  GLR  approximation  (x  =  R2/a^) 


.  /...  \  const 
A  (x ,  p )  ~ 

r  WT 


(1-e)  exp 

(^ -  >^)2 

,  e  ) 

4-  —  flvn  { _ 

-  /SD2 

2 

“2V7 

(1-e)  exp|-x/2j  +  exp  j-x/2V2j 


(4.1-3) 


This  approximation  is  expected  to  be  valid  [23]  for  /2ab  >  3.75,  or 
p  >  7V4  /x. 
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Examination  of  the  numerator  of  (4.1-3)  suggests  that  for  large 
SNR  the  GLR  will  exhibit  a  peak  near  the  normalized  squared  envelope  value 
x  =  2p0.  This  phenomenon  is  confirmed  for  p0  =  20  dB  and  V2  =  100  or  1000 
by  the  previous  Figures  3.3-2,  3.3-3,  and  3.3-6.  Thus  we  observe  that  in 
general  the  GLR  characteristic  acts  as  a  "window",  permitting  high  output  values 
only  in  the  vicinity  of  input  values  near  the  anticipated  SNR  (if  the  signal 
is  present),  and  suppressing  the  output  for  higher  input  values  which  are  more 
likely  to  be  due  to  noise. 

Implementation  of  the  high  SNR  approximation  to  the  GLR  is  only 
slightly  less  complicated  than  that  of  the  GLR  itself.  More  important,  this 
approximation  is  still  parametric  in  SNR,  offering  no  advantage  in  terms  of 
a  priori  information  requirements. 

4.1.2  Small  SNR  Approximation  and  Locally  Optimum  Detector. 


For  small  values  of  SNR  (b<<l),  we  have 

f(a,  b)  «  1  +  b  -  lj  ,  b<<l, 


(4.1-4) 


with  the  resulting  GLR  approximation  (xsR2/af) 

(1-e)  e”x/2  V**(x-2)  +  jz  e~x/2yZ  (x-2V2) 


Ar(x;p)  »  1  +  — 
r  2V4 


•  (4.1-5) 


(1-e)  e  x^2  +  irr  e 


■x/2V2 


This  approximation  is  expected  to  be  valid  for  /2ab  <  3.75  or  p  <  7/x. 

For  detection  purposes,  we  may  ignore  the  additive  constant  in  (4.1-5)  and 
also  the  constant  factor  p/2,  and  use  the  weak  signal  "locally  optimum 
detector"  (LOD), 
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Z(x)  =  [l-VJ(x )]  (x-2)  +  W(x)  (x-2V2)/V4  (4. l-6a) 


where 


(1-e) 


e 

w 


e~x/2V2 


(4. l-6b) 


The  form  of  (4.1-6)  is  rather  easily  interpreted  as  a  combination 
of  LOD's  for  the  Gaussian  case,  Zg(x)  =  x-2 o2/o\  .  This  interpretation  can 
also  be  understood  from  the  plots  of  (4.1-6)  given  in  Figures  4.1-1  to  4.1-3 
for  various  combinations  of  e  and  V2. 

We  note  that  the  LOD  is  not  parametric  in  the  SNR,  although  it 
still  requires  a  priori  information  in  the  form  of  e  and  V2,  of  course. 

It  is  often  argued  that  the  weak  signal  case  is  the  most  interesting  one. 
However,  historically  the  LOD  has  been  studied  primarily  because  it  usually 
involves  a  simpler  detector  structure  than  the  GLR  and  therefore  is  more 
amenable  to  analysis. 

For  Gaussian  and  other  monotonic  GLR's,  the  LOD  performs  well  at 
high  SNR.  As  will  be  demonstrated  below,  this  is  not  the  case  for  the 
Gaussian-Gaussian  mixture  LOD. 


4.1.3  Performance  of  the  Locally  Optimum  Detector 

Computations  of  the  false  alarm  and  detection  probabilities  for 
the  bandpass  Gaussian-Gaussian  mixture  noise  locally  optimum  detector  proceed 
in  a  way  similar  that  described  for  the  GLR  in  Section  3.3,  with  the  important 
exception  that  the  false  alarm  thresholds  no  longer  depend  on  an  a  priori 
or  assumed  value  of  SNR. 
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Figure  4.1-1  Locally  optimum  detector  characteristic  for  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.01)  for  several 
values  of  variance  ratio 
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Figure  4.1-2  Locally  optimum  detector  characteristic  for  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.1)  for  several 
values  of  variance  ratio. 
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NORMALIZED  ENVELOPE  SQUARED 

Figure  4.1-3  Locally  optimum  detector  characteristic  for  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.5)  for  several 
values  of  variance  ratio. 
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4. 1.3.1  False  Alarm  Probability 

Single-sample  false  alarm  probabilities  for  the  LOD  given  by 
(4.1-6)  are  plotted  in  Figures  4.1-4  to  4.1-6  as  a  function  of  the  first 
threshold  on  the  normalized  squared  envelope,  m/af  •  It  may  be  observed 
from  these  figures  that  the  false  alarm  probability  is  very  sensitive  to 
variation  in  the  threshold  for  PpA  <  e,  because  at  approximately  this  point 
the  threshold  x  on  the  LOD  starts  being  above  the  local  maximum.  For 
example,  in  the  plot  of  the  LOD  characteristic  for  e  =  0.01  (Figure  4.1-1), 
the  peaks  for  V2  =  10,  100,  and  1000  occur  at  normalized  squared  envelope 
values  of  12.5,  15.2,  and  19,  respectively.  In  Figure  4.1-4  we  see  that  the 
PpA  curves  for  these  cases  have  very  steep  slopes  in  the  vicinity  of  these 
values. 

A  similar  effect  is  seen  for  V2  =  100  and  1000  in  Figure  4.1-6, 
in  which  the  PpA  is  seen  to  decrease  suddenly  near  the  threshold  value  of  2; 
this  behavior  corresponds  to  the  X  threshold  on  the  LOD  in  Figure  4.1-3 
starting  to  rise  above  the  local  minimum  of  the  LOD. 

False  alarm  thresholds  for  the  LOD  are  given  in  Table  4.1-1 
for  various  values  of  e,  V2,  and  PpA- 

4. 1.3. 2  Detection  probability  for  a  single  sample. 

Using  the  false  alarm  thresholds  in  Table  4.1-1,  the  detection 
probability  for  the  locally  optimum  weak  signal  detector  in  bandpass 
Gaussian-Gaussian  mixture  noise  was  computed  according  to 

Pq(p )  =  Pi (p ;  m)  -  ?i(p;  ^2)  +  pi(p‘»  ^3),  (4.1-7) 
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Figure  4.1-4.  False  alarm  probability  vs.  first  threshold  for  locally 
optimum  detector  in  bandpass  Gaussian-Gaussian  mixture 
noise  (e  =  0.01). 
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Table  4.1-1. 


False  alarm  thresholds  for  locally  optimum  detector  in 
bandpass  Gaussian-Gaussian  mixture  noise. 


EPSILON 

V2 

TARGET  PFA 

ETA1 

ETA2 

0.01 

1.0 

1.000D-01 

4.60517D+00 

0.01 

1.0 

1.0000-02 

9 . 21034D+00 

0.01 

1.0 

1.000D-03 

1 . 38155D+01 

0.01 

1.0 

1.000D-04 

1.84207D+01 

0.01 

10.0 

1.000D-01 

4.66694D+00 

1.909520+01 

0.01 

10.0 

1.000D-02 

9.40157D+00 

1.504550+01 

0.01 

10.0 

1. 0000-03 

1 . 20191D+01 

1.280400+01 

0.01 

10.0 

1.000D-04 

1 . 23762D+01 

1.245520+01 

0.01 

100.0 

1.0000-01 

4.60190D+00 

2.24671D+01 

0.01 

100.0 

1.000D-02 

9.27215D+00 

1.92562D+01 

0.01 

100.0 

1.0000-03 

1 . 35694D+01 

1.651750+01 

0.01 

100.0 

1.0000-04 

1.49714D+01 

1 . 53022D+01 

0.01 

1000.0 

1.000D-01 

4.58750D+00 

2.73836D+01 

0.01 

1000.0 

1.0000-02 

9.204850+00 

2.45440D+01 

0.01 

1000.0 

1.0000-03 

1 . 38548 D+01 

2.24278D+01 

0.01 

1000.0 

1.000D-04 

1 . 78277D+01 

2.003690+01 

0.10 

1.0 

1.000D-01 

4.60517D+00 

0.10 

1.0 

1.000D-02 

9.21034D+00 

0.10 

1.0 

1.000D-03 

1 . 38155D+01 

0.10 

1.0 

1.000D-04 

1.84207D+01 

- - 

0.10 

10.0 

1.000D-01 

4.90189D+00 

1 . 23322D+01 

0.10 

10.0 

1.0000-02 

7 . 98925D+00 

9.01454D+00 

0.10 

10.0 

1.000D-03 

8.44989D+00 

8.55314D+00 

0.10 

10 » 0 

1.000D-04 

8.49635D+00 

8.50667D+00 

0.10 

100.0 

1.0000-01 

4.51264D+00 

1.69902D+01 

0.10 

100.0 

1.000D-02 

9 . 1 1 103D+00 

1.304710+01 

0.10 

100.0 

1.0000-03 

1 . 0947  5D+01 

1 . 14108D+01 

0.10 

100.0 

1.000D-04 

1 . 115740+01 

1.120390+01 

0.10 

1000.0 

1.000D-01 

4.412160+00 

2 . 22325D+01 

0.10 

1000.0 

1.000D-02 

9  08521D+00 

1.89195D+01 

0.10 

1000.0 

1.000D-03 

1.33478D+01 

1 . 61705D+01 

0.10 

1000.0 

1.000D-04 

1.46877D+01 

1 . 50013D+01 

0.50 

1.0 

1.000D-01 

4.605170+00 

0.50 

1.0 

1.000D-02 

9 . 21034D+00 

0.50 

1.0 

1.000D-03 

1 . 38155D+01 

0.50 

1.0 

1.000D-04 

1.842070+01 

0.50 

10.0 

1.0000-01 

4.39572D+00 

7.425940+00 

0.50 

10.0 

1.0000-02 

5.62999D+00 

5.933800+00 

0.50 

10.0 

1.000D-03 

5.76669D+00 

5.794580+00 

0.50 

10.0 

1.000D-04 

5.78044D+00 

5 . 78062D+00 

0.50 

100.0 

1.000D-01 

3.693220+00 

1.28308D+01 

0.50 

100.0 

1.000D-02 

7.255220+00 

8.64562D+00 

0.50 

100.0 

1.0000-03 

7.87922D+00 

8.020080+00 

0.50 

100.0 

1.000D-04 

7.942600+00 

7.95668D+00 

0.50 

1000.0 

1.000D-01 

3.29672D+00 

1.88056D+01 

0.50 

1000.0 

1.000D-02 

8.04259D+00 

1.401850+01 

0.50 

1000.0 

1.000D-03 

1.08155D+01 

1 . 16785D+01 

0.50 

1000.0 

1.0000-04 

1 . 12088D+01 

*  1.129590+01 

ETA3 


2 .843860+02 
7 . 11396D+02 
8. 360350+02 
8.387230+02 
2.61932D+04 
7.22034D+04 
1.06968D+05 
1 . 11320D+05 
2.58924D+06 
7 . 19963D+06 
1 . 17367D+07 
1.47332D+07 


2.82218D+02 
4.42917D+02 
4.46757D+02 
4.467960+02 
2.506650+04 
6.476830+04 
7 . 17003D+04 
7 . 18029D+04 
2.4 1 173D+06 
7.01433D+06 
1.04380D+07 
1.08404D+07 


1.52506D+02 
1 . 72626D+02 
1.72842D+02 
1.728440+02 
1.61288D+04 
3.87099D+04 
3 . 94852D+04 
3.94933D+04 
1.29202D+06 
5.72649D+06 
7.216060+06 
7.252020+06 
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where  P]^  )  is  given  by  (3.3-15).  The  results  are  shown  in  Figures  4.1-7 
to  4.1-10. 

The  cases  of  e  =  0.01  and  V2  =  1,  10,  100,  and  1000  are  shown 
in  Figure  4.1-7  for  false  alarm  probabilities  of  10-1,  10"2,  and  10"3. 

The  V2  =  1  cases  represent  the  performance  of  the  Gaussian  detector  in 
Gaussian  noise.  We  observe  that  the  LOD's  performance  closely  follows 
that  of  the  V2  =  1  case  for  low  SNR  (when  the  signal  actually  is  weak,  as 
assumed)  and  for  higher  Pp^.  As  the  required  Pp^  is  decreased,  the  LOD 
detection  probability  departs  more  from  the  V2  =  1  case.  The  most  sig¬ 
nificant  fact  made  apparent  by  the  figure  is  the  severe  degradation  in 
performance  when  the  SNR  is  greater  than  6  dB  or  7  dB;  P^  actually  falls 
below  Pp^  before  rising  again  to  unity  for  very  high  SNR. 

A  comparison  of  Figure  4.1-7  with  Figure  3.3-13  reveals  that 
the  LOD  detection  performance  is  quite  close  to  the  optimum  for  SNR  values 
less  than  zero  dB. 

It  is  interesting  to  note  that  the  LOD  PD  is  generally  higher 
for  higher  values  of  the  variance  ratio,  V2  when  V2  is  greater  than  10. 
This  behavior  is  due  to  the  wider  "window"  or  acceptance  region  for 
higher  V2  which  was  apparent  from  the  plots  of  the  LOD  characteristic. 

Increasing  e  to  0.1  yields  the  curves  shown  in  Figure  4.1-8, 
and  for  e  =  0.5  the  Pp  curves  of  Figure  4.1-9  are  produced.  The  P^  per¬ 
formance  is  seen  to  deteriorate  with  increasing  values  of  the  mixture 
parameter,  c,  except  at  high  SNR's.  The  fact  that  eventually,  for  very 
high  SNR,  the  P^  approaches  unity  is  demonstrated  by  Figure  4.1-10. 
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SIGNAL  TO  NOISE  RATIO  CdB) 


Figure  4.1-8.  Receiver  operating  characteristics  for  locally  optimum 
detector  in  bandpass  Gaussian-Gaussian  mixture  noise 
(e  =  0.1). 


133 


DETECTION  PROBABILITY 


134 


DETECTION  PROBABILITY 


Figure  4.1-10.  Receiver  operating  characteristics  for  locally 

optimum  detector  in  bandpass  Gaussian-Gaussian  mixture 
noise  for  e  =  0.1,  V2  =  10,  and  Pp^  =  10_1. 
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4.1.4  Performance  of  the  Weak  Signal  Locally  Optimum  Detector  for 
Multiple  Samples 

The  LOD  statistic  Z(x)  given  by  (4.1-6)  pertains  to  a  single, 
normalized  squared  envelope  sample,  x=R 2/of.  For  independent  multiple 
samples  it  is  easy  to  show  that  the  locally  optimum  detector  is 

K 

zK(x)  =  £  Z(xk).  (4.1-8) 

k=l 

Since  Z(x)  is  nonlinear,  the  distribution  of  Zk(x.)  is  very  difficult  to 
obtain  analytically,  if  possible.  Instead,  the  numerical  technique  given 
in  Appendix  3B  was  used  to  evaluate  the  detection  performance  of  Z^U)  • 

Figures  4.1-11  and  4.1-12  show  the  performance  of  the  statistic 
Z^  for  bandpass  Gaussian-Gaussian  mixture  noise  (e=0.1  and  V2=  100)  for 
P^IO"1  and  10"2,  respectively.  The  numerical  technique  used  32  levels  to 
represent  the  pdf  for  a  single  sample.  We  observe  from  these  figures  that 
the  detection  performance  is  much  improved  at  low  signal  levels  when  the 
number  of  independent  samples  is  increased.  The  same  "dip"  in  performance 
for  moderate  SNR's  is  observed  as  for  one  sample,  however.  Thus,  unless 
one  takes  the  view  that  we  are  interested  only  in  signals  with  SNR  <  6  dB, 
say,  then  we  cannot  be  satisfied  with  using  the  locally  optimum  detector. 

In  the  next  sections  we  consider  other,  suboptimum  detectors 
which  do  not  perform  as  well  as  the  LOD  for  small  SNR's,  but  on  the 
other  hand  perform  better  uniformly  as  the  SNR  increases. 
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SNR  (dB) 


Figure  4.1-11.  Multiple-sample  performance  of  locally  optimum 

detector  for  bandpass  signals  in  Gaussian-Gaussian 
mixture  noise,  for  e=0.1,  V2=100  and  P_.=10_1. 


137 


PflOB ABILITY  OF  DETECTION 


Figure  4.1-12.  Multiple-sample  performance  of  locally  optimum 

detector  for  bandpass  signals  in  Gaussian-Gaussian 
mixture  noise,  for  e=0.1,  V2=10Q,  and  PF.=10-2. 
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4.2 


SERIAL  NORMALIZATION  DETECTORS 


Under  various  assumptions  concerning  the  time  dependency  among 


samples  of  the  bandpass  Gaussian-Gaussian  mixture  process,  we  can  develop 
suboptimal  detector  approaches  and  evaluate  them.  In  this  section,  we 
consider  detectors  based  on  an  assumption  of  slowly-varying  noise  variance. 
In  Section  4.3  this  assumption  is  relaxed  but  it  is  assumed  that  several 
channels  of  data  have  noise  variances  which  change  in  the  same  way  with 
time. 

4.2.1  Conceptual  formulation 


If  K  samples  of  a  Gaussian  process  with  unknown,  fixed  variance 


are  to  be  used  to  test  the  hypothesis  H0:m=0  against  the  composite  alter¬ 
native  Hoim^O,  where  m  is  an  unknown  but  fixed  mean  value,  then  an 
appropriate  test  statistic  [27]  is  the  ratio  (Student's  t-statistic) 


x 


(4.2-1) 


where  7  is  the  sample  mean  of  the  K  samples  (xx,  x2,  ...,  x^,). 

For  the  bandpass  case,  this  test  statistic  generalizes  to  the 
quantity  [25] 


xJ  = 


(4.2-2) 


rA:’  -s 


139 


J.  S.  LEE  ASSOCIATES,  INC. 


where  the  {x^,  x^}  are  the  in-phase  and  quadrature  samples  of  the 
data. 

For  example,  for  K=2  samples  (the  minimum  number),  the  test 
statistic  can  be  expressed  as 


(x„  +x„  )2  +  (x  +x  )2 


T2(x„  ,  x„  ,  xc  ,  xe  )  = 


Si  s2 


Cl'  C2’  $1 


(x  -X  )2  +  (x  -  X  )2 
C1  C2  S1  S2 


(4.2-3) 


This  statistic  may  be  interpreted  as  an  estimate  of  the  SNR. 
Detection  of  a  signal  is  said  to  have  occurred  when  this  estimate  exceeds 
a  threshold. 


4.2.2.  Analysis  of  the  detector  performance 

It  can  be  shown  that  under  the  assumption  that  the  K  successive 
pairs  of  quadrature  samples  are  indeed  jointly  Gaussian,  conditioned  on 
a  random,  unknown  variance  which  is  either  a2  or  a2  for  all  K  samples, 
then  T^,  is  distributed  as  a  noncentral  F-statistic: 


Tk  ^  F(2,  2K-2,  X=2KP). 


(4.2-4) 


Because  of  the  randomness  of  p,  we  can  express  the  distribution 

of  T.,  in  this  "slowly  varying"  bandpass  Gaussian-Gaussian  mixture  noise 
l\ 

unconditionally  by 


T 


K 


/ 

F(2,  2K-2;  2KP)  with  probability  (1-e) 

< 

F(2,  2K-2;  2Kp/V2)  with  probability  c; 


(4.2-5) 
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since  the  ratio  form  of  the  statistic  in  effect  cancels  out  the  variance, 
the  random  distributional  parameter  becomes  the  SNR. 

4. 2. 2.1  False  Alarm  Probability 

Under  the  hypothesis  H0:p=0,  the  distribution  of  is  simply 

Tk  "  F(2,  2K-2;  0),  (4.2-6) 

or  central-F.  Thus  the  false  alarm  probability  does  not  depend  on  either 
e  or  V2,  the  parameters  of  the  mixture: 


pfa"  Pr  {Tk  >  ^p-0> 
=  I^K-1,  1] 

with 

,  A  K-l  . 

K-l+n  . 


where  I^[a,  b]  is  the  (normal i zed) incompl ete  Beta  function: 


I5(a,  b)  =  /  dt  ta_1  (l-t)5"1  x  r(a+  b)/r(a)r(b) 


'o 


b-1 

a  V'  (1-g) 


■  r  I 

r=0 


,  (a)r>  b  an  integer. 


Therefore  the  false  alarm  probability  has  the  simple  form 


PFA  =  sK' 


■i  rjwJ 

[K-l+nJ 


K-l 


(4.2-7a) 


(4 . 2 -7  b ) 


(4.2-7c) 

(4.2-7b) 


(4.2-8) 
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this  is  easily  inverted  to  yield  the  threshold  expression 

n0  =  ( K-l )  Uo1-!]  =  (K-U  [pFA"1/(K'1)  -11.  (4.2-9) 

For  example,  for  K=2  samples, 

n0  ^o1  -1  »  Pp^1  >1  •  (4-2.10) 

Because  no  and  £0  are  equivalent  thresholds,  we  can  express 
our  results  in  terms  of  either  one. 

4. 2. 2. 2  Detection  Probability 

From  (4.2-5)  it  follows  that  the  detection  probability  is 

PD  =  Pr  {T^  >  n  |  P5*0} 

=  (1-e)  Pr  {F(2,  2K-2;  2Kp)  >  n) 

+  e  Pr {F(2,  2K-2;  2Kp/V2)  >  n) 

=  (1-e)  Pi(ri;2Kp)  +  ePjtn;  2Kp/V2)  (4.2-11) 

where  we  define 

Pi(n;  X)  =  Pr {F(2,  2K-2;  X)  >  nK  (4.2-12) 
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Various  expressions  can  be  developed  for  Pi  .  Among  the  most  convenient 
ones  are 

Pi(n;  x)  =  e~x/Z 

=  (K-l) 

with  5  as  given  previously  by  (4.2-7b).  The  expression  given  by  (4.2-13a) 
is  based  on  the  noncentral-F  statistic  probability  integral  in  [23]; 
(4.2-13b)  follows  from  the  development  shown  in  Appendix  4A. 

4.2.3.  Numerical  results 

Figures  4.2-1  through  4.2-9  give  the  detection  performance  of  a 
serial  normalization  detector  for  bandpass  Gaussian-Gaussian  mixture 
noise  and  variance  ratios  of  1,  10,  and  100.  The  figures  use  the 
following  parameters: 


I 

r=0 

K-l 


I?  (K-l,  r+1) , 


(4.2-13a) 


i  (?) 

r=0 


XV 


r! 


dt  tK-2  (l-t)r  e"U/2’ 


(4.2-13b) 


e,  mixture  parameter 


0.01 

0.1 

0.5 

K,  number 

2 

Fig  4.2-1 

Fig  4.2-2 

Fig  4.2-3 

of  samples 

3 

Fig  4.2-4 

Fig  4.2-5 

Fig  4.2-6 

4 

Fig  4.2-7 

Fig  4.2-8 

Fig  4.2-9 
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Although  the  number  of  samples  is  small,  we  are  reluctant  to 
use  larger  numbers  because  the  assumption  is  made  that  the  noise  is 
stationary  Gaussian  for  at  least  K  samples. 

From  these  figures,  we  observe  that  the  parameter  V2  does  not 
affect  the  detection  performance  significantly  for  small  values  of  e. 
Increasing  the  number  of  samples  does  improve  the  performance,  but  the 
performance  is  subject  to  a  (temporary)  levelling-off,  for  high  values 
of  V2,  which  becomes  more  evident  as  K  increases.  For  example,  in 
Figures  4.2-7  to  4.2-9,  we  observe  that  the  detection  probability  rises 
to  a  value  slightly  greater  than  (1-e)  before  levelling  off.  This 
behavior  is  attributable  directly  to  the  two  terms  in  the  expression 
(4.2-11),  which  rise  at  the  same  rate  with  SNR,  but  at  a  10  log^Q  (V2)  dB 
separation, and  with  different  rates. 


144 


DETECTION  PROBABILITY 


EPSILON-O.  01 
K=2 


Figure  4.2-1  Performance  of  serial  normalization  detector  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.01)  for  two  time 
samples 
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-7.00  T 00  TOO  TOO  13.00  78.00  23. 00 

SIGNAL  TO  NOISE  RATIO  CdB) 

Figure  4.2-2  Performance  of  serial  normalization  detector  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.1)  for  two  time 
samples 
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SIGNAL  TO  NOISE  RATIO  (dB) 

Figure  4.2-3  Performance  of  serial  normalization  detector  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.5)  for  two  time 
samples 


147 


DETECTION  PROBABILITY 


SIGNAL  TO  NOISE  RATIO  (dB) 

Figure  4.2-4  Performance  of  serial  normalization  detector  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.01)  for  three 
time  samples. 
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Figure  4.2-5  Performance  of  serial  normalization  detector  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.1)  for  three 
time  samples 
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Q.0001— j - - 1 - 1 - 1 - 1  i  n 

-7.00  -2.00  3.00  8.00  13.  00  18.00  23.00 

SIGNAL  TO  NOISE  RATIO  (dB) 

Figure  4.2-6  Performance  of  serial  normalization  detector  in 

bandpass  Gaussian-Gaussian  mixture  noise  (e  =  0.5)  for 
three  time  samples 
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Figure  4.2-7  Performance  of  serial  normalization  detector  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.01)  for  four 
time  samples 
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SIGNAL  TO  NOISE  RATIO  (dB) 

Figure  4.2-8  Performance  of  serial  normalization  detector  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.1)  for  four  time 
samples 
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SIGNAL  TO  NOISE  RATIO  (dB) 


Figure  4.2-9  Performance  of  serial  normalization  detector  in 

bandpass  Gaussian-Gaussian  mixture  noise  (e  =  0.5)  for 
four  time  samples 
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4.3  PARALLEL  NORMALIZATION  DETECTORS 

The  Gaussian-Gaussian  mixture  noise  we  have  been  treating  as 
a  particular  case  of  non-Gaussian  noise  can  be  considered  as  arising  from 
random  time  variation  in  the  power  or  variance  of  a  Gaussian  noise  process. 

In  this  view,  the  noise  variance  of  is  that  of  a  stationary,  background 
Gaussian  noise  process,  and  >  o\  is  the  combined  variance  of  the 
background  process  plus  a  "switched  burst"  [26]  of  higher  variance  noise. 

The  mixture  parameter  e  then  corresponds  to  the  probability  of  having 
received  the  total  noise  of  variance  o|,  modelled  in  our  work  as  bandpass 
Gaussian. 

For  weak  signals  in  lowpass,  "switched  burst"  or  Gaussian-Gaussian 
mixture  noise,  it  was  shown  in  [26]  that  a  detector  of  the  form  illustrated 
in  Figure  4.3-1  can  perform  detections  more  reliably  than  a  fixed  detector 
structure  in  high  kurtosis  Arctic  under-ice  noise  data.  In  effect,  given 
K  samples,  the  data  are  classified  into  two  groups  on  the  basis  of  their 
magnitudes;  the  larger  data  samples  are  processed  assuming  the  noise  is 
background  plus  the  higher  variance  noise,  while  the  remainder  of  the  samples 
are  processed  as  if  only  background  noise  is  present,  in  addition  to 
a  possible  signal.  The  detector  nonlinearities  are  the  appropriate 
locally  optimum  detectors. 

The  extension  of  the  detector  concept  of  Figure  4.3-1  to  the 
bandpass  Gaussian-Gaussian  mixture  case  would  take  the  form  shown  in 
Figure  4.3-2.  Since  the  locally  optimum  or  weak  signal  detector  for 
noncoherent  signals  in  bandpass  Gaussian  noise  is  the  square-law  envelope 
detector,  the  switching  operation  selects  whether  normalization  is  per¬ 
formed  by  af  or  by  af  =  a^V2.  In  part  (b)  of  the  figure,  we  show  that 
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smoothing:  =  1  if  more  than  half  of  (j=0,  1,  m)  are  equal  to  one 

t  =  1.282  oj,  Oj  estimated  from  data. 


Figure  4.3-1  Detector  for  weak  signal  in  lowpass,  switched  burst  noise. 


BPF 


smoothing 


Figure  4.3-2  Switched  gain  detector  for  bandpass  Gaussian-Gaussian  mixture  noise. 
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if  the  switching  decision  is  not  smoothed  over  ±m  samples  but  is 
instantaneous,  then  the  resultant  detector  characteristic  is  a  nonlinearity 
which  is  a  simplified  form  of  the  locally  optimum  detector  discussed 
previously. 

4.3.1  Parallel  Normalization  Concept 


In  Section  4.2  we  based  a  detector  on  the  assumption  that 


successive  time  samples  of  the  Gaussian-Gaussian  mixture  process  were 
jointly  Gaussian  with  the  same  variance  for  short  periods  of  time.  Now 
we  consider  a  detector  based  on  the  following  assumptions: 

(a)  one  or  more  "parallel"  channels  are  available,  not  con¬ 
taining  a  signal,  for  a  total  of  M  channels 

(b)  these  channels  contain  bandpass  Gaussian-Gaussian  mixture 
noise,  with  the  joint  pdf  at  the  same  time 


(4.3-1) 


Without  loss  of  generality,  we  assume  that  the  channel  of  interest  is 
Channel  1  (m=l);  (4.3-1)  implies  that,  given  whether  all  M  channels 
have  variance  of  or  of,  they  are  conditionally  independent  bandpass 
Gaussian  processes. 
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Using  these  assumptions,  a  noncoherent  detector  can  be  formulated 
as  follows:  at  time  t^,  the  square-law  envelope  detector  sample  in  the 
channel  to  be  tested  for  a  signal , 

Rik  =  °ik  +  e?k  •  (4-3‘2) 

can  be  normalized  by  the  sum  of  the  samples  in  the  M  channels  to  form  the 
statistic. 


IR.2k 

m=l 


(4.3-3) 


where  the  sum  acts  as  a  quantity  proportional  to  an  estimate  of  the 
noise  variance  common  to  the  channels. 

Detection  then  consists  of  comparing  the  accumulation  of 
K  such  statistics  to  perform  the  test 

K  Hl 

2  =  ^2^  >  threshold.  (4.3-4) 

k=l  H0 


4.3.2  Distribution  of  Test  Statistic 

It  is  relatively  easy  to  show  that  for  K=1  and  M=2,  then  the 

pdf  of  the  test  statistic  given  by  (4.3-4)  is 

/ 


P2(«) 


{ 1— e )  e"p(1_a)  (1  +  Po) 
i  +  e  e  p^  (1  +  pa/V2),  0  1 ", 


(4.3-5) 


0,  elsewhere. 
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Note  that  when  p=0  (H0  is  true),  then  z  is  uniformly  distributed 
between  0  and  1. 

Convolution  may  be  used  to  obtain  results  for  K>1.  For  M>2, 
it  may  be  shown  that 


where 


Pz(a)  =  (1-e)  p)  +  e  PM(a;  p/V2) 


pm(«;  P)  = 


(M-l)  e-fO-0)  (1-„)H-2  (.pa), 


< 


0  < 


a  s  1; 


0,  otherwise. 


(4.3-6a) 


(4.3-6b) 


In  (4.3-6b),  • )  is  the  Laguerre  polynomial  of  order  M-l. 

4. 3. 2.1  False  Alarm  Probability 


For  p=0,  (4.3-7)  reduces  to 

Pz(oc  i  H0)  =  (M-l)  (l-a)M"2  ,  0  $  a  <  1.  (4.3-7) 


Thus  for  one  sample,  the  false  alarm  probability  is 
PpA  =  Pr  {z  >  n  |  p=0} 


/ 

i 


n  <  0 
0  $  n  <  1 
n  >  1 


(4.3-8) 
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For  two  samples  (K=2)  and  two  channels  (M=2),  the  false  alarm 
probability  is 


n  <  0 


P 


FA 


< 


(4.3-9) 


|(2-n)2,  1  <  n  <  2 


We  observe  that  the  false  alarm  probability  does  not  depend  on  the 
Gaussian-Gaussian  parameters  e  or  V2;  thus  the  false  alarm  threshold  can 
be  chosen  without  knowing  or  estimating  them. 

4. 3. 2. 2  Detection  Probability 

Integration  of  (4.3-5)  gives  the  K=M=2  case  of  the  detection 
probability: 


PD  =  Pr  {z  >  n  |  P  t  0} 


-  (l-e)  [1-n  e_p(1"n) 


(4.3-10) 


+  e  [1-n  e 


-p(l-n)/V2 

J  9 


0  <  n  < 1. 


Since  Pp^  =  1-n,  (4.3-10)  can  be  written 


PD  =  1~(1'PFA)  [(l“e)  6 


+  e  e 


(4.3-11) 
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Generally,  for  K>2  and/or  M>2  it  is  very  tedious  although 
straightforward  to  find  the  detection  probability.  Numerical  procedures 
such  as  convolution  and  integration  are  recommended.  For  example,  the 
Pp  for  K=M=2  can  be  expressed  by 


where 


(4.3-12) 


f(a;  p)  =  (1-e)2  g(a"»  p,  p) 


+  2e(l-e)  g(a',  p,  p/V2) 


+  e2  g(a;  p/V2,  p/V2) 


(4.3-13) 


and  for  1  s  a  <  2, 


g(a;  px,  p2)  =  e  ^p+pa[(l+p)2  (2-a)  -p(l+p)(2-a)2 


(P1-P2)3 


+  p2(2-a)3/6],  Pl=  p2  =  p  ;  (4.3-14) 

p2(a-2)  [_2pip2  +  p^(Pl-  p2) 

+p2(Pl-  p2)(pf-Pi  p2  -  p2) (a-1 ) ) 


+  ePl^a  ^  [ 2p xp2  +  p^(pi-  P2) 


+  p^pi-  p2)(p^-  P1P2-  Pi )  (a_l)  ] 

(4.3-15) 
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4.3.3  Numerical  Results 


The 

detection  in 
to  4.3-8  for 
used  are  the 


performance  of  a  two-channel  (M=2)  parallel  normalization 
bandpass  Gaussian-Gaussian  noise  is  shown  in  Figures  4.3-3 
variance  ratios  of  V2=l,  10,  and  100.  The  other  parameters 
following: 


e  =  mixture  parameter 


0.01 

0.1 

0.5 

K=  #  of  time 

1- 

Fig  4.3-3 

Fig  4.3-4 

Fig  4.3-5 

samples 

2 

Fig  4.3-6 

Fig  4.3-7 

Fig  4.3-8 

The  behavior  of  this  detector  with  the  various  parameters  is 
very  similar  to  that  of  the  serial  normalization  detector  treated  in 
Section  4.2: 


(a)  V2  has  little  influence  for  small  e. 

(b)  Pp  uniformly  increasing  with  SNR. 

(c)  Pp  improving  for  K  increasing. 

Considering  that  the  assumptions  in  the  analysis  of  this  detector 
are  probably  less  restrictive  than  those  of  the  serial  normalization 
detector  (depending  upon  what  constitutes  a  "parallel"  channel),  we  observe 
that  reasonable  comparable  performances  are  achieved  and  that  therefore 
parallel  normalization  techniques  may  be  preferred  because  of  their  simpler 
impl ementation. 
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EPSILON  =  O.Ol 


SIGNAL  TO  NOISE  RATIO  (dB) 

Figure  4.3-3  Performance  of  two-channel  parallel  normalization 

detector  in  bandpass  Gaussian-Gaussian  mixture  noise 
(e  =  0.01)  for  a  single  time  sample. 
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SIGNAL  TO  NOISE  RATIO  CdB) 

Figure  4.3-4  Performance  of  two-channel  parallel  normalization 

detector  in  bandpass  Gaussian-Gaussian  mixture  noise 
(e  =  0.1)  for  a  single  time  sample. 
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Figure  4.3-5  Performance  of  two-channel  parallel  normalization 
detector  in  bandpass  Gaussian-Gaussian  mixture 
noise  (e  =  0.5)  for  a  single  time  sample. 
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SIGNAL  TO  NOISE  RATIO  (dB) 

Figure  4.3-6  Performance  of  two-channel  parallel  normalization 

detector  in  bandpass  Gaussian-Gaussian  mixture  noise 
(e  =  0.01)  for  two  time  samples. 
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EPSILON  =  0. 10 
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Figure  4.3-7  Performance  of  two-channel  parallel  normalization 
detector  in  bandpass  Gaussian-Gaussian  mixture 
noise  (e  =  0.1)  for  two  time  samples. 
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SIGNAL  TO  NOISE  RATIO  (dB) 

Figure  4.3-8  Performance  of  two-channel  parallel  normalization  detector 
in  bandpass  Gaussian-Gaussian  mixture  noise  (e  =  0.5)  for 
two  time  samples. 
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4.4  COMPARISON  OF  DETECTOR  PERFORMANCES 

Figures  4.4-1  and  4.4-2  provide  graphical  comparisons  of  the 
performances  of  the  weak  signal  LOD,  the  serial  normalization  detector, 
and  the  parallel  normalization  detector.  It  is  seen  that  the  LOD  is  the 
best  choice  when  it  can  be  assumed  that  the  signal  is  weak  (SNR  <  5dB) , 
but  it  does  not  perform  well  if  the  signal  is  strong. 

For  strong  signals  both  the  serial  normalization  detector 
and  the  parallel  normalization  detector  are  superior  to  the  LOD,  with  the 
parallel  one  tending  to  be  more  effective  than  the  serial  one.  In 
addition,  as  discussed  previously,  these  detectors  are  attractive  in  that 
they  can  be  designed  to  maintain  a  given  false  alarm  probability  without 
a  priori  information  on  the  noise  distribution  parameters. 
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Figure  4.4-1 


SIGNAL-TO-NOISE  RATIO.P(dB) 

Comparison  of  suboptimum  detector  performances  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.1,  V2  =  100)  for 


FA 


=  0.1  and  K  =  2 
170 


samp!  es . 


PROBABILITY  OF  DETECTION 


SIGNAL- TO-NOISE  RATIO,  p(dB) 

Figure  4.4-2  Comparison  of  suboptimum  detector  performances  in  bandpass 
Gaussian-Gaussian  mixture  noise  (e  =  0.1,  V  -  100)  for 
p  =  0.01  and  K  =  2  samples. 
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APPENDIX  3-A 


FORTRAN  program  to  calculate  false  alarm  and  detection 
probabilities  for  optimum  detector  in  bandpass  Gaussian- 
Gaussian  mixture  noise. 
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PDP-11  FORTRAN-77  V4.0-1  10:35:01  7-May-85  Paqe  1 

DETECT. FTN;3  /F77/TR : BLOCKS/WR 

0001  PROGRAM  DETPRB 

C 

C  DETECTION  PROBABILITY  CURVES 
C 

C  ANALYSIS:  L.E.  MILLER  IP  APR  85 

C  PROGRAM:  R.H.  FRENCH  19  APR  85 

C  LATEST  VERSION  7  MAY  85 

C 
C 

C  THE  PROGRAM  COMPUTES  THE  DETECTION  PROBABILITY  FOR 

C  THE  OPTIMUM  RECEIVER  FOR  THE  GAUSSI AN-GAUSSIAN  MIXTURE 

C  NOISE  CHANNEL. 

C 

C  THE  COMPUTATIONS  MUST  BE  DONE  IN  DOUBLE  PRECISION  TO 
C  AVOID  UNDESIRABLE  EFFECTS  OF  ROUNDOFF  DUE  TO  THE  SENSITIVITY 
C  OF  DETECTION  PROBABILITY  TO  SMALL  ERRORS  IN  SETTING  THE 
C  THRESHOLDS.  HOWEVER,  THE  PLOT  ROUTINES  REQUIRE  SINGLE 
C  PRECISION  INPUTS;  THUS  EVERYTHING  IS  DONE  IN  DOUBLE  PRECISION 
C  UP  TO  THE  FINAL  ANSWER  ARRAY,  WHICH  THEN  TRUNCATES  THE 
C  RESULT  TO  SINGLE  PRECISION  FOR  THE  PLOT  PACKAGE. 

C 

C  THE  PROGRAM  IS  RATHER  COMPLICATED,  AS  IT  IS  IMPERATIVE  THAT 
C  THE  PROGRAM  AVOID  STARTING  THE  SEARCH  FOR  ROOTS  OF  THE 
C  EQUATION  LAMBDA(X)  =  CONST  IN  THE  REGION  BETWEEN  THE 
C  LOCAL  MAXIMUM  AND  LOCAL  MINIMUM  (IF  THEY  EXIST)  OF  THE 
C  LIKELIHOOD  RATIO.  ALSO,  THE  PROGRAM  MUST  HANDLE  THE  CASE 
C  OF  THE  LIKELIHOOD  RATIO  BECOMING  MONOTONIC  FOR  V**2=l  (THE 
C  GAUSSIAN  CASE)  OR  V**2  .NE.  1  AND  HIGH  SNR.  MUCH  CARE  MUST 
C  BE  TAKEN  THAT  THE  ROOT  SEARCH  NEVER  GOES  BEYOND  A  VALID  POINT 
C  DUE  TO  EVEN  1  LSB  ROUNDOFF. 

C 

C 

C  EVEN  WITH  ALL  THIS  CARE,  A  FEW  IRREGULARITIES  IN  THE  SMOOTH 
C  CURVES  WILL  APPEAR  FOR  SOME  COMBINATIONS  OF  PARAMETERS.  THESE 
C  CAN  BE  MANUALLY  SMOOTHED  OUT  ON  THE  PLOTTED  OUTPUT,  USING  PEN 
C  AND  INK. 

C 

C 

C  BECAUSE  OF  TIME  CONSIDERATIONS,  AS  EACH  CURVE  IS  COMPLETED  THE 
C  DATA  IS  WRITTEN  TO  AN  ARCHIVE  FILE  FOR  POSSIBLE  REUSE  IF  THE 
C  PROGRAM  IS  RESTARTED  AT  A  LATER  DATE. 

C 

C 

0002  IMPLICIT  DOUBLE  PRECISION( A-H,0-Z) 

0003  REAL  EOR I G , EUP I ,PORIG,RTEMP 

C 

C  PLOT  PARAMETERS  xORIG  =  VALUE  AT  ORIGIN,  xUPI  =  UNITS/INCH 
C  WHERE  x  =  E  FOR  SNR  AXIS  AND  x  =  P  FOR  PROBABILITY  AXIS 
C 

0004  PARAMETER  (E0RIG=0. ,EUPI=0.8,P0RIG=-4. ,PUPI=0.5) 

C 

C 
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0006 

0007 
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0009 
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0012 


FORTRAN-77  V4.0-1  10:35:01  7-May-85  Paqe  2 

FTN ; 3  /F77/TR : BLOCKS/WR 

C 

C  PAGE  EJECT  CONTROL  VARIABLES 
C 

CHARACTER* 1  KICKER, ONE 
C 

C  MONO,  M0N02  -  MONOTONIC  CASE  FLAGS  (MONO  FOR  GAUSSIAN  CASE, 

C  M0N02  FOR  HIGH-SNR  CASE) 

C  REUSE  -'  FLAG  FOR  WHETHER  ARCHIVE  FILE  IS  BEING  USED 
C  GOOD  -  FLAG  FOR  REUSE  OF  PREVIOUS  THRESHOLD  (SAVES  TIME) 

C  FIRST  -  FIRST-POINT-OF-CASE  FLAG 

C  ARCFIL  -  ARCHIVE  FILE  NAME  BUFFER.  USE  L0GICAL*1  RATHER  THAN 
C  CHARACTER  ARRAY  BECAUSE  THE  COMPILER  DOES  NOT  PROPERLY 

C  PASS  CHARACTER  ARRAY  TO  OPEN  PROCESSOR  (THE  PROCESSOR 

C  ONLY  SEES  THE  FIRST  WORD  OF  THE  ARRAY  AS  THE  COMPLETE 

C  FILE  NAME) 

C  NUMBRS  -  ARRAY  OF  CHARACTER  CONSTANTS  FOR  CONSTRUCTING  NAME  IN 

C  ARCFIL 

C 

LOG I CAL*1  MONO, M0N02, GOOD, REUSE, FIRST, 

S  ARCFIL(IO) ,NUMBRS(0:3) 

C 

C  PLOT  ARRAYS  AND  TEMPS  WHICH  MUST  BE  SINGLE  PRECISION 
C 

REAL  DBRH0( 153) ,PD( 153) ,DPD,PART,DSNR 

C 

C  LIST  OF  EPSILON  VALUES  TO  RUN 
C 

DIMENSION  EL I ST ( 3 Y 
C 

C  FLAG  FOR  WHETHER  A  GOOD  THRESHOLD  IS  AVAILABLE  FROM  THE 
C  PREVIOUS  POINT 
C 

COMMON  /VALID/  GOOD 
C 

C  PASSES  LN( V**2)  TO  LIKELIHOOD  RATIO  FUNCTION  TO  AVOID 
C  UNNECESSARY  RECOMPUTATION  OF  A  CONSTANT 
C 

COMMON  /LOGCOM/  VSOLN 
C 

C  PARAMETERS  OF  THE  CURVES:  EPSILON,  1-EPSILON,  SNR,  V**2 
C 

COMMON  /PARMS/  EPS,  OME,  RHO,  VSO 

C 

C  SWITCHES  FOR  LIKELIHOOD  RATIO  SHAPE,  MONOTONIC  VS.  NONMONOTONIC. 

C  OEMAX  IS  THE  CRITERION  FOR  DECIDING  THAT  A  MAXIMUM  WILL  NOT  BE 
C  FOUND.  IT  MUST  BE  INITIALIZED  BEFORE  THE  FIRST  LIKELIHOOD  RATIO 
C  IS  EXAMINED,  THUS  IT  IS  MADE  ACCESSIBLE  TO  THE  MAIN  DRIVER  PROGRAM. 
C 

COMMON  /SWITCH/  OEMAX ,M0N02 , FIRST 
C 
C 
C 
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C 

C  PARAMETERS  OF  THE  LIKELIHOOD  RATIO'S  SHAPE: 


C 

ETAMAX  -  THE  VALUE  OF  ETA  FOR  WHICH  LAMBDA(ETA)  ATTAINS 

C 

ITS  LOCAL  MAXIMUM 

C 

ALMAX  -  LAMBDA( ETAMAX) 

c 

ETAMIN  -  THE  VALUE  OF  ETA  FOR  WHICH  LAMBDA(ETA)  ATTAINS 

c 

ITS  LOCAL  MINIMUM 

c 

ALMIN  -  LAMBDA( ETAMIN) 

c 

E3MAX  -  THE  VALUE  OF  ETA,  GREATER  THAN  ETAMIN,  FOR 

c 

WHICH  LAMBDA! E3MAX)  =  LAMBDA( ETAMAX) 

c 

E3PLUS  -  A  SMALL  AMOUNT  BEYOND  E3MAX  TO  INSURE  THAT 

c 

LAMBDA! E3PLUS) -ALMAX  >  0,  REGARDLESS  OF  THE 

c 

SIGN  OF  THE  ERROR  IN  SOLVING  FOR  E3MAX 

c 

r 

MONO  -  MONOTONIC  FLAG 

0013 

L 

COMMON  /FLECT/  ETAMAX, ALMAX, ETAMIN, ALMIN, E3MAX,E3PLUS, 

c 

S 

MONO 

c 

c 

ARCHIVE  FILE  PROTOTYPE  NAME 

c 

THE 

000  IS  CHANGED  TO  THE  INDICES  IEPS, IPFA, IVSQ 

c 

c 

IN  THAT  ORDER,  TO  IDENTIFY  INDIVIDUAL  FILES 

0014 

c 

DATA  ARCFIL/ ' P ' , 1 D 1 , ' 0 ' , ' 0 1 , ' 0' , 1 . 1 , ' D ' , 1 A' , ' T' ,0/ 

c 

CONSTANTS  FOR  USE  IN  ARCFIL 

0015 

DATA  NUMBRS/'O'  ,’1'  ,'2'  ,'3'/ 

c 

LIST 

OF  EPSILONS  AT  IRREGULAR  INCREMENTS 

0016 

DATA  ELI ST/O. 01D0,0.1 DO, 0.5 DO/ 

c 

PAGE- 

-EJECT  CONSTANT  AND  INITIAL  VALUE 

0017 

DATA  ONE, KICKER  /'l' ,'  7 

c 

CALL 

INSTALLATION-STANDARD  RUN  IDENTIFICATION  ROUTINE 

0018 

r 

CALL  JSLGGO 

L 

c 

GIVE 

OPERATOR  A  CHANCE  TO  SELECT  PLOTTER  PEN  TO  USE 

c 

WITH 

RIGHT-HAND  STALL  (#2)  THE  DEFAULT 

0019 

4099 

WRITE(5,4100) 

0020 

4100 

FORMAT!’  ENTER  PEN  TO  USE:  1  FOR  LEFT,  2  FOR  RIGHT  T2l: 

0021 

READ(5,4101)  IPEN 

0022 

4101 

FORMAT! 11) 

0023 

IF( IPEN . EQ.O) IPEN=2 

0024 

IF( IPEN.NE.l. AND . IPEN.NE.2)  GOTO  4099 

C 

INSTALLATION-STANDARD  PLOT  INITIALIZATION  FOR  DIRECT  DRIVE 

C 

OF  THE  HP7470A  PLOTTER 

0025 

CALL  SETUP(KODE) 

0026 

IF(KODE.NE.O)  STOP  'UNABLE  TO  ATTACH  PLOTTER1 

C 

SELECT  THE  DESIRED  PEN 

0027 

CALL  NEWPEN(IPEN) 

C 

C 

C 
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0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 


0045 

0046 

0047 


C 

C  LOOP  ON  MIXTURE  PARAMETER  EPSILON  IS  OUTERMOST  SINCE 
C  IT  IS  A  CONSTANT  PARAMETER  FOR  EACH  SHEET  OF  THE  PLOTS 
C 

DO  900  IEPS=1 ,3 

C  SET  UP  FIRST  DIGIT  OF  ARCHIVE  FILE  NAME 
ARCFIL( 3) =NUMBRS( IEPS) 

EPS=ELIST( IEPS) 

0ME=1 . DO-EPS 

C  INSTALLATION  STANDARD  PLOT  IDENTIFICATION  LABEL 
CALL  PLOTID 

C  REDEFINE  THE  ORIGIN  NEAR  LOWER  LEFT  CORNER  OF  PAGE  IN  NORMAL 
C  REPORT  ORIENTATION,  ALLOWING  ROOM  FOR  ADDING  FIGURE  CAPTION 
CALL  PLOT( 1.25,6.25,-3) 

C  PLOT  A  PROBABILITY-SCALE  AXIS  AS  THE  ORDINATE 
CALL  PR0BAX(0.,0.,' DETECTION  PROBABILITY', 

$  LEN( ' DETECTION  PROBABILITY1 ) ,8. ,0.) 

C  PLOT  A  LINEAR  AXIS  AS  THE  ABSCISSA 

CALL  AXIS(0.,0.,' SIGNAL  TO  NOISE  RATIO  (dB)', 

$  -LENC SIGNAL  TO  NOISE  RATIO  (dB) 1 ) ,6. ,270. ,-10. ,5. ) 
C  ANNOTATE  VALUE  OF  EPSILON  IN  UPPER  LEFT  CORNER  OF  PLOT  AREA 
CALL  SYMB0L(7. 75, -0.25, 0.14, 'EPSILON  =  ',270., 

$  LEN(' EPSILON  =  ')) 

RTEMP=EPS 

CALL  NUMBER(999. ,999. ,0.14,RTEMP,270. ,2) 

C 

C  LOOP  ON  FALSE  ALARM  RATE 
C 

DO  800  IPFAT-1,3 

C  FILL  IN  SECOND  DIGIT  OF  ARCHIVE  FILE  NAME 
ARCFIL(4)=NUMBRS( IPFAT) 

C  TARGET  FALSE  ALARM  RATE 
PFAT=10 . D0**(- IPFAT) 

C 

C  LOOP  ON  VARIANCE  RATIO 
C 

DO  700  IVSQ=0,2 

C  FILL  IN  THIRD  DIGIT  OF  ARCHIVE  FILE  NAME 
ARCFIL(5)=NUMBRS( IVSQ) 

C 

C  TRY  TO  OPEN  ARCHIVE  FILE 
C 

OPEN ( UN IT=2 , NAME=ARCFIL , FORM= ' UNFORMATTED ' , ERR=3000 , 

S  STATUS=' OLD' .READONLY) 

C 

C-  OPEN  WAS  SUCCESSFUL,  JUST  READ  THE  DATA  FROM  ARCHIVES 

C 

REUSE=. TRUE. 

READ( 2)  ISUB,( DBRHO( I) , 1=1 , ISUB+2) , ( PD( I) ,1=1 , ISUB+2^ 
GOTO  3001 


C 

C 
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0048 

0049 


0050 

0051 

0052 


0053 


0054 

0055 

0056 

0057 

0058 

0059 

0060 

0061 

0062 

0063 

0064 

0065 

0066 


C 

C  OPEN  WAS  NOT  SUCCESSFUL,  SO  WE  MUST  COMPUTE  THE  DATA 

C 

3000  REUSE®. FALSE. 

VSQ=10.D0**IVSQ 

C  PASS  LOG  OF  VARIANCE  RATIO  TO  LIKELIHOOD  RATIO  FUNCTION  TO 
C  AVOID  UNNECESSARY  RECOMPUTATION  OF  THE  LOGARITHM 
VSQLN=DLOG(VSO) 

C  PAGE  HEADERS  FOR  PRINTED  OUTPUT  TABLES 
WRITE(6 ,4000)  KICKER, EPS, PFAT,VSO 
4000  F0RMAT(A1,' DETECTION  PROBABILITY  -  PERFECT  RECEIVER'/// 

S  '  EPSILON  =  1  ,F4 . 2 , 1  PFA  =  ' ,1PD9 . 2 , '  V**2  =  ', 

$  0PF6.1///1  RHO  (dB) 1 ,9X , 1  PD 1 ) 

C  KICK  ALL  PAGES  AFTER  FIRST 
KICKER=ONE 
C 

C  LOOP  ON  SIGNAL-TO-NOISE  RATIO 

C  WE  NEED  ABOUT  30  POINTS  PER  INCH  OF  ABSCISSA  TO  GET  A 
C  SMOOTH-APPEARING  CURVE  FROM  THE  PLOTTER 
C 

DO  600  IRH0=1 , 151 

C  INITIALLY  SAY  THE  LIKELIHOOD  RATIO  IS  NONMONOTONIC 
M0N02=. FALSE. 

C  CRITERION  FOR  DECIDING  WE  HAVE  SEARCHED  FAR  ENOUGH 
C  AND  THE  LIKELIHOOD  RATIO  IS  MONOTONIC  IS  THE  LOCATION 
C  OF  THE  LOCAL  MAXIMUM  FOR  THE  PREVIOUS  CASE 
IF( IRHO.GT.l)  OEMAX=ETAMAX 
C  FIRST-POINT  FLAG 

FIRST= IRHO . EQ . 1 

C  THRESHOLD  MIGHT  BE  GOOD  IF  NOT  THE  FIRST  POINT 
GOOD=IRHO.GT.l 

C  STEP  SNR  IN  DECIBELS  -10(0.2)20 
DRTEMP=( IRH0-l)/5. D0-10. DO 

C  BUT  THE  FIRST  POINT  IS  -9.95  DB  RATHER  THAN  -10  DB  TO 
C  AVOID  PROBLEMS  WHEN  SNR  GETS  TOO  SMALL 
I F ( FIRST)  DRTEMP=-9 .95D0 

C  AND  SAVE  SINGLE-PRECISION  VERSION  FOR  PLOT  SOFTWARE 
DBRHO(IRHO)=DRTEMP 
RH0=10 . D0**( DRTEMP/10 . DO) 

ISUB=IRHO 

C 

C  FIND  CRITICAL  POINTS  OF  LIKELIHOOD  RATIO'S  SHAPE 
C 

CALL  CRITIC 
C 

C  SET  UP  INTERPOLATION  CONSTANTS  FOR  THRESHOLDS  NEAR  ETAMAX 
C 

IF( ( .NOT.  MONO)  .AND.  (.NOT.  M0N02))  CALL  TERCON 
C 

C  FIND  THE  THRESHOLD  FOR  SPECIFIED  FALSE  ALARM  PROBABILITY 
C 

CALL  GETETA(PFAT,ETA) 
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0067 


0068 

0069 


0070 

0071 

0072 


0073 

0074 

0075 

0076 

0077 


0078 


0079 

0080 

0081 

0082 

0083 


C 

C  COMPUTE  THE  DETECTION  PROBABILITY 
C 

PD( IRHO)=PDET(ETA) 

C 

C  OUTPUT  RESULTS  TO  PRINT  FILE 
C 

WRITE(6,4001)  DBRHO( IRHO) ,PD( IRHO) 

4001  F0RMAT(1X,F8.4,4X,1PE11.4) 

C 

C  TRUNCATE  THE  PLOT  AT  THE  UPPER  EDGE  OF  THE  PLOT  AREA 
C 

IF(PD( IRHO) .GT. 0.9999)  GOTO  602 
600  CONTINUE 

GOTO  603 
C 

C  INTERPOLATE  TO  EDGE  OF  GRAPH  WHEN  CURVE  GOES  OFF-SCALE 
C 

602  DPD=PD( I SUB) -P  D ( I SUB- 1 ) 

PART=0.9999-PD( I  SUB-1) 

DSNR=DBRHO(  ISUB)  -DBRHO(  ISIJB-1) 

DBRHO( ISUB) =DBRHO( ISUB-1)+PART*DSNR/DPD 
PD( ISUB) =0 .9999 
C 

C  TAKE  INVERSE  GAUSSIAN  DISTRIBUTION  FUNCTION  OF  DATA  FOR 
C  PLOTTING  ON  THE  PROBABILITY-SCALE  AXIS  AND  SET  UP  THE 
C  SCALING  PARAMETERS  (LOCALLY  GENERATED  LIBRARY  ROUTINE) 

C 

603  CALL  PROBSC( PD, ISUB ,8. ) 

C 

C  SET  SCALING  PARAMETERS  FOR  SNR  AXIS;  THE  UNITS/INCH  IS 
C  NEGATIVE  BECAUSE  WE  ARE  PLOTTING  IN  THE  -X  DIRECTION  OF 
C  THE  PLOTTER  HARDWARE  TO  GET  A  GRAPH  ORIENTED  UPRIGHT  ON 
C  THE  PAGE  OF  THE  REPORT 
C 

DBRHO( ISUB+1) =-10 . 

DBRH0( ISUB+2) =-5 . 

C 

C  COME  HERE  IMMEDIATELY  IF  DATA  WAS  READ  FROM  AN  ARCHIVE 
C  FILE 

C 

3001  CONTINUE 
C 

C  DRAW  THE  LINE  BETWEEN  DATA  POINTS 

C 

CALL  LINE(PD,DBRHO, ISUB, 1,0,0) 

C 

C  RAISE  THE  PEN  WHILE  COMPUTING  NEXT  CURVE 
C 

CALL  PENUP 
C 
C 
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0084 

0085 

0086 

0087 

0088 


0089 

0090 


0091 

0092 

0093 


0094 

0095 

0096 


0097 

0098 


C 

C  IF  NO  ARCHIVE  WAS  AVAILABLE,  CREATE  ONE  FOR  THE  NEW  DATA 
C 

IFf .NOT. REUSE)  THEN 

0PEN(UNIT=2,STATUS=I NEW1 ,FILE=ARCFIL ,FORM= 1  UNFORMATTED1 ) 
WRITE( 2)  I SUB , ( DBRHO( I ) , 1=1 , ISU8+2) , ( PD( I ) ,1=1, ISUB+2) 
END  IF 

CLOSE (UNIT=2) 

C 

C  LOOPS  END  HERE 
C 

700  CONTINUE 

800  CONTINUE 

C 

C  START  NEW  PAGE  WHEN  LOOP  ON  EPSILON  INCREMENTS 
C 

IFf IEPS.EQ.3)  THEN 

C  LAST  PAGE,  JUST  FLUSH  PLOT  BUFFER 
CALL  PL0T(0. ,0. ,999) 

ELSE 

C  INSTALLATION-STANDARD  NEW  PLOT  ROUTINE  ASKS  FOR  NEW 
C  PAGE,  WAITS  FOR  SIGNAL  FROM  TERMINAL,  AND  RESETS  THE 
C  PLOT  ORIGIN  TO  INITIAL  LOCATION 
CALL  NEWPLT 
END  IF 

900  CONTINUE 
C 

C  ALL  DONE 
C 

STOP  'DONE' 

END 
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DOUBLE  PRECISION  FUNCTION  PDETl(ETA) 

C 

C  ELEMENTARY  DETECTION  PROBABILITY  FORMULA 

C  FOR  ENVELOPE  SQUARED  GREATER  THAN  A  THRESHOLD  ETA 

C 

C  THIS  USES  MARCUM’S  Q-FUNCTION  FUNCTION  FROM  LIBRARY. 
C  THE  LIBRARY  ROUTINE  IS  BASED  ON  SHNIDMAN'S  ALGORITHM 
C  FROM  IEEE  TRANS.  ON  INFORMATION  THEORY,  VOL.  IT-22, 
l  NO.  6,  NOV.  1976,  PP  746-751. 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

COMMON  /PARMS/  EPS,OME,RHO,VSQ 
QA1=DSQRT( 2 . DO*RHO) 

QA2=DSQRT(2.D0*RH0/VSQ) 

QB1=DSQRT( ETA) 

QB2=DSQRT(ETA/VSQ) 

PDET1=0ME*Q(QA1,QB1)+EPS*Q(0A2,0B2) 

RETURN 

END 
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DOUBLE  PRECISION  FUNCTION  PDET(ETA) 

C 

C  OVER-ALL  DETECTION  PROBABILITY  COMPUTATION 
C 

C  IF  THERE  IS  ONLY  ONE  SOLUTION  TO  LAMBDA(ETA)  =  CONST, 

C  PDET(ETA)  =  PDETHETA) 

C  BUT  IF  THERE  ARE  THREE  SOLUTIONS  TO  LAMBDA(ETA)  =  CONST, 

C  PDEHETA!  =  PDETHETA)  -  PDET1(R00T2)  +  PDET1(R00T3) 

C 

IMPLICIT  DOUBLE  PRECISION( A-H.O-Z) 

L0GICAL*1  UNIQUE 

C  FLAG  FOR  UNIQUE  ROOT  REGION  ON  A  NONMONOTONIC  LIKELIHOOD 
C  RATIO  IS  PASSED  FROM  THE  ROOT  FINDER  WHEN  SETTING  THRESHOLD 
C  TO  MEET  FALSE  ALARM  CRITERION 
COMMON  /WUNRUT/  UNIQUE 
IF! UNIQUE)  THEN 
PDET=PDETlf ETA) 

ELSE 

C  FIND  THE  OTHER  TWO  ROOTS  LAMBDA! R00T2) =LAMBDA( R00T3) =LAMBDA( ETA) 
CALL  ROOTER( ETA, R00T2.R00T3) 

PD1HPDETKETA) 

PD12=PDET1(R00T2) 

PD13=PDET1(R00T3) 

PDET=PD11-PD12+PD13 
END  IF 
RETURN 
END 
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DOUBLE  PRECISION  FUNCTION  PF(ETA) 

C 

C  OVER-ALL  FALSE  ALARM  PROBABILITY  COMPUTATION 
C 

C  IF  THERE  IS  ONLY  ONE  SOLUTION  TO  LAMBDA(ETA)  =  CONST, 

C  PF(ETA)  =  FAP(ETA) 

C  BUT  IF  THERE  ARE  THREE  SOLUTIONS  TO  LAMBDA(ETA)  =  CONST, 

C  PF(ETA)  =  FAP(ETA)  -  FAP(R00T2)  +  FAP(R00T3) 

C 

IMPLICIT  DOUBLE  PRECISI0N(A-H,0-Z) 

LOGICAL*!  MONO 

COMMON  /PARMS/  EPS,OME,RHO,VSQ 

COMMON  /TERROO/  ETA1 , TER001 ,TER002 , PFTER1 , PFTERX 

COMMON  /FLECT/  ETAMAX , ALMAX , ETAMIN , ALMIN , E3MAX , E3PLUS , 

$  MONO 

C 

C  *  *  *  CAUTION*** 

C 

C  THIS  FUNCTION  ASSUMES  ETA  IS  IN  THE  DOMAIN  OF  THE  FUNCTION... 
C  ...  SO  DON'T  CALL  IT  WITH  ETA  BETWEEN  ETAMAX  AND  E3MAX 
C 

R1=RATI0( ETA) 

C 

C  THIS  TEST  FOR  WHETHER  THE  SPECIFIED  ETA  LIES  IN  THE  3-ROOT 
C  REGION  MUST  MAKE  USE  OF  LOCATION  OF  RIGHT-MOST  ROOT  AS  WELL 
C  AS  THE  VALUE  OF  THE  LIKELIHOOD  RATIO  TO  AVOID  PROBLEMS  DUE 
C  TO  ROUND-OFF  AND  ERROR  TOLERANCE  IN  SOLVING  FOR  THE  THIRD 
C  ROOT 
C 

IF(MONO  .OR.  (Rl.LE. ALMIN  .OR.  Rl.GT. ALMAX 
$  .OR.  ETA.GE.E3MAX))  THEN 

C 

C  HAVE  ONLY  ONE  ROOT 
C 

PF=FAP(ETA) 

C 

ELSE  IF (ETAMAX -ETA. GT.l.D-5  .AND. 

$  Rl.GT. ALMIN+l.D-5)  THEN 

C 

C  HAVE  THREE  ROOTS,  NOT  NEAR  THE  LOCAL  EXTREMA  -  FIND  THEM 
C 

CALL  ROOTER ( ETA, R00T2 ,R00T3) 

C 

C  FACTOR  EQUATION  FOR  OVERALL  PF  TO  DEFER  UNDERFLOWS  A  BIT 
C 

PF1=DEXP( -ETA/2. DO)*(1.DO-DEXP((ETA-ROOT2)/2.DO)+ 

S  DEXP( ( ETA-R00T3) /2 . DO) ) 

ETV=ETA/VSQ 

PF2=DEXP( -ETV/2 . D0)*( 1 .DO-DEXP( (ETV-ROOT2/VSQ)/2.DO)+ 
$  DEXP ( ( ETV-R00T3/VSQ) /2 . DO) ) 

PF=0ME*PF1+EPS*PF2 

C 
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C 

C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 

C 


C 

C 

C 

C 

C 


ELSE 

NEAR  EITHER  OF  THE  LOCAL  EXTREMA,  WE  MUST  USE 
SPECIAL  CARE  TO  AVOID  ROUND-OFF  PROBLEMS 
IF( ETAMAX-ETA.LE .  1 .D-5)  THEN 
VERY  CLOSE  TO  THE  LOCAL  MAXIMUM 
IF(ETA.NE.ETAMAX)  THEN 

NEAR  (BUT  NOT  EXACTLY  AT)  THE  LOCAL  MAX,  INTERPOLATE 
THE  PROBABILITY.  THIS  AVOIDS  SEARCHING  FOR  ROOTS  OF 
THE  EQUATION  LAMBDA(ETA)  =  CONST  IN  A  REGION  WHERE 
ROUND-OFF  MAY  MAKE  THE  SOLUTION  IMPOSSIBLE 
EDEL=ETA-ETA1 
FRAC=EDEL/( ETAMAX-ETA1) 

PDEL=PFTERX-PFTER1 

PART=FRAC*PDEL 

PF=PFTER1+PART 

ELSE 

EXACTLY  AT  THE  LOCAL  MAXIMUM  AS  FOUND  BY  SUBROUTINE  CRTTIC; 
TREAT  THIS  AS  A  SPECIAL  CASE  TO  AVOID  ANY  POSSIBILITY  OF 
STARTING  TO  RUN  DOWN  THE  CURVE  TO  THE  LOCAL  MINIMUM  WHEN 
WE  SHOULD  REALLY  JUMP  TO  THE  THIRD  ROOT  AND  BEYOND  WHEN 
SEARCHING  FOR  THE  THRESHOLD  TO  SATISFY  THE  FALSE  ALARM 
CRITERION. 

PF=PFTERX 
END  IF 
ELSE 

NEAR  THE  LOCAL  MIN,  TREAT  AS  IF  AT  LOCAL  MIN 
(HERE  ROUNDOFF  IS  NOT  AS  SEVERE  A  PROBLEM, 

AS  IT  WON'T  CAUSE  THE  SEARCH  TO  GET  STUCK) 

PF=FAP(ETA) 

END  IF 
END  IF 
RETURN 
END 
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DOUBLE  PRECISION  FUNCTION  FAP(ETA) 

C 

C  ELEMENTAL  FALSE  ALARM  PROBABILITY  COMPUTATION 
C 

C  PROB( ENVELOPE  SQUARED  >  THRESHOLD) 

C 

IMPLICIT  DOUBLE  PRECISION( A-H,0-Z) 

COMMON  /PARMS/  EPS,OME ,RHO,VSQ 

FAP=OME*DEXP( -ETA/2. D0)+EPS*DEXP(-0.5D0*ETA/VSQ) 

RETURN 

END 
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DOUBLE  PRECISION  FUNCTION  RATIO(X) 

C 

C  COMPUTE  THE  LIKELIHOOD  RATIO  FUNCTION 
C 

IMPLICIT  DOUBLE  PRECISION A-H.Q-Z) 

COMMON  /PARMS/  EPS,OME,RHO, VSQ 

C  L0G(V**2)  PASSED  IN  FROM  MAIN  PROGRAM  TO  AVOID  RECOMPUTATION 
COMMON  /LOGCOM/  VSQLN 
00VM1=1 . DO/VSQ-1 . DO 

DENOM=OME*DEXP( VSQLN+O . 5D0*X*00VM1)  +  EPS 
BARG=DSQRT(2.D0*RH0*X) 

C 

C  *  *  *  CAUTION  *  *  * 

C 

C  THE  SUBROUTINE  DXBESI  COMPUTES 
C 

C  EXP(-Z)  *  In(Z) 

C 

C  WHERE  In( Z)  IS  THE  MODIFIED  BESSEL  FUNCTION  OF 
C  ORDER  n.  THIS  IS  DONE  TO  AVOID  OVERFLOWS. 

C  THUS  WE  MUST  REMOVE  THE  EXPONENTIAL  WEIGHTING 
C  BY  ADDING  THE  ARGUMENT  OF  THE  BESSEL  FUNCTION  TO 
C  THE  ARGUMENT  OF  THE  EXPONENTIAL  IN  THE  LIKELIHOOD 
C  RATIO'S  FORMULATION.  THIS  AVOIDS  (AS  LONG  AS  POSSIBLE) 

C  THE  SITUATION  OF 

C  (UNDERFLOWING  EXPONENTIAL)  *  (OVERFLOWING  BESSEL  FUNCTION) 
C  =  (GOOD  ANSWER) 

C 

C  ARGUMENT,  ORDER,  RESULT,  ERROR  CODE 

CALL  DXBESI (  BARG,  0,  BESSEL,  KODE) 
IF(KODE.NE.O)  THEN 
WRITE(5,1)  KODE 

1  FORMATC  DXBESI  -  1-  KODE  =  ',12) 

I F ( KODE .NE.3)  STOP  'FATAL  ERROR  FROM  DXBESI' 

END  IF 

PART1=0ME*VS0*DEXP(0.5D0*X*00VM1-RH0+BARG)* BESSEL 
BARG=BARG/VSQ 

CALL  DXBESI(BARG,0, BESSEL, KODE) 

IF(KODE.NE.O)  THEN 
WRITE(5 ,2)  KODE 

2  FORMATC  DXBESI  -  2-  KODE  =  ',12) 

IF( KODE .NE.3)  STOP  'FATAL  ERROR  FROM  DXBESI' 

END  IF 

PART2= ( EPS*DEXP( BARG-RHO/VSO) )* BESSEL 

RATIO=( PART1+PART2) /DENOM 

RETURN 

END 
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DOUBLE  PRECISION  FUNCTION  ROOTF(X) 

C 

C 

C  THE  FUNCTION 
C 

C  RATIO ( X )  -  (TARGET  VALUE) 

C 

C  TO  PASS  TO  THE  ROOT-FINDING  ROUTINES  WHICH  SOLVE 
C  THE  EQUATION  F(X)  =0  FOR  X 
C 
C 

IMPLICIT  DOUBLE  PRECISION( A-H.O-Z) 

C 

C  PASS  IN  THE  TARGET  VALUE  OF  THE  FUNCTION 
C 

COMMON  /TARGET/  FIND 
C 

C  COMPUTE  IT 
C 

ROOTF=RATIO(X) -FIND 

RETURN 

END 
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SUBROUTINE  CRITIC 
C 

C  SUBROUTINE  TO  DETERMINE  SHAPE  OF  THE  LIKELIHOOD 
C  RATIO  AND  ITS  CRITICAL  VALUES 
C 

C  DETERMINES  MONOTONIC  VS.  NONMONOTONIC.  IF  NONMONOTONIC 
C  THIS  ROUTINE  GOES  ON  TO  FIND: 

C 

C  THE  LOCAL  MAXIMUM  (ALMAX)  AND  VALUE  OF  ARGUMENT  FOR 

C  WHICH  IT  OCCURS  (ETAMAX) 

C 

C  THE  LOCAL  MINIMUM  ( ALMIN)  AND  VALUE  OF  ARGUMENT  FOR 

C  WHICH  IT  OCCURS  ( ETAMIN) 

C 

C  THE  ARGUMENT  ( E3MAX) ,  GREATER  THAN  ETAMIN,  FOR  WHICH 

C  LAMBDA( E3MAX)  =  LAMBDA( ETAMAX) 

C 

C  A  POINT  ( E3PLUS)  SLIGHTLY  BEYOND  E3MAX  FOR  WHICH  IT 

C  IS  GUARANTEED  THAT  LAMBDA( E3PLUS)  >  LAMBDAf ETAMAX) , 

L  REGARDLESS  OF  ERRORS  IN  THE  SOLUTION  AND  ROUNDOFF 

C 

C 

C  THE  ROUTINE  ALSO  SETS  UP  THE  MONOTONIC  FLAGS  MONO  AND  M0N02 

C 

C 

IMPLICIT  DOUBLE  PRECISION A-H,0-Z) 

EXTERNAL  ROOTF 
L0GICAL*1  MONO, M0N02, FIRST 
COMMON  /SWITCH/  0EMAX,M0N02, FIRST 
COMMON  /PARMS/  EPS,OME,RHO,VSQ 

COMMON  /FLECT/  ETAMAX, ALMAX, ETAMIN, ALMIN, E3MAX,E3PLUS, 
$  MONO 

COMMON  /TARGET/  FIND 
C 

C  IF  PREVIOUS  CASE  WAS  MONOTONIC,  THEN  THIS  ONE  FOR  HIGHER 
C  SNR  WILL  ALSO  BE  MONOTONIC.  THE  GAUSSIAN  CASE  (V**2  =  1) 

C  IS  ALWAYS  MONOTONIC. 

C 

IF(M0N02 .OR. VSQ . EQ. 1 .DO)  THEN 
C 

C  IF  MONOTONIC,  SAY  IT  IS  SO  AND  WE  ARE  DONE 
C 

MONO= . TRUE . 

RETURN 

ELSE 

C 

C  SAY  IT  IS  NOT  MONOTONIC,  GO  ON  TO  FIND  EXTREMA,  ETC. 

C 

MONO=. FALSE. 

END  IF 
C 
C 
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C 

C  LOCATE  THE  LOCAL  MAXIMUM 

C  (THE  SERIAL  SEARCH  IS  NOT  NECESSARILY  FAST,  BUT 
C  THE  TECHNIQUE  IS  ROBUST  PROVIDED  THE  INITIAL  STEP 
C  IS  SMALL  ENOUGH) 

C 

EINC=0 .IDO 
El=l .DO 

91  R1=RATI0( El) 

1  E2=E1+EINC 

R2=RATI0( E2) 

IF(R2.LT.R1)  GOTO  10 

R1=R2 

E1=E2 

IF ( FIRST)  THEN 

C  IF  THE  FIRST  TIME  AROUND  FOR  THIS  PARAMETER  SET, 

C  SEARCH  FOR  MAX  UNCONDITIONALLY 

C  (THIS  PRESUMES  STARTING  AT  LOW  SNR  WHERE  IT  IS  SURE 
C  TO  BE  NONMONOTONIC) 

GOTO  1 
ELSE 

C  STOP  SEARCH  AT  THREE  TIMES  LOCATION  OF  MAX  OF 
C  PREVIOUS  SNR  AND  DECLARE  IT  MONOTONIC  IF  WE  GET  THIS 
C  FAR  WITHOUT  FINDING  A  LOCAL  MAXIMUM 
IF( E2 . LE . 3 . DO*OEMAX)  THEN 
C  ...  KEEP  ON  SEARCHING 

GOTO  1 
ELSE 

C  ...  TOO  FAR,  IT  MUST  BE  MONOTONIC  (AS  WILL 

C  THE  REST  OF  THIS  RUN  OF  SNR'S) 

M0N02= .TRUE . 

RETURN 
END  IF 
END  IF 

C  STOPPING  CRITERION  IS  1  PART  IN  100,000,000 
10  IF( EINC . LE . 1 .  D-8*E1)  GOTO  20 

C 

C  BACK  UP  ONE  STEP  JUST  IN  CASE  WE  HAVE  THE  FOLLOWING  SITUATION 

C 

C 

C 

C 

c 

C  I  . 

c  .  I 

C  El-EINC  El  E2 

C 

C  WHERE  F( El-EINC)  <  F( El)  AND  F( El)  >  F(E2)  BUT  TRUE  MAX 
C  LIES  BETWEEN  El-EINC  AND  El 
E1=E1-EINC 
EINC=EINC/10 .DO 
GOTO  91 


186 


PDP-11  FORTRAN-77  V4.0-1  10:35:35 

DETECT. FTN;3  /F77/TR : BLOCKS/WR 


7-May-85 


Page  17 


C 

C  HAVE  LOCATED  THE  MAXIMUM 
C 


0037 

20 

ETAMAX=E1 

0038 

ALMAX=R1 

C 

NOW  FIND  THE  MINIMUM 

0039 

EINC=1 . DO 

0040 

2 

E2=E1+EINC 

0041 

R2=RATI0( E2) 

0042 

IF(R2.GT.R1)  GOTO  30 

0043 

R1=R2 

0044 

E1=E2 

0045 

GOTO  2 

0046 

30 

IF(EINC.LE.1.D-8*E1)  GOTO  40 

C 

SIMILAR  REASON  FOR  BACKING  UP  A  STEP 

0047 

E1-E1-EINC 

0048 

Rl-RATIO(El) 

0049 

EINC=EINC/10.D0 

0050 

GOTO  2 

C 

HAVE  LOCATED  THE  MINIMUM 

0051 

40 

ETAMIN=E1 

0052 

ALMIN=R1 

C 

NOW  TO  FIND  SOLUTION  FOR  LAMBDA( E3MAX)  =  LAMBDA( ETAMAX) 

C 

WHERE  E3MAX  >  ETAMIN 

C 

SET  UP  THE  TARGET  VALUE  FOR  ROOT  FINDER 

0053 

FIND=ALMAX 

C 

FIND  BRACKETING  VALUES  El  AND  E2  SUCH  THAT 

C 

LAMBDA(El)  <  LAMBDAf ETAMAX)  <  LAMBDA(E2) 

0054 

EINC=VSO 

0055 

E2=E1*2 . DO 

0056 

3 

R2= RAT 1 0  f E2) 

0057 

IF(R2 .LE .ALMAX)  THEN 

0058 

E1=E2 

C 

NOT  CRITICAL  HOW  CLOSE  THE  VALUES  ARE  TO  THE 

SOLUTION, 

C 

SO  DOUBLE  AT  EACH  STEP  FOR  SPEED 

0059 

E2=E2*2 .DO 

0060 

GOTO  3 

0061 

END  IF 

C 

DO  A  SERIAL  SEARCH  BETWEEN  BRACKETING  VALUES 

0062 

CALL  SERETA( ROOTF , E3MAX , El , E2-E1 , E2) 

C 

GUARANTEE  A  POINT  WHERE  ERROR  IS  POSITIVE 

C 

BUT  NOT  TOO  MUCH  BEYOND  THE  TRUE  ROOT.  THIS 

IS 

C 

REQUIRED  TO  KEEP  THE  SEARCH  FOR  THRESHOLD  TO 

MEET 

C 

FALSE  ALARM  CRITERION  FROM  GETTING  STUCK  GOING  IN 

C 

THE  WRONG  DIRECTION. 

0063 

E3PLUS=E3MAX 

0064 

699  IF(R00TF(E3PLUS) .GT.O.DO)  GOTO  670 

0065 

E3PLUS=E3PLUS*1 .001 

0066 

GOTO  699 

0067 

670  RETURN 

0068 

END 
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SUBROUTINE  R00TER(ETA,R00T2,R00T3) 

C 

C  SUBROUTINE  TO  FIND  THE  SECOND  AND  THIRD 
C  ROOTS  OF  EQUATION  LAMBDA(ETA)  =  CONST, 

C  GIVEN  THE  FIRST  (SMALLEST)  SOLUTION  ETA 
C 

IMPLICIT  DOUBLE  PRECISI0N(A-H,0-Z) 

EXTERNAL  ROOTF 
LOGICAL* 1  MONO 
COMMON  /TARGET/  FIND 

COMMON  /FLECT/  ETAMAX , ALMAX , ETAMIN ,ALM IN , E3MAX , E3PLUS , 

$  MONO 

C  SET  UP  TARGET  VALUE  FOR  ROOT  FINDER 
FIND=RATIO(ETA) 

C 

C  FIND  THE  ROOT  BETWEEN  ETAMAX  AND  ETAMIN 
C  USING  MUELLER'S  ITERATION  SUBROUTINE  FROM  DEC'S  SCIENTIFIC 
C  SUBROUTINE  PACKAGE,  MODIFIED  FOR  DOUBLE  PRECISION 
C 

CALL  DRTM I ( R00T2 , ERROR , ROOTF , ETAMAX , ETAM I N , 5 . D-8 , 150 , KODE) 
IF(KODE.NE.O)  THEN 
WRITE(5,100)  1  ,KODE 

100  FORM AT ( '  DRTM I  ERROR  IN  ROOTER  -  ',11,'  -  CODE  =  ',12/ 

$  '  USING  SERIAL  SEARCH') 

C  IF  MUELLER'S  ITERATION  FAILS,  FALL  BACK  TO  THE  SLOW  BUT 
C  FAIRLY  CERTAIN  SERIAL  SEARCH.  THIS  CAN  NOT  BE  THE  SAME 
C  ROUTINE  AS  USED  TO  FIND  THE  FALSE  ALARM  THRESHOLD,  AS  WE 
C  CAN  NOT  ALLOW  RECURSION  TO  OCCUR  (FALSE  ALARM  THRESHOLD 
C  SEARCH  CALLS  PF  WHICH  IN  TURN  CALLS  ROOTER) 

CALL  SEROOT ( ROOTF , R00T2 , ETAMAX , ETAM I N-ETAMAX ) 

END  IF 
C 

C  FIND  THE  ROOT  BETWEEN  ETAMIN  AND  E3MAX 
C 

C  FIRST  TRY  MUELLER'S  METHOD. 

C  HERE  WE  MUST  USE  THE  POINT  E3PLUS  TO  BE  SURE  THE  BASIC 
C  CRITERION  FOR  MUELLER'S  METHOD  IS  SATISFIED. 

C 

CALL  DRTM I ( R00T3 , ERROR , ROOTF , ETAM IN , E3PLUS , 5 . D-8 , 150 , KODE ) 
IF(KODE.NE.O)  THEN 
WRITE(5,100)  2, KODE 
C  IF  IT  FAILS,  SERIAL  SEARCH 

CALL  SEROOT( ROOTF , R00T3 , ETAMIN , E3PLUS-ETAMIN) 

END  IF 
RETURN 
END 
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DOUBLE  PRECISION  FUNCTION  FINDPF(ETA) 

C 

C  COMPUTE  THE  FUNCTION 
C 

C  PF(ETA)  -  (TARGET  FALSE  ALARM  RATE) 

C 

C  FOR  USE  BY  ROOT  FINDERS  IN  SEARCHING  FOR  FALSE 
C  ALARM  THRESHOLD 
C 

IMPLICIT  DOUBLE  PRECISION! A-H.O-Z) 

C 

C  PASS  IN  THE  TARGET  VALUE 
C 

COMMON  /PFATAR/  FALSEA 
C 

C  COMPUTE  THE  VALUE 
C 

FINDPF=PF( ETA) -FALSEA 

RETURN 

END 
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SUBROUTINE  GETETA( PFAT, ETA) 

C 

C  SUBROUTINE  TO  FIND  THE  THRESHOLD  FOR  WHICH 
C  THE  FALSE  ALARM  PROBABILITY  IS  A  SPECIFIED 
C  VALUE 
C 

IMPLICIT  DOUBLE  PRECISION( A-H.O-Z) 

EXTERNAL  FINDPF 
LOGICAL*!  MONO, GOOD, UNIQUE 
COMMON  /WUNRUT/  UNIQUE 
COMMON  /VALID/  GOOD 
COMMON  /PFATAR/  FALSEA 

COMMON  /FLECT/  ETAMAX , ALMAX , ETAMI N , ALM IN , E3MAX , E3PLUS , 
$  MONO 

C  SET  UP  TARGET  VALUE  FOR  ROOT  FINDERS 
FALSEA=PFAT 
C 

C  IF  THE  PREVIOUS  CASE  WAS  A  SINGLE  ROOT,  THEN  IF  WE 
C  STILL  HAVE  A  SINGLE  ROOT  THE  THRESHOLD  WILL  BE  THE 
C  SAME  AND  WE  DON'T  HAVE  TO  WASTE  TIME  RECOMPUTING  IT. 

C 

'  IF(GOOD. AND. UNIQUE)  THEN 
ALAM=RATIO(ETA) 

C  IF  LAMBDA( PREVIOUS  ETA)  IS  NOT  BETWEEN  THE  LOCAL  MINIMUM 
C  AND  THE  LOCAL  MAXIMUM  VALUES  OF  LAMBDA,  THEN  THERE  IS  ONLY 
C  THE  ONE  SOLUTION  AND  THE  VALUE  OF  ETA  REMAINS  GOOD 
C  FOR  THIS  CASE. 

IF(ALAM.LE .ALMIN  .OR.  ALAM.GE .ALMAX)  RETURN 
END  IF 

1111  IF(MONO)  THEN 
C  IF  MONOTONIC,  ONLY  ONE  ROOT  TO  FIND. 

C  SO  WE  SEARCH  FROM  1.0  ON  UP  TO  NEAR  OVERFLOW  UNTIL 

C  WE  FIND  THE  ROOT,  USING  SERIAL  SEARCH  TECHNIQUE. 

CALL  SERETA( FINDPF, ETA, 1. DO, 1. DO, 1.D38) 

C  AND  FLAG  IT  AS  A  UNIQUE  ROOT 
UNIQUE= .TRUE . 

ELSE 

C  WE  MIGHT  HAVE  3  ROOTS 

C  FIRST  TEST  TO  SEE  IF  BEYOND  THE  3-ROOT  AREA 
C 

PFEMAX=FINDPF( ETAMAX) 

IF( PFEMAX .GT .0. DO)  THEN 
C 

C  WE  HAVE  ONLY  ONE  ROOT  BEYOND  E3MAX  IF  FALSE  ALARM 

C  RATE  IS  TOO  HIGH  AT  E3MAX.  START  SEARCHING 

C  FROM  E3MAX  ON  TOWARDS  AN  OVERFLOWING  VALUE. 

C 

CALL  SERETA( FINDPF, ETA, E3MAX.1. DO, 1.D38) 

C  AND  FLAG  IT  AS  A  UNIQUE  ROOT 
UNIQUE=.TRUE . 

RETURN 
END  IF 
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C 

C  ELSE  WE  ARE  IN  THE  3-R00T  REGION 
C 

IF(GOOD)  THEN 
C 

C  IF  WE  HAVE  JUST  STEPPED  SNR  WITHOUT  CHANGING  V**2  AND  PFA, 

C  THE  PREVIOUS  THRESHOLD  MAKES  A  GOOD  STARTING  POINT  FOR  THE 
C  SEARCH 
C 

XLEFT=ETA 

ELSE 

C 

C  BUT  OTHERWISE  WE  MUST  START  FROM  THE  BEGINNING 
C 

XLEFT=1 . DO 
END  IF 

C  SERIAL  SEARCH  TO  BRACKET  THE  ROOT 
555  YLEFT=FINDPF(XLEFT) 

5  IF( YLEFT.GT.O.DO)  THEN 

C  LEFT  ENDPOINT  FUNCTION  >  0,  SO  INCREASE  RIGHT  HAND  LIMIT 
XRIGHT=XLEFT*2 .DO 

IF (XRIGHT.GT.ETAMAX .AND. XLEFT.LT. ETAMAX)  THEN 
C  RIGHT  END  POINT  MIGHT  BUMP  INTO  ETAMAX,  BE  SURE  WE 
C  DON'T  EXCEED  IT  OR  THE  SEARCH  WILL  GO  WILD 
XRIGHT=ETAMAX 
END  IF 
ELSE 

C  LEFT  ENDPOINT  FUNCTION  <0,  SO  WE  NEED  TO  TRY  A 
C  SMALLER  THRESHOLD  FIRST;  CUT  LEFT  ENDPOINT. 

XLEFT=XLEFT/2.D0 
GOTO  555 
END  IF 

YRIGHT=FINDPF(XRIGHT) 

C  IF  ROOT  LIES  BETWEEN  THE  PROSPECTIVE  ENDPOINTS  WE  ARE  READY 
C  TO  CALL  THE  ROOT  FINDER 

IF( YRIGHT*YLEFT.GE .O.DO)  THEN 

C  BUT  IF  NOT  THEN  MUST  MOVE  THE  INTERVAL  UNTIL  WE  FIND  ONE  THAT 
C  BRACKETS  THE  ROOT.  TO  REDUCE  SEARCH  REGION,  THROW  OUT  THIS 
C  WHOLE  INTERVAL  IN  WHICH  WE  KNOW  THE  ROOT  DOES  NOT  LIE. 
XLEFT=XRIGHT 
YLEFT=YRIGHT 
GOTO  5 
END  IF 

C  DO  THE  SERIAL  SEARCH  FOR  THE  ROOT 

CALL  SERETA( FINDPF,ETA,XLEFT,XRIGHT-XLEFT,XRIGHT) 
ALAM=RATIO(ETA) 

C  DETERMINE  IF  THE  ROOT  IS  UNIQUE 

UNIQUE=ALAM.LE .ALMIN  .OR.  ALAM.GE .ALMAX 
END  IF 
RETURN 
END 
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SUBROUTINE  SERETA(FUNC, ETA, XI, DX, SULIM) 

C 

C  SERIAL  SEARCH  ROUTINE  TO  USE  WHEN  SOLVING 
C  FOR  THRESHOLD  TO  MEET  FALSE  ALARM  CRITERION 
C 

IMPLICIT  DOUBLE  PREC I S I ON ( A-H , 0-Z ) 

E1=X1 

P1=FIJNC(  El) 

I F ( P 1 .EQ.O.DO)  THEN 

C  HAPPENED  TO  HIT  IT  EXACTLY  AT  LEFT  END  POINT 
ETA=E1 
RETURN 
END  IF 
EINC=DX 

1  E2=E1+EINC 

C  AVOID  THE  ROUNDOFF  IN  THE  MAX  SEARCH  LIMIT 
IF(E2.GT. SULIM)  E2=SULIM 
P2=FUNC( E2) 

IF(P2. EQ.O.DO)  THEN 
C  JUST  HAPPENED  TO  HIT  IT  EXACTLY 
ETA=E2 
RETURN 

ELSE  IF(DSIGN(1 .D0,P1) .EQ.DSIGN(1 ,D0,P2) )  THEN 

C 

C  KEEP  STEPPING  (IF  HAVEN'T  BUMPED  INTO  LIMIT;  THIS  IS 
C  NEEDED  TO  AVOID  GOING  OUT  OF  THE  REGION  OF  DEFINITION 
C  OF  THE  FUNCTION  WHEN  SEARCHING  TO  THE  LEFT  OF  ETAMAX) 

C 

IF( E2 . EQ. SULIM!  STOP  'UPPER  SEARCH  LIMIT  TOO  LOW' 

E1=E2 
P1=P2 
GOTO  1 
ELSE 

C  WE  HAVE  IT  BRACKETED 

C  ...  IF  VERY  STEEP  SLOPE,  JUMP  TO  THE  INTERPOLATION  STEP 
I F( ABS( P1-P2)/EINC.GE. 2000 . DO)  GOTO  100 
C  ...  OR  IF  HAVE  IT  TO  ENOUGH  PLACES  JUMP  TO  INTERPOLATION  STEP 
IF( El/EINC. GT .2 . D8)  GOTO  100 
C  ...  OTHERWISE  CUT  INCREMENT  AND  TRY  AGAIN 
EINC=EINC/10.D0 
GOTO  1 
END  IF 
C 

C  INTERPOLATE  BETWEEN  FINAL  TWO  BRACKETING  VALUES  TO  OBTAIN 
C  A  MORE  PRECISE  APPROXIMATION  TO  THE  ROOT 

C 

100  DP=P2-P1 
PP=-P1 

ETA=E1+EINC*PP/DP 

RETURN 

END 
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SUBROUTINE  SEROOTf FUNC ,ETA,X1 ,DX) 

C 

C  SERIAL  SEARCH  FOR  ROOT  OF  EQUATION  FUNC(X)  =  0 
C  FOR  USE  IN  FINDING  ROOTS  OF  LIKELIHOOD  RATIO.  THIS 
C  SECOND  ROUTINE  IS  NEEDED  TO  AVOID  RECURSION. 

C 

C  THE  SEARCH  LOGIC  IS  THE  SAME  AS  ROUTINE  SERETA. 

C 

IMPLICIT  DOUBLE  PRECISION( A-H.O-Z) 

E1=X1 

P1=FUNC( El) 

IF(Pl.EQ.O.DO)  THEN 
ETA=E1 
RETURN 
END  IF 
EINC=DX 

1  E2=E1+EINC 

P2=FUNC( E2) 

IF(P2.EQ.O. DO )  THEN 
C 

C  JUST  HAPPENED  TO  HIT  IT  EXACTLY 
C 

ETA=E2 

RETURN 

ELSE  IF( DSIGN( 1 . DO ,P1) .EO .DSIGN( 1  .D0,P2) )  THEN 
C 

C  KEEP  STEPPING 
C 

E1=E2 
P1=P2 
GOTO  1 
ELSE 

C 

C  WE  HAVE  IT  BRACKETED 
C 

C  ...  IF  VERY  STEEP  SLOPE,  JUMP  TO  THE  INTERPOLATION  STEP 
IF(ABS(P1-P2)/EINC.GE.2000.D0)  GOTO  100 
C  ...  OR  IF  HAVE  IT  TO  ENOUGH  PLACES  JUMP  TO  INTERPOLATION  STEP 
IF(E1/EINC.GT.2.D8)  GOTO  100 
C  ...  OTHERWISE  CUT  INCREMENT  AND  TRY  AGAIN 
EINC=EINC/10.D0 
GOTO  1 
END  IF 
C 

C  INTERPOLATE  BETWEEN  FINAL  TWO  BRACKETING  VALUES 
C 

100  DP=P2-P1 
PP=-P1 

ETA=E1+EINC*PP/DP 

RETURN 

END 
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SUBROUTINE  TERCON 
C 

C  SET  UP  INTERPOLATION  CONSTANTS  FOR  USE  WHEN  THRESHOLD  IS 
C  NEAR  ETAMAX 
C 

C  THIS  IS  A  TIME-SAVER.  DO  IT  ONLY  ONCE  RATHER  THAN 
C  RECOMPUTE  THINGS  WITHIN  A  SUBPROGRAM  WHICH  IS  CALLED 
C  MANY  TIMES  IN  THE  COURSE  OF  SETTING  THE  THRESHOLD. 

C 

IMPLICIT  DOUBLE  PRECISIONf A-H.O-Z) 

LOGICAL*!  MONO 

COMMON  /FLECT/  ETAMAX, ALMAX,ETAMIN,ALM I N,E3MAX,E3PLUS, 
$  MONO 

COMMON  /TERROO/  ETA1,TER001,TER002,PFTER1,PFTERX 
C 

C  BACK  OFF  A  SMALL  DISTANCE  FROM  THE  MAXIMUM 
C 

ETA1=ETAMAX-1 .0-5 

.C 

C  SOLVE  THE  EQUATION  FOR  ROOTS  LAMBDA(ETA)  =  LAMBDA(ETAl) 

C 

CALL  R00TER(ETA1,TER001 ,TER002) 

C 

C  FALSE  ALARM  PROBABILITY  WHEN  ETA=ETAMAX 
C 

PFTERX=FAP(E3MAX) 

C 

C  FALSE  ALARM  PROBABILITY  WHEN  ETA  =  ETAMAX-0 .00001 
C 

PFTER1=FAP ( ETA1 ) -FAP ( TER001 ) +FAP ( TER002) 

RETURN 

END 
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APPENDIX  3B : 

NUMERICAL  TECHNIQUE  FOR  EVALUATING  MULTIPLE  SAMPLE  DETECTOR  PERFORMANCE 


Given  the  detector  form 


U 


Kl— >  =  Z  U  (xk)  S  Z  uk  5  U’ 

k  k 


(3B-1) 


the  probabilities  of  false  alarm  and  detection  are  given  by 


and 


PFA(n)  =  Pr 


PD(n)  =  Pr 


Z-, 


>  n|H0 


k 


Euk > 

k 


(3B-2) 

(3B-3) 


For  independent  jx^j,  and  therefore  independent  |u^|,  the  characteristic 
function  of  the  sum  is 


Vv)  Y\  (3B_4) 

k  K 

which  implies  that  the  probability  density  function  (pdf)  for  the  sum  U  is 
the  K-fold  convolution  of  the  pdf's  for  the  individual  lu^ 


P.|(ct)  =  P„  (<*)  *  P„  (a)  *---*pit  (a).  (3B-5) 

u  U1  u2  UK 

Our  technique  is  based  on  using  a  discrete  pdf  to  approximate  the 
pdf  of  the  individual  detection  statistics,  assumed  to  be  identically  dis¬ 
tributed,  with 
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Pu  (a)  =  Pu(a)  •  £  Pn  6(ct-nA), 
k  n=0 


where 


p  =  PrjnA  <  u  <  (n+l)Aj,  n  =  0, 
n  (  ) 


and  the  number  of  terms  N  and  quantization  step  A  are  traded  off 
the  requirement  that 

N-l 

I  H-i 

n=0 

with  good  precision.  The  K-fold  convolution  of  the  discrete  pdf 
accomplished  iteratively,  with 


min(N-l ,n) 


,(2)  = 


I 


p  p  ,  n  =  0 ,  1 , . . .  ,2N-2 
Fm  Fn-m 


m=max(0,n-N+l) 


,(3)  = 


min(N-l,n) 


I 


p  p^,  n  =  0,  l,...,3N-3 
Fm  Kn-m 


m=max(0,n-2N+2) 


and,  ultimately. 


,00  = 


min(n,N-l) 

I 


p  p(k_^  ,  n  =  0,  1,. 
pm  pn-m 


m=max(0,n-r(N-l) ) 

Using  this  convolved  discrete  density  jpWj  ,  the 
approximately  calculated  using 


(3B-6a) 

(3B-6b) 
to  satisfy 

(3B-6c) 

i  s 

(3B-7 ) 

(3B-8) 

,K(N-1) . 

(3B-9) 

and  Pq  are 
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K(N-l) 


FA 


I 


n=[n/A] 


i(K) 


K(N-l) 


PD  = 


I 


n=[n/A] 


»(K) 


The  different  hypotheses  Ho  and  Hi  are  taken 
calculating  the  original  discrete  set  of  probabilities 


=  PrjnA  <  u^  <  ( n+1 ) a | H - | 
*  Prk  E  RnlHiS' 


where  R  is  the  region  of  x.  which  corresponds  to  nA  < 
n  k 


(3B-6) 

(3B-7 ) 

into  account  by 
according  to 

i  =  0 ,  1 

(3B-8) 

uk=  u(xk)  <  (n+1)A- 
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APPENDIX  4A:  Development  of  alternate  expression  for  noncentral  F 
probability  integral 


ftil  X)  =  e 


-X/2 


r(K+r) 

r(K-l)r(r+l) 


dt  tK"2  (l-t)r 


r(K) 

Hi^ry 


(4A-1) 


where  iF1(a;  b;  x)  is  the  confluent  hypergeometric  function.  Using  (Cummer's 
transformation 


1F1  (a;  b;  x)  =  ex  L Fx  (b-a;  b;  -x) 


(4A-2) 


we  obtain 


f(5.  X)  =  (K-l)  J  dt  tK‘2  e“U/2  jFj  [l-K;  1;  -  \  (1-t)] 

5 


o 

K-l 


=  (K-l)  ^  f 


dt/'2  (l-t)r  <TJt/2 


r=0 


=  (K-l) 


dttK“2  (l-t)r  e‘U/2, 

(4A-3) 
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PART  II: 

AN  INVESTIGATION  OF  CANONICAL  CORRELATION 

AS  AN  AUTOMATIC  DETECTION  AND  BEAMFORMING  TECHNIQUE 

1.0  INTRODUCTION 

1.1  BACKGROUND 

Detection  of  signals  in  noise  can  be  enhanced  by  the  use  of  arrays 
of  sensors.  If  the  N  sensor  outputs  are  delayed  (phased)  such  that  the  signals 
at  each  output  are  in  phase,  the  signal  power  in  the  sum  will  be  proportional 
to  N2.  For  independent  noise  at  each  sensor,  the  noise  power  in  the  sum  will 
be  proportional  to  N.  Thus  an  output  SNR  (signal -to-noise  ratio)  gain  pro¬ 
portional  to  N.  At  the  same  time,  information  on  the  direction  of  the  signal's 
arrival  is  embedded  in  the  time  delays  (phase  shifts).  Figure  1-1  illustrates 
the  array  sum  concept;  it  is  understood  that  the  operations  performed  are 
valid  at  a  given  center  frequency  and  bandwidth  of  interest. 

For  pairs  of  sensors  as  opposed  to  arrays,  the  concept  of  correla¬ 
tion  is  well  understood  when  the  noise  is  Gaussian.  Figure  l-2(a)  shows  how 
the  joint  operations  of  detection  and  time  delay  estimation  can  be  performed 
using  a  correlator.  Although  the  same  operations  can  be  performed  by  squaring 
the  sum  of  the  two  sensor  outputs,  the  sensitivity  of  the  output  to  the  delay 
is  not  as  great.  If  the  background  noise  contains  an  impulsive  component, 
the  usual  assumptions  (isotropic,  uncorrelated  noise)  are,  however,  no  longer 
val id. 

For  more  than  two  sensors,  attention  has  been  given  to  the  relative 
merits  in  Gaussian  noise  of  "standard"  processing  (summing  all  sensors  and 
squaring)  and  of  "multiplicative"  processing,  in  which  partial  array  sums 
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FIGURE  1-1  ARRAY  SUM  CONCEPT 
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(a)  CORRELATION  PROCESSING  OF  TWO  SENSORS 


(b)  CORRELATION  PROCESSING  OF  TWO  ARRAYS 
FIGURE  1-2  CORRELATION  PROCESSING  OF  SENSORS  AND  OF  ARRAYS 
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are  formed,  then  multiplied  [1],  In  general,  the  multiplicative  (correlator) 
configurations,  such  as  shown  in  Figure  l-2(b),  provide  the  same  directivity 
for  fewer  sensors  but  are  slightly  less  sensitive  in  detection.  In  either 
case,  the  task  of  "steering",  or  assigning  delays  or  complex  weights  to  the 
individual  sensors  must  be  made.  When  multiple  sources  are  present,  the 
success  of  both  the  standard  and  multiplicative  approaches  usually  depends 
on  the  notion  of  sweeping  the  "look"  angle,  and  resolution  is  proportioned  to 
the  number  of  sensors. 

Under  the  usual  assumptions  of  Gaussian,  isotropic  ambient  noise, 
it  is  well  known  how  to  combine  two  sensors  for  detection  of  a  signal  of  known 
or  assumed  form;  however,  in  the  presence  of  impulsive  noise  components,  the 
performance  of  this  sensor  application  is  less  well  understood.  How  to  operate 
two  arrays  jointly  is  an  open  subject  with  or  without  Gaussian  noise  assumptions. 

1.2  MOTIVATION  FOR  STUDY 

When  two  arrays  are  separately  located,  the  question  of  combining 
their  sensor  outputs  to  achieve  the  effects  of  a  "super  array"  arises. 

Resolution  is  expected  to  improve  because  of  the  large  baseline.  However, 
the  two  arrays  must  be  time-aligned  or  steered  in  the  correct  initial  directions 
before  any  fine  tuning  can  be  done,  and  it  needs  to  be  ascertained  whether 
they  are  observing  the  same  or  different  sources.  Therefore  correlation  (or 
coherence,  in  the  frequency  domain)  serves  the  purpose  of  confirming  that  the 
arrays  are  looking  at  the  same  source,  as  well  as  actually  performing  the 
detection  and  localization  functions. 

A  new  method  for  cross-correlating  arrays  is  believed  to  be 

l 

promising  .  The  approach  is  to  exploit  the  principles  of  "canonical 
correlation"  [2],  which  may  be  explained  as  follows: 
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Denoting  the  N  =  N  +  N  sensor  outputs  for  the  two  arrays  as  x_  and 
x  y 

the  covariance  matrix  of  the  total  sensor  vector  is 


(1-1) 


*9c 

where  (  )  stands  for  conjugate  transpose  and  i  is  an  N  x  N  =  (N  +  N  )  matrix. 

a  Jr 

Two  steering  vectors  a  (N  x  1)  and  j3  (N  x  1)  can  be  defined;  the  array  sums 

x  y 

become 


w  -  Y.  “ix1  ■  2*£ 

i=l 


2 


Z  6iyi 

i=l 


* 

ir 


(1-2) 


(1-3) 


Canonical  correlation  procedures  find  vectors  a  and  so  that  the  correlation 


E{W2}  =  a*  E(>CX*)jg.  =  a*Zxyi  (1“4) 

is  maximized.  In  fact,  as  many  as  N  or  N  (whichever  is  smaller)  correlations 

X  j 

of  w,z  combinations  can  be  determined  and  ranked. 

Since  the  canonical  correlation  procedure  is  based  on  the  covariance 
matrix  (in  practice,  its  estimate),  its  outputs  (the  steering  vectors  and  the 
correlations)  seem  to  constitute  an  adaptive  solution  to  simultaneously  steering 
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in  multiple  directions.  Thresholding  of  the  correlations  would  then  correspond 
to  detection  of  multiple  targets.  A  multiple  target  capability  would  of  course 
be  of  considerable  utility  in  ocean  surveillance  applications,  especially  in 
view  of  the  consideration  that  the  correlations  can  be  performed  on  data  in 
a  common  spectral  band,  rather  than  at  different  frequencies. 

Another  advantage  of  the  canonical  correlation  method  would  seem 
to  be  that  the  solutions  for  the  steering  vectors  from  information  embedded 
in  the  sensor  covariance  matrix  does  not  require  knowledge  of  sensor  positions. 
Thus  the  method  is  a  form  of  automatic  beamforming. 


1.3  REVIEW  OF  CANONICAL  CORRELATION  METHOD  FOR  REAL  DATA 


As  set  forth  by  Anderson  [2],  the  canonical  correlation 
method  is  the  solution  to  the  following  problem:  Given  the  m-component  zero- 
mean  random  vector  x  with  m  *  m  covariance  matrix  S,  suppose  that  x_  is 
partitioned  into  two  vectors. 


;) 


where  x^  contains  mx  components  and  x 
The  covariance  matrix,  assumed  positive 


contains  m2  components  (mi+m2=  m) • 
definite,  can  be  partitioned  correspondingly 
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where  is  *  m2,  E12  is  irij  x  m2,  E21  is  m2  x  mi >  and  z22  is  m2  x  m2* 

Now  consider  the  arbitrary  linear  combinations  U  and  V,  where 

(1-7) 

and  we  require  that  U  and  V  have  unit  variance,  that  is, 

E  {U2}  =  E  {aTX(1)X(1)Ta}  =  aTElia  =  1  (l-8a) 

and 

E  {V2}  =  E  dV2)x(2)Ti}  =  /e22I  =  L  ( 1  -8b) 

In  (1-7)  and  following  equations  we  use  the  notation  (  )'  to  indicate  the 
transpose  of  a  vector  or  matrix. 

The  correlation  between  U  and  V,  which  have  zero  means,  is 

E  {UV}  =  E  {aV^x'2^Tl>  =  Jzi2  1  .  (1-9) 

The  problem  is  to  find  a  and  §_  such  that  the  correlation  (1-9)  is  maximized 
subject  to  the  constraints  (1-8).  In  the  usual  manner,  the  solution  is  obtained 
by  first  defining  the  function 

y  J 

*  =  a  S12i.  -  kxia  Zu  or  -  1)  -  %(e.TE22_6  -  1),  (1-10) 

where  x  and  y  are  Lagrange  multipliers.  Next,  t  is  differentiated  with 
respect  to  the  elements  of  a  and  _g_,  and  the  derivatives  are  set  to  zero: 
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lit 

3a 


X 1 2  B_  ■  X  X  i  ]  a  =  0 


(l-lla ) 


TT  ~  S12£L  "M2 22 JL  =  £•  ( 1  ~  1 1  b ) 

op_ 

Multiplication  of  (l-lla)  on  the  left  by  aj  and  (1-llb)  on  the  left  by 
yields  the  equations 

^  l  i£L  =  0  ( 1  -12a ) 

BTs|2a  -  ye_TE22i  =0  ( l-12b ) 

Using  the  constraints  (1-8)  in  (1-12),  we  see  that 

X  =  y  =  aTZ12£  ,  (1-13) 


and  (1-11)  can  be  combined  to  form  the  matrix  equation 


(1-14) 


Nontrivial  solutions  of  this  equation  require  that  the  matrix  be  singular: 


det 


XZU 


*12 


-XZ22 


*  0  , 


(1-15) 


\  E21 

or  equivalently,  using  v  =  X2  , 
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det  12^22^2 1  “  vSijJ  —  0  j 


v  =  X" 


(1-16) 


It  can  be  shown  that  if  <  m2  ,  there  are  mx  roots  to  the 
polynomial  in  v  generated  by  (1-16).  The  vectors  a  and  _£  corresponding 
to  these  roots  generate  nij  uncorrelated  linear  combinations. 

If  the  roots  x  =  AT  are  ranked,  then  x.|  (i  =  1,  2,  . . . ,  m  )  is  the 
i  :th  canonical  correlation  between  x.^  and  x.  ^  •  The  vectors  a..  and 
defining  the  linear  combinations  U.  =  alx^  and  V.  =  jjlx'^  satisfy  (1-15) 


for  X  =  X- 


The  conditions  on  the  x  's,  a  's,  and  js  's  can  be  summarized  as 


A  euA  =  I 


^22^1  “  ^ 


A  S  j 2B i“  A  , 


(l-17a) 

( 1 -17b ) 

(l-17c) 


where 


A  A 

•  t 

SLl  * 

ot2 

. a  >  ( m x  m, ) 

.  .-■*  j 

(l-18a) 

<11 

i— H 

CO 

t 

* 

K  * 

e2' 

•  • • • • s  I  >  (m2  x  m, ) 

J 

( 1 -18b ) 

A  A 

diag 

( xi , 

^2  )  *  *  *  A  ) . 

’  mi 

( 1 -18c ) 

210 


J.  S.  LEE  ASSOCIATES,  INC. 


In  addition 


where  B2  is 


if  m1  <  m2,  we  have  the  conditions 


BT  e22B  =  o  ( 1 -1 9a ) 

BT2E22B2=I  (l-19b) 


an  auxiliary  matrix  of  "extra"  J31 s  given  by 


m2  *  (m2 


(1-20) 


211 


J.  S.  LEE  ASSOCIATES,  INC. 


2.0  EXTENSION  OF  CANONICAL  CORRELATION  TO  COMPLEX  DATA 

We  consider  now  extending  the  canonical  correlation  concept  and 
procedure  to  complex  data,  in  order  to  treat  bandpass  signals. 

2.1  FORMULATION  OF  THE  PROBLEM 

Figure  2-1  illustrates  a  situation  in  which  M=2rr.  sensors, 
arbitrarily  located  within  some  area,  receive  one  or  more  signals  arriving 
from  different  directions.  We  assume  that  the  signal  wavefronts  are 
adequately  represented  as  planar,  and  for  simplicity  consider  only  a  two- 
dimensional  case  in  which  the  sensors  lie  in  a  plane.  Each  sensor  also 
samples  the  ambient  noise  background  and/or  generates  within  itself  a  noise 
background. 

The  direction  of  each  signal's  arrival  is  embedded  in  relative 
delays  of  the  arrival  at  the  sensor.  That  is,  for  sensor  i  and  signal  k, 
the  received  waveform  is 

x.j(t)  =  sk(t-T.^)  +  n.(t)  (2-1) 

and  the  direction  of  arrival  information  is  present  in  the  set  of  relative 
(k) 

delays  {t where 

t  00  4  T(k)  .  T.(k> 

V'z  '2  ’»  .  (2-2) 
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Signal  2  arrival  direction 


Signal  1  arrival  direction 


Sensor 


FIGURE  2-1  ILLUSTRATED  SENSOR  PLACEMENT  AND  RECEIVED  SIGNAL  ARRIVALS 
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2.1.1  Complex  representation  of  the  data.  (2-3) 

Bandpass  filtering  of  the  sensor  outputs  is  assumed,  so  that  they 
have  the  form 


x.(t)  =  R.(t)  cos(u)ot  +<f>.(t)}  (2-3a) 

=  x  .(t)  cosw  t  -  x  .(t)  sinto  t  (2-3b) 

Cl  o  SI  o 

=  Re  {(x.c-  +jxsi>  ( 2-3c) 

In  this  formulation  R  is  an  envelope  and  <f>  is  a  phase  with  respect  to  a  center 
frequency  f  =  oj.  / 2tt.  The  quadrature  components  x  .  =  R.cos<j>.  and  x  .  =  R.  sini>- 
are  an  alternate  representation;  samples  of  these  components,  each  at  the  rate 
B  Hz,  where  B  is  the  bandpass  bandwidth,  are  sufficient  data  to  reconstruct 
x  (t)  on  a  given  time  interval.  The  third  form,  (2-3c),  indicates  how 
x-j  (t )  is  related  to  its  complex  envelope  xci(t)  +  jx$..(t). 

In  place  of  x^(t),  we  may  consider  samples  of  the  complex 
envelope  as  constituting  the  data.  Together,  the  M  sensors  produce  the  data 
vector 


x(t) 


Xci(t>  +  j*sl(t) 
xc2(t)  + 


^(t)  +  j>^(t). 


x  (t) 
cm  ; 


jx  (t) 
0  sm 


We  assume  that  the  means  of  ^.(t)  and  (t)  are  zero. 


(2-4) 
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The  covariance  matrix  for  the  sampled  complex  data  vector  is, 
after  Goodman  [4] 

%E{xx*}  =  %E{ (x  +  jx  )  (xl  -  jxl)} 


-c 


s  -c  J-s' 


*  +  4  s 


with  elements 


°ik  *  «xcixck  +  xsixsk  +  jxsixck  •  jxcixsk) 

From  ( 2- 5b )  we  see  that  a^..  =o"ik  the  complex  conjugate  of  cr^  .  Thus 
* 

E  =  Z,  that  is,  E is  Hermitian. 


2.1.2 


Correlation  measure 


Let  the  number  of  sensors  be  an  even  number  M  =  2m,  and  let  the 
data  vector  be  partitioned  into  two  m  *  1  vectors. 


r  d)n 

X 


X  = 


(2) 


the  covariance  matrix  then  can  be  partitioned  as 


E  = 


Z11  e12 
E2i  £22 


where,  since  is  Hermitian,  each  submatrix  is  in  x  m,  and 


(2-5a) 


(2-5b) 


(2-6) 


(2-7a) 


£ll*  =  £n>  £22*  =  ^22»  3nd  £21  =  £12' 


( 2-7  b) 
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A  complex  correlation  between  the  linear  combinations  U  =  a*  x_  and  V  =  x 
may  be  defined  as 


(1)  (2)  , 

E{UV}  =  E{a*  x  x  3> 

=  a*E 321- 

In  general,  this  quantity  is  complex.  The  variances  of  U  and  V,  of  course, 

—  2 

are  real  when  defined  in  the  inner  product  sense  of  EUIU}  =E{|u(  }  and 

2 

E{ | V I  },  respectively. 

2.1.3  Maximum  correlation 

Since  the  correlation  (2-8)  can  be  complex,  the  concept  of  a 
maximum  is  still  ambiguous.  However,  as  we  shall  demonstrate  below,  the 
solution  to  the  maximization  problem  provides  that  the  correlation,  when 
maximized,  is  real. 

The  function  to  be  maximized,  with  constraints,  is 

"k  "k  k 

ip  =  a_  Si2_§.  ~  [a.  ^lia.  “  1]  “  %i[J3  ^22^.  “  • 

First,  consider  differentiation  with  respect  to  the  elements  of  a  =  0  +  ja  , 

R  I 

where  a  and  a  are  the  real  and  imaginary  parts  of  the  vector  a. 

“R  I 

This  operation  yields  the  equations 


M  -  £>2s 


92r 


T_ 

xtina  +  ^na]  =  0 


3iii 

3a 


=j  E3.2JI 


=  0. 


(2-8) 


(2-9) 


(2-10a) 

(2-10b) 
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Premultiplying  (2-10a)  by  aj  and  premultiplying  (2-10b)  by  ,  and  then 

R  I 

adding  the  two  resulting  equations  yields  the  single  equation 


* 

a  £ 1 2_6 


*  T  T  — 

[a_  Ei  ia_  +  a_  Z 1 1  a]  «  0_. 


But  since  the  quantity 


*  *  j  j  — 

a  Ena  =  a  I! ]_ ia_  =  Zn  a_  -  1 


is  real,  (2-11)  reduces  to 


* 

X  =  Z  1 2  S » 

This  demonstrates  that  the  correlation  is  a  real  quantity  if  the  Lagrange 

multiplier  X  is  taken  to  be  real. 

By  a  similar  process,  differentiation  with  respect  to  6  =  B  + j B_ 

R  I 


yields  the  result 

T  T-  * 

vi  =  B_  2*12  ^L  =  5L  ^12^.  = 

Therefore  the  maximization  problem  reduces  to  the  matrix  equation 


S12 

£2i  ~XZZ2_ 

Nontrivial  solutions  to  this  equation  require,  as  in  the  real 

data  case. 


det 


-XEn  e12 
£21  -XE22 


=  0 


5 


(2-11) 


(2-12) 


(2-13) 


(2-14) 


(2-15) 


(2-1 6a ) 
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2 

or  equivalently,  with  v  =  X  , 

-l 

det  [^12^22  ^21  —  =  ®  ’  (2-  16b ) 

2.1.4  Sample  solution  :  single  signal . 

Before  considering  more  complicated  cases,  we  demonstrate  the 
canonical  correlation  method's  solution  for  the  simple  case  of  uncorrelated 
sensor  noises  and  a  single,  sinusoidal  signal.  For  this  case,  the  waveforms 
at  the  sensors  are 

x..(t)  =  s(t  -  t.)  +  n.(t) 

=  RelAe’W*  “V  +  J4>s  +  (nci  +  jnsi)eju)ot  } 

=  Re{^-J*1  *  nc.  +  jns.)e3“ot}  (2-17) 

This  formulation  neglects  any  attenuation  of  the  signal  from  sensor  location 
to  sensor  location.  The  covariance  matrix  for  complex  samples  of  these 
waveforms  is 

* 

Z  =  ^Z{XX  } 

★  T  T 

=  ^E{(Av  +  n  +  jnj  (Av  +  n  -  jn  )> 

—  — c  — s  —  — c  s 

=  h[&  vv  +  Zo  I] 

=  /[I  +  pyv*]  .  (2-18) 
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In  this  development  we  have  used  the  notation 


v 


(2-19) 


and  assume  the  noise  components  have  equal  variances, 

E{nci2}=  E{nsi2}  =  *  '  (2'20) 

with  p  =  A2/2a2  denoting  the  signal -to-noise  ratio  (SNR). 

For  M  =  2m  sensors,  we  partition  v  into  two  m  *  1  vectors  v 

i 

and  v  so  that  the  partitioned  covariance  matrix  has  the  form 
— 2 


z 


2 


a 


I  +  pv  v  * 

— i-i 
* 

pv  V 
-1~1 


pv  V  * 
— I  “2 


I  +  PV  V 


(2-21) 


The  determinant  to  be  solved  for  the  maximum  correlation  is  of  the 


matrix 


-l 

Zi2^22  ^2 l  "  vZll 


=[  pv  v  *( 
“I-! 2 


I  + 


pv  V 

— : 2— 2 


_1 

*)  p 


*2 


■£  *  2 
v  -v(I+Py^y^  )]u 


(2-22) 
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The  needed  matrix  inverse  is 

1  pv  v  * 

(I  +  pv  v  *)  =  I  -  — I2Z2 - 

~2~2  2 

1  +  p|v  | 

2 


Substituting  (2-23)  in  (2-22)  yields  the  determinant 


det 


[<- 


p  v 


2  * 
v)pa  v  v 
-1“1 


2 

va  I 


]■ 


1  +  P  v 

2 


or,  using  the  fact  that  |v  |  =  m, 

det  [(v‘  i+ip)  Qa\^* +  v*2']  =  °- 


It  can  be  shown  by  induction  that 


5 


det  (al  +  bvv_  )  =  am”^  (a  +  mb). 

Applying  this  fact  to  (2-24b)  results  in  the  following  equation  to  be 
solved  for  v: 


(o2v) 


m-i 


2  2 
a  v+  mpo  (v 


.J2£_  )]  =  o. 

1  +  mP  J 


(2-23) 


(2-24a) 


(2-24b) 


(2-25) 


(2-26) 


220 


J.  S.  LEE  ASSOCIATES,  INC. 


The  solution  is  immediately  seen  to  be  that  there  are  m-1  zero  correlations 
(v=0),  and  one  maximum  correlation  (normal  ized)  given  by 


\  =  sfj  = 


mp 


1  +  mp 


Substitution  of  this  value  of  v  in  the  equation 


-i 


E 12^22^2. 1 ®  ~~ 


yields  the  information  that 


,V  1  a 
—  m  '  “i 


kv 


“i> 


where  k  is  a  constant  scale  factor.  The  appropriate  scale  factor  is  found 
from  the  requirement 


or  finally. 


*  2*2  2  * 
a  Zna  =  1  =  k  Vi  (la  +  pa  V  V  )V.i 


=  k  a  m  (1  +  mp) 


Vl 

V  ma2  (1  +  mp) 


2-27) 


2-28a) 


(2-28b) 


(2-29a) 


(2-29b) 
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In  a  similar  manner,  6  is  found  to  be 


v 

“2 


^  V  ma2  (1  +  mp ) 


(2-30) 


For  this  simple  case,  the  solution  is  easily  interpreted.  Maximum 


correlation  between  a  linear  combination  (complex  weighted  sum)  of  signals 
from  half  of  the  sensors  and  a  linear  combination  of  signals  from  the  other 
half  is  achieved  when  the  weights  are  chosen  to  remove  the  relative  propa¬ 
gation  delays  (i.e.,  steer  a  beam  in  the  direction  of  signal  arrival).  For 
example,  a  reconstruction  of  time  samples  of  the  linear  combination  U  would 
produce  the  waveform 


(t)  =  Re{U(t)ejuJ0t} 


u 


(2-31) 


kmAcosojgt  +  Z  n.(t) 


i=l 


In  the  more  general  case,  this  ideal  solution  will  be  affected 
variously  by  complicating  factors,  including 


(a)  the  presence  of  more  than  one  signal 

(b)  non-zero  noise  correlations  between  sensors 


222 


J.  S.  LEE  ASSOCIATES,  INC. 


(c)  attenuation  of  the  signal  as  it  propagates  through  the 
array  of  sensors 

(d)  non-tonal  (modulated)  signals. 
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2.2  SOLUTION  FOR  TWO  SIGNALS 

We  consider  now  how  the  presence  of  more  than  one  signal  affects 
the  canonical  correlation  solution.  The  model  for  the  data  in  this  situation 
is 

N  j6 

-  =  H  Ane  n  vn  +  nc  +  jns  ,  (2-34) 

n=l 


in  which  (An,0n)  are  the  sampled  amplitude  and  phase  of  the  nth  of  N  signals. 
The  covariance  matrix  for  this  data  model  is 


%E(x  x^) 

N 

ozU  +2 

n=l 


p  v  v  *}. 
rt-n— n 


(2-35) 


This  formulation  assumes  that 

je 


E{Ane  nvn  Ake_J’ekvk*}  =  0, 


(2-36) 


that  is,  the  signals  are  uncorrelated. 

2.2.1  General  formulation 

Analytically  we  can  pursue  a  solution  conveniently  for  two  signals. 
For  ease  of  notation,  let  px  =  p,  p2  =  r,  v^v,  and  v  =  w. 


Then  the  covariance  matrix  is 

Y  =  a2{I  +  pv_  v *  +  rww*}.  (2-37) 
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Suitably  partitioned,  this  matrix  has  the  component  matrices  (  all  m  x  m) 


£ll  =  a2(I  + 

pv  v  *  +  rw  w  *) 
l  “l  ~i  “l 

(2-38a) 

Z12  =  a2(pv_i 

*  ,  » 
v  +  rw  w*) 

2  -1—2 

( 2 - 38  b ) 

E2i  =  o2(p12 

i 

11 

(2-38c) 

Z22  =  ct2(I  + 

pv  v*  +  rw  w*) . 

_2—2  —2—2 

( 2-38d) 

The  matrix  whose  determinant  is  to  be  found  is 


-l 


.  vZn 


=  a2(pv  v*  +  rw  w*)  (I  +  pv  v*  +  r w  w*)-^  (pv  v*  +  rw  w *) 

-1  — 2  ~L  “2  “2  “2  ~2  ~2  -2“1  _2“1 


-vo2 (I  +  pv  v*  +  rw  w*) . 

-1— l  -1“1 


(2-39) 


The  inverse  matrix  has  the  form [3] 

(I  +  pv_  v*  +  rw  w*)  ^ 
i  i  i  i 


=  I 


pr[z.v.w.*  +  z.w.v_.*]  -  p(l+mr)v.v.*  -  r(l+mp)w  w* 
+  ill  *1  'll  ii 

1  +  m(p+r)  -  prA. 


using  -  v?  Ej  >  A-j  =  | z .  |2-m2. 


(2-40a) 

(2-40b) 
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After  substituting  (2-40),  the  matrix  in  (2-39)  becomes 
(neglecting  the  factor  a  2) 

[A( z2 )  -  vp]v^*  +  [B(z2)  -  vr]  w^w* 

* 

+  C(z2)  [ZoV  w*  +  Z2W  V  ]  -  vl ,  (2-41) 

1”1  -1“1 

where  z2  is  the  complex  number 


z2  =  v*w 

“2  2 

and  the  coefficients  A,  B,  and  C  are 


A(z2)  = 


B(z2)  = 


p  (m-rA2) 


1  +  m(r+p)  -  prA2 
2 

r(m  -  PA2) 

1  +  m(r+p)  -prA2 


pr 


l  +  m(r+p)-prA2 . 


(2-42a) 


(2-42b) 

(2-42c) 

(2-42d) 


As  a  check  on  our  algebraic  manipulations,  we  note  that  if  r=0,  (2-39)  reduces 
to  the  single  signal  matrix,  (2-24). 


2.2.2  Formulation  for  linear  array. 

When  the  sensors  are  aligned  to  form  a  linear  array  with  equal  spacing 
d,  the  delay  vector  for  the  n:th  signal  assumes  the  special  form 


-je  /2  jms  j(m-l )en 
■  e  [e  ,  e 


-j(m-l)e  r 
e  ] 


(2-43a) 


where 

©n  =  cos  bn  ,  bn  =  bearing. 


(2-43b) 
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using  the  geometry  of  Figure  2-2.  When  partitioned,  we  find  that 


and 


using 


(2) 


-jme  (1) 

e  v 
— n 


W  =  l"Dmt9r%)  exP{— 2 

0  (a)  4  D  (-a)  . 

m  sin(a/2) 


jm(eu-  ej 


}, 


~m '  '  m  sin(a/2)  m 
In  terms  of  the  quantities  introduced  in  Section  2.2.1,  the  linear 
array  geometry  produces  the  relations 


(2-44) 

(2-45a) 
( 2 -4  5  b ) 


=  E22 

(2-46a) 

Zl  =  22  (2-46b) 

~  A2  ( 2-46c ) 

The  overall  covariance  matrix  has  the  form,  for  m  =  2, 

£  =  II  5iic  +  pexp{j(k-i)e1 }  +  r  exp{ j ( k- i ) e2 }  jj  .  (2-46d) 

These  simplify  the  solution  somewhat,  but  further  simplification  does  not  seem 
possible. 

2.2.3  Special  case  of  four  sensors. 

In  order  to  pursue  numerical  cases,  we  consider  a  small  array  size 
of  M  =  2m  =  4.  In  this  case,  we  can  write 
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-(m+1/2)d 


FIGURE  2-2  LINEAR  ARRAY  CONFIGURATION 


m 
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I  S  (1  +  p  +  r) 
11 


I  =  (1  +  p  +  r) 

22 


=1  1 


1 


'2  1 


and 


z  -  (1  +  p  +  r) 
12 


C3 

C5 


c 

c 


4 

6- 


The  canonical  correlations  are  then  the  solutions  to  the  equation 

0  =  X4  (1-|c112)(1-|c2|2) 

2  _  _  _  _ 

+X  {2Re[c2c5Cg  +  c-jC^Cg  +  c^Cg  +  c2c3c4 

-  c1c2c3Cg  -  c1c2c4c5] 

-|c3i2  -  |c4|2  -  |C5|2  -  |c6 I2} 

+X  {  c3|  | c g |  +  |c4|  1  c 5 1  -  c3c4c5cg] L 


(2-47a) 


( 2-47b) 


(2-47c) 


(2-48) 


2 

This  equation  is  quadratic  in  X  =v.  Having  calculated  X,  the  a  steering 
vector  is  found  to  have  components  which  satisfy 


^2 

al 


t c 3 1 2  +  |c4|2-2-Retc2c3c4]  -  X2(l  -|c2|2) 

CT  X  ( 1  -  1 c2 1 ^  “  C3C5  +  C2C3C6  "  C4C6  +  C2C4C6 


and 


a  I  +  la- I  +  2Re[C,a  a  ]  -  1 . 
1  2  1  1  2 


(2-49a) 

(2-49b) 
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These  constraints  determine  a_  to  within  a  complex  factor  with  unit  magnitude, 
so  conveniently  we  may  take  to  be  real. 

Having  found  a,  £  is  determined  by 


3 


X.  £22  ^21  ®  » 

X 


(2-50a) 


or 


and 


Si  =  [(c3-c2c4)a1  +  (cg  -  c2c6)a2]/x(l-|c2|  ) 

_ _  _ _  2 

32  —  f  (c4  "  v-2c3)al  (cg  "  c2c5'a2l  /  ^-(l  -  !G2  I  )• 


( 2- 50b ) 


(2-50c) 


For  the  special  case  of  linear  arrays  with  four  elements,  we  have 

c,  =  c„  =  cr  and  c_  =  c,,  resulting  in  the  simplified  equations 

I  c  b  j  b 


C1 


2  _  3  _ 

+  X  {2Re[2c.|  (c-jC^  +  c^c^)  -  c^  c4] 

-2|c3|2-  2|Ci|2|c3|2  -  |c-|  1 2- i c4 j 2 > 

X  { | c3 1  +  | |  |c4|  -  2Re[c^c32c4'j } 

(2-51) 

a,  <c3|2  +  1C4|2  ~  ^Re^c^l-x'd-lc-,!2) 

ai  C-J - IC1  I2)  -  c.jC3  +  c^Cgl2  -  c3c4-c1c3c4 
_  _  _  _  2 

3i  =  [ (c3-cic4^ai  +  (c]  ~  c^c3)a2]/x(l  -  |c-|  |  )  (2-53a) 

B2  =  [(c^-CjCg)®!  +  (c3-c12)a]  /x(l-|c1|  ).  (2-53b) 
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3.0  NUMERICAL  STUDIES 

Since  an  easily  interpreted  analytical  solution  for  the 
canonical  correlations  and  their  corresponding  steering  vectors  has  not 
been  found,  we  resort  to  selected  numerical  studies  to  explore  the 
dependence  of  these  quantities  upon  various  parameters.  In  all  the 
numerical  cases  presented,  four  sensors  are  assumed  (M=2),  and  the 
array  configurations  illustrated  in  Figure  3-1  were  used.  As  indicated 
in  that  figure,  the  direction  of  arrival  of  the  planewave  signals  is 
given  in  terms  of  the  bearing  n  relative  to  the  x-axis. 

The  parameters  used  in  the  numerical  studies  are  listed  in 

Table  3-1. 

3.1  RESULTS  FOR  SINGLE  SIGNAL 

Although  in  Section  2.1  it  was  shown  that  the  canonical 
correlation  method  yields  an  accurate  steering  vector  for  a  single  signal, 
irrespective  of  array  configuration  or  spacing,  a  number  of  single  signal 
cases  were  calculated  in  order  to  verify  the  analysis  and  also  investigate 
effects  of  noise  correlation. 

According  to  the  analysis  presented  earlier,  the  maximum 
canonical  correlation  value  for  a  single  signal  should  be 

i  -  mP  =  i£ 

A  '  1+mp  1+2P  .  (3-1) 
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►  X 


(a)  Array  configuration  1:  antiparallel  pairs 


(b)  Array  configuration  2:  crossed  pairs 


(c)  Array  configuration  3:  partitional  linear  array 


FIGURE  3-1  ARRAY  CONFIGURATIONS  USED  IN  NUMERICAL  STUDIES 
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Case 

Array 

d/X 

P 

r 

/*"■ 

^2 

correlation 

1 

1 

0.25 

10 

0 

0( .  1)1 

- 

none 

2 

0.25 

1 

0.75*77 

3 

0.25 

5 

0.75*77 

4 

2 

0.25 

10 

0 

0(.1)1 

- 

none 

5 

0.25 

1 

0.75*77 

6 

0.225 

1 

0. 75*77 

7 

0.175 

1 

0.75*77 

8 

0.075 

1 

0.75*77 

9 

0.075 

5 

0.75*7 7 

10 

3 

0.075 

10 

1 

0(  .1)1 

0. 50*77 

none 

11 

0.075 

10 

1 

0 ( .05)1 

Til  _7r/2 

none 

12 

0.25 

10 

1 

0( .05)1 

Til  -77/2 

none 

13 

0,125 

10 

1 

0( . 05)1 

TI-,-7772 

none 

14 

0.075 

10 

5 

0( .05)1 

T|  1  -*77/2 

none 

15 

0.075 

10 

5 

0( . 05) 1 

Tli -*77/4 

none 

16 

0.075 

10 

1 

0( .05)1 

ni-^/4 

none 

17 

0.25 

10 

1 

0(.  1)1 

0.75*77 

none 

18 

0.075 

10 

1 

0(. 1)1 

■  0.75*77 

none 

19 

0.075 

10 

1 

0(. 1)1 

0.75*77 

0.01 

20 

0.075 

10 

1 

0(.  1)1 

0.75*77 

0.10 

21 

0.075 

10 

0 

0( - 1)1 

- 

0.10 

22 

0.075 

10 

0 

0(  .1)1 

- 

0.20 

23 

0.075 

1 

0 

0(  .1)1 

- 

0.20 

24 

0.075 

1 

.1 

0( . 1)1 

0.75*77 

0.20 

25 

0.075 

.1 

.05 

0( .1)1 

0.7577 

0.20 

26 

0.075 

10 

0 

0( .1)1 

- 

none 

27 

0.075 

1 

0 

0( . 1)  1 

- 

0.10 

28 

0.075 

1 

o 

J — 1 

0( . 1)1 

0.75*77 

0.10 

29 

0.075 

10 

1 

0(  .1)1 

0.75*77 

0.20 

Arrays:  1  :  antiparallel  pairs 

2  :  crossed  pairs 

3  :  partitioned  linear 
P  =  SHR1  ,  r  =  SNR2 

TABLE  3-1  PARAMETERS  USED  IN  NUMERICAL  STUDIES 
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In  order  to  asses  the  agreement  of  the  steering  vector 
solution  with  the  actual  direction  of  arrival,  we  define  the  agreement 
metric 

Agreement  -  1  +  cos  (02-<j>2+<(>i ) 

+  cos  (e3-«|)3+<j)1) 

+  cos  )  »  (3-2) 

where  the  steering  solution  vector  is 

(ai,o2>S,,62)  -  (MeJ*1,|a2|eJ^|s1|ej*3jB2|eJN,  (3-3) 


(3-4) 


and  the  actual  received  delay  vector  is 

(ejei!eJ»2,eje3;ejeu).  ■ 

This  metric  ignores  any  constant  phase  difference  between  the  actual 
vector  and  the  solution  vector,  as  well  as  any  differences  in  the  magnitudes 
of  the  components. 

3.1.1  Cases  without  noise  correlation. 

The  elements  of  the  normalized  covariance  matrix  (2-47)  were 
generated  using 

Ci  =  pexp{j(ie1~e2)}/(l+p) 
c2  =  pexp{j("93_94)}/(l+p) 
c3  =  pexp{j(0i-03)}/(l+p) 
c4  =  pexp{j(0i-04)}/(l+p) 
c5  =  pexp{j(02-03)}/(l+p) 

Cg  =  pexp{j(02-04)}/(l+p)  (3-5) 


234 


J.  S.  LEE  ASSOCIATES,  INC. 


and 


(3-6) 


In  the  one-signal  cases  with  no  noise  correlation  (Nos.  1,4,26), 


without  exception  the  canonical  correlation  solution  yielded  one  non-zero 
root  (  X  =  20/21  =  .95238  for  p=  10),  and  the  agreement  metric  was  always 
equal  to  4,  indicating  a  correct  numerical  solution  for  the  delay  vector. 


3.1.2  Cases  with  noise  correlation. 


Several  cases  of  noise  correlation  between  sensors  were  examined 


for  the  linear  array  configuration.  For  that  array  type,  the  correlation 
is  easily  modelled  by  assuming  that  the  covariance  matrix  is,  in  the 
absence  of  signals. 


(3-7) 


I! 


where  a  is  the  correlation  coefficient  (0<a<l)  between  nearest  sensors. 

Table  3-2  gives  the  canonical  correlation  results  as  a  function 
of  signal  bearing,  for  an  SNR  of  p  =  10  and  1,  and  a  noise  correlation 
coefficient  of  a  =  0.1  and  a  =  0.2.  It  is  seen  from  this  data  that  the 
second  canonical  correlation  now  is  nonzero,  but  that  the  noise  correla¬ 
tion  has  increased  the  first  canonical  correlation  values  from  their 
zero-noise  correlation  values  (.95238  for  p  =  10  and  .66667  for  p  =  1). 

The  vector  agreement  metric  indicates  that  the  vector  solution  is  slightly 
degraded  from  a  perfect  value  of  4.0000,  in  proportion  to  the  noise 
correlation,  except  when  the  signal  is  broadside  to  the  array  (n  =  m/2). 
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ti/tt 

a  =  0.1 

a 

=  0.2 

p 

*1 

x2 

agreement 

A.2 

agreement 

10 

0.0 

.95216 

.00722 

3.9750 

.95468 

.01184 

3.8382 

0.1 

.95213 

.00720 

3.9772 

.95462 

.01177 

3.8516 

0.2 

.95207 

.00716 

3.9830 

.95447 

.01158 

3.8886 

0.3 

.95200 

.00710 

3.9908 

.95427 

.01135 

3.9384 

0.4 

.95193 

.00705 

3.9974 

.95412 

.01116 

3.9823 

0.5 

.95191 

.00703 

4.0000 

.95406 

.01109 

4.0000 

0.6 

.95193 

.00705 

3.9974 

.95412 

.01116 

3.9823 

0.7 

.95200 

.00710 

3.9908 

.95427 

.01135 

3.9384 

0.8 

.95207 

.00716 

3.9830 

.95447 

.01158 

3.8886 

0.9 

.95213 

.00720 

3.9772 

.95462 

.01177 

3.8516 

1.0 

.95216 

.00722 

3.9750 

.95468 

.01184 

3.8382 

1 

0.0 

.66716 

.03117 

3.9724 

.67851 

.05038 

3.8663 

0.1 

.66713 

.03109 

3.9748 

.67849 

.05007 

3.8773 

0.2 

.66704 

.03090 

3.9813 

.67846 

.04928 

3.9079 

.  0.3 

.66694 

.03065 

3.9898 

.67844 

.04829 

3.9491 

0.4 

.66686 

.03045 

3.9971 

.67843 

.04748 

3.9854 

0.5 

.66682 

.03037 

4.0000 

.67844 

.04717 

4.0000 

0.6 

.66686 

.03045 

3.9971 

.67843 

.04748 

3.9854 

0.7 

.66694 

.03065 

3.9898 

.67844 

.04829 

3.9491 

0.8 

.66704 

.03090 

3.9813 

.67846 

.04928 

3.9079 

0.9 

.66713 

.03109 

3.9748 

.67849 

.05007 

3.8773 

1.0 

.66716 

.03117 

3.9724 

.67851 

.05038 

3.8663 
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3.2  RESULTS  FOR  TWO  SIGNALS 

For  more  than  one  signal,  our  analysis  did  not  reveal  the 
effects  of  the  various  parameters  upon  the  canonical  correlation  solution; 
only  the  method  for  obtaining  the  solution.  In  this  section,  we  consider 
from  numerical  results  the  effects  of  array  configuration  and  size, 
relative  strengths  of  two  signals,  relative  bearing  separations  of  two 
signals,  and  noise  correlation. 

3.2.1  Effect  of  array  configuration  and  size. 

Assuming  the  array  sensor  spacing  parameter  d  in  Figure  3-1 
is  A/4,  one-quarter  wavelength  of  the  frequency  of  interest,  we 
compare  the  canonical  correlation  solutions  for  the  three  array  config¬ 
urations  in  Figure  3-2  and  3-3,  when  the  first  signal  has  an  SNR  of  p  =  10  dB 
and  the  second  signal  (at  bearing  r\2  =  3ir/4)  has  an  SNR  of  r  =  1  =  0  dB. 

Figure  3-2  shows  that  the  presence  of  the  second  signal, 
though  relatively  weak,  manifests  itself  in  there  being  two  nonzero 
canonical  correlations,  Ax ,  and  A2  ,  except  when  the  bearings  of  the 
two  signals  coincide.  There  is  significant  variation  in  these  correlation 
values  as  a  function  of  the  strong  signal's  bearing,  m  . 
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CANONICAL  CORRELATION  VALUE 


0  .2  .4  .6  .8  1.0 

FIRST  SIGNAL  BEARING,  ^ /tt 

FIGURE  3-2  EFFECT  OF  ARRAY  CONFIGURATION  ON  CANONICAL  CORRELATIONS  VS. 

10dB  SIGNAL  BEARING  WHEN  A  SECOND,  OdB  SIGNAL  IS  PRESENT  AT  A 
BEARING  OF  3*^/4  (d  -  X/4) 
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VECTOR  AGREEMENT  VALUE 


FIRST  SIGNAL  BEARING,  T11  /tt 

FIGURE  3-3  EFFECT  OF  ARRAY  CONFIGURATION  ON  AGREEMENT  OF  STEERING  VECTOR 
SOLUTION  WITH  10dB  SIGNAL  DELAY  VECTOR  WHEN  A  SECOND,  OdB  SIGNAL 
IS  PRESENT  AT  A  BEARING  OF  3t/4  (d  ~X/4) 
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For  the  crossed  pairs  array,  in  Figure  3-2,  the  first  canonical  correlation 
appears  to  be  oscillating  about  the  value 


m(p+r) 


=  0.95 


l+m(p+r) 

while  the  second  canonical  correlation  seems  to  vary  about  the  value 
\2  ~  0.05.  Thus  it  does  not  appear  that  for  this  array  spacing  we 
can  infer  the  strengths  of  the  signals  from  the  canonical  correlation 
val ues. 


(3-8) 


Figure  3-3  reveals  that  for  the  d  =  A/4  spacing  parameter, 
none  of  the  array  configurations  yields  a  very  good  vector  solution', 
except  when  the  two  signals  appear  as  one  signal,  although  the  quality 
of  the  vector  solution  for  the  linear  array  is  high  for  a  wider  range  of 
signal  bearings  than  for  the  other  two  arrays. 

When,  the  array  size  is  decreased,  the  solution  for  the  first 
signal  improves  greatly.  Figures  3-4  and  3-5  demonstrate  for  the  crossed 
pair  array,  that  the  solution  for  canonical  correlations  is  of  more 
consistent  quality  as  bearing  varies  when  the  array  spacing  is  small  (0.3 
times  a  quarter  wavelength).  However,  the  values  of  x2  indicate 
that  this  parameter  does  not  reflect  the  relative  strength  of  the  second 


signal  since 

a2Ai  f  r/p  =  0.1 . 


(3-9) 


240 


CANONICAL  CORRELATION  VALUE 


FIRST  SIGNAL  BEARING,  H-j  Itt 

FIGURE  3-4  EFFECT  OF  CROSSED  PAIR  ARRAY  SIZE  (d/X)  ON  CANONICAL  CORRELATIONS 
VS.IOdB  SIGNAL  BEARING  WHEN  A  SECOND,  OdB  SIGNAL  IS  PRESENT  AT  A 
BEARING  OF  3*^/4 
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VECTOR  AGREEMENT  VALUE 


0  .2  .4  .6  .8  1.0 


FIRST  SIGNAL  BEARING.  n1  lir 

FIGURE  3-5  EFFECT  OF  CROSSED  PAIR  ARRAY  SIZE  (d/X)  ON  AGREEMENT  OF  STEERING 
VECTOR  SOLUTION  WITH  10dB  SIGNAL  DELAY  VECTOR  WHEN  A  SECOND, 

OdB  SIGNAL  IS  PRESENT  AT  A  BEARING  OF  3ir/4 
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The  size  of  the  linear  array  is  varied  in  Figures  3-6  and 
3-7,  and  the  second  signal  is  oriented  so  that  =  -  ir/2  • 

As  for  the  crossed  pair  array,  the  most  consistent  results  occur  for  the 
smallest  array  spacing  tested  (d  =  0.3A/4).  The  least  agreement  of  the 
vector  solution  with  the  first  signal's  delay  vector  occurs  when  m  =  3tt/4 
and  i2  =  tt/4;  at  these  values  the  two  signals  are  located  symmetrically 
with  respect  to  array  broadside  (n=Y2),  and  their  respective  delay 
vectors  are 

V;=  v_  (ex ) ,  w  =  V.  (-0!)  =  v  ,  (3-10) 

that  is,  the  delay  vector  of  the  second  signal  is  the  complex  conjugate 
of  that  of  the  first  signal.  This  relationship  causes  the  covariance 
matrix  for  the  array  to  have  the  elements 

<5jk  +  p  exp{(k-i)0i>  +  r  expl-(k-i )ei^ 

=  6..  W p 2+r2+2rpcos[ 2 ( k-i ) 0 1 ] exp{ j  tan  1  (p~r)sin;t(k-i)9l]_  } 
lk  (p+r)cos  [( k-i ) 0 ! ] 

“  6ik  +  (p+r)exp{j(^-)(k-i)e1},  0!  «  1.  (3-11) 

Thus  when  0i  =  2ird  cos  m/A  is  small,  the  solution  will  approach  that  for 
a  single  signal  at  the  bearing  of  the  strong  signal;  otherwise,  the  bearing 
solution  will  be  distorted  from  the  correct  value. 
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CANONICAL  CORRELATION  VALUE 


FIGURE  3-6  EFFECT  OF  LINEAR  ARRAY  SIZE  (d/X)  ON  CANONICAL  CORRELATIONS  VS. 

1 0dB  SIGNAL  BEARING  WHEN  A  SECOND,  OdB  SIGNAL  IS  PRESENT  AT 
-•7T/2  BEARING  RELATIVE  TO  FIRST  SIGNAL 
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VECTOR  AGREEMENT  VALUE 


0  .2  .4  .6  .8  1.0 


FIRST  SIGNAL  BEARING,  "H -j  /-jr 

FIGURE  3-7  EFFECT  OF  LINEAR  ARRAY  SIZE  (d/A,)ON  AGREEMENT  OF  STEERING  VECTOR 
SOLUTION  WITH  10dB  SIGNAL  DELAY  VECTOR  WHEN  A  SECOND,  OdB  SIGNAL 
IS  PRESENT  AT  -  ir! 2  BEARING  RELATIVE  TO  FIRST  SIGNAL 
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In  Figures  3-8  and  3-9  the  crossed  pair  and  partitioned  linear 
array  solutions  are  compared  for  d  =  0.3A/4  and  the  second  signal  at  n2  =  3^/4. 
Little  difference  exists  between  the  canonical  correlation  values, 
while  the  difference  in  vector  agreement  values  tends  to  reflect  the 
response  patterns  of  the  array  when  the  signal  bearing  is  near  n  =  0. 

3.2.2  Effect  of  relative  signal  strengths  and  angular  spacings. 

The  previous  curves  all  were  for  a  first  signal  strength  of  p  =  10 
and  a  second  signal  strength  of  r  =  1 ,  or  a  10  dB  SNR  ratio  of  p/r.  Now 
we  consider  numerical  cases  for  which  the  SNR  is  varied.  These  cases  are 
presented  in  Figures  3-10  to  3-13  for  the  linear  array. 

From  Figures  3-10  and  3-11  we  observe  that  the  canonical 
correlation  solution  for  the  first  signal  is  degraded  in  proportion  to  the 
strength  of  the  second  signal.  In  these  figures  the  bearing  of  the  second 
signal  was  assumed  to  be  H2  8  ii  -  tt/2 ,  so  that  at  m  =  ir/4  the  two 
signals  appear  as  one  signal  to  the  array  (02  =  0i)  and  at  m  =  3ir/4,  they 
appear  in  opposition,  as  noted  in  the  previous  subsection. 

A  smaller  angular  separation  between  the  signals  was  assumed  in 
Figures  3-12  and  3-13,  namely  n2  =  -  tt/4.  For  these  cases,  the  two 

signals  appear  as  one  signal  for  m  =  ^/8  and  are  opposed  for  ni  =  5V8. 

Again,  the  interfering  effect  of  the  second  signal  is  seen  to  be  in 
proportion  to  its  strength. 
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FIRST  SIGNAL  BEARING,  T-j  Itr 

FIGURE  3-8  EFFECT  OF  ARRAY  CONFIGURATION  ON  CANONICAL  CORRELATIONS 
VS.IOdB  SIGNAL  BEARING  WHEN  A  SECOND,  OdB  SIGNAL  IS  PRESENT 
AT  A  BEARING  OF  3tt/4  (d  =0.3X/4) 
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VECTOR  AGREEMENT  VALUE 


FIRST  SIGNAL  BEARING,  ^  Irr 

FIGURE  3-9  EFFECT  OF  ARRAY  CONFIGURATION  ON  AGREEMENT  OF  STEERING  VECTOR 

SOLUTION  WITH  10dB  SIGNAL  DELAY  VECTOR  WHEN  A  SECOND,  OdB  SIGNAL 
IS  PRESENT  AT  A  BEARING  OF  3^/4  (d  =0.3  A/4) 
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0  .2  .4  .6  .8  1.0 

FIRST  SIGNAL  BEARING,  T11/w 

FIGURE  3-10  EFFECT  OF  SECOND  SIGNAL  STRENGTH  ON  CANONICAL  CORRELATION 
SOLUTIONS  VS.IOdB  SIGNAL  BEARING  WHEN  THE  SECOND  SIGNAL  IS 
PRESENT  AT  -TTI2  BEARING  RELATIVE  TO  FIRST  SIGNAL 
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FIGURE  3-11  EFFECT  OF  SECOND  SIGNAL  STRENGTH  ON  AGREEMENT  OF  STEERING  VECTOR 
SOLUTION  WITH  10dB  SIGNAL  DELAY  VECTOR  WHEN  THE  SECOND  SIGNAL  IS 
PRESENT  AT  -  -jt/2  BEARING  RELATIVE  TO  FIRST  SIGNAL 
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-1 
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0  .2  .4  .6  .8  1.0 
FIRST  SIGNAL  BEARING, Tl1  /tt 
FIGURE  3-12  EFFECT  OF  SECOND  SIGNAL  STRENGTH  ON  AGREEMENT  OF  STEERING 

VECTOR  SOLUTION  WITH  10dB  SIGNAL  DELAY  VECTOR  WHEN  THE  SECOND 
SIGNAL  IS  PRESENT  AT  -  tt/4  BEARING  RELATIVE  TO  FIRST  SIGNAL 
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VECTOR  AGREEMENT  VALUE 


FIRST  SIGNAL  BEARING,  T11 /7r 

FIGURE  3-13  EFFECT  OF  SECOND  SIGNAL  STRENGTH  ON  AGREEMENT  OF  STEERING 

Vector  solution  with  iodB  signal  delay  vector  when  the  second 

SIGNAL  IS  PRESENT  AT  -  -ttM  BEARING  RELATIVE  TO  FIRST  SIGNAL 
(EXPANDED  SCALE) 
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3.2.3  Effect  of  noise  correlation. 

The  results  of  canonical  correlation  solutions  including 
noise  correlation  among  sensors  in  the  partitioned  linear  array  are 
presented  in  Figure  3-14  and  in  Tables  3-3  and  3-4. 

We  observe  from  Figure  3-14  that  the  agreement  between  the 
steering  solution  for  the  first,  stronger  signal  and  its  delay  vector  is 
degraded  in  proportion  to  the  noise  correlation,  as  in  the  case  of  a  single 
signal,  when  the  signals  are  relatively  strong  (10  dB  and  0  dB).  However, 
when  both  signals  are  relatively  weak  (0  dB  and  -10  dB,  -10  dB  and  -13  dB), 
the  vector  agreement,  though  still  affected  adversely  by  the  noise 
correlation,  is  consistently  good. 

The  data  in  Tables  3-3  and  3-4  show  that  the  stronger 
signal's  canonical  correlations  (^i)  are  slightly  increased  by  increasing 
noise  correlation,  while  the  second  canonical  correlation  (^2)  is 
affected  very  little,  being  more  influenced  by  the  relative  SNR  values. 
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VECTOR  AGREEMENT  VALUE 


0  .2  -4  .6  .8  1.0 

FIRST  SIGNAL  BEARING,  V^/ir 

FIGURE  3-14  EFFECT  OF  SENSOR  NOISE  CORRELATION  ON  AGREEMENT  OF  VECTOR 

SOLUTION  WITH  DELAY  VECTOR  FOR  ONE  SIGNAL  WHEN  A  SECOND  SIGNAL 
IS  PRESENT  AT  A  BEARING  OF3ir/4  (LINEAR  ARRAY) 
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TABLE  3-3  CANONICAL  CORRELATIONS  FOR  P  =  10 ,  r  =  1 ,  AND  DIFFERENT  SENSOR  NOISE  CORRELATIONS 
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TABLE  3-4  CANONICAL  CORRELATIONS  FOR  WEAK  SIGNALS  AND  DIFFERENT  SENSOR  NOISE  CORRELATIONS 
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3.3  INTERPRETATIONS  OF  THE  NUMERICAL  RESULTS 

From  our  numerical  studies  in  some  cases  we  may  draw  conclusions, 
while  in  other  cases  the  trends  observed  stimulate  further  questions. 

Among  the  conclusions  which  seem  appropriate  are  the 

following: 

(a)  For  a  single  signal,  the  canonical  correlation  concept 

is  sound.  The  technique  allows  detection  of  the  signal  to  be  based  on  the 
largest  of  the  roots  ( X i ) ,  and  the  vector  solution  corresponding  to 
constitutes  a  beamforming  solution  for  the  direction  of  the  signal's  arrival. 
Thus  detection  can  take  place  without  actually  forming  a  beam,  and  if  the 
relative  sensor  positions  are  known,  the  direction  of  arrival  can  be 
determined. 

(b)  Inter-sensor  noise  correlation  as  high  as  20  per  cent  has 
little  effect  on  the  solution  for  a  single  signal,  for  signals  as  weak  as 
zero  dB  relative  to  the  noise. 

(c)  For  relatively  strong  (10  dB  SNR)  first  signals,  the 
canonical  correlation  solution  is  degraded  severely  when  there  is  a  second 
signal  present,  unless  the  array  spacing  is  small  (less  than  a  quarter 
wavelength)  or  the  second  signal  is  quite  weak  (10  dB  below  first  signal). 

(d)  In  general  inter-sensor  noise  correlation  degrades  the 
canonical  correlation  solution  for  two  signals. 
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Trends  observed  which  are  so  far  inconclusive  with  the 
amount  of  data  generated  include  the  following: 

(a)  Possible  superiority  of  certain  sensor  patterns  over 
others  for  canonical  correlation  solutions;  more  data  is  needed  to 
distinguish  between  the  effects  of  array  configuration  and  sensor  spacing. 

(b)  Possible  good  vector  solutions  for  the  important  case 
of  weak  signals  (see  Figure  3-14);  further  data  is  needed  to  clarify 
the  dependence  of  multiple  signal  solutions  on  signal  strength. 

(c)  Possible  solutions  for  second  signal;  the  second  canonical 
correlation  ( X2)  value  is  proportional  to  the  strength  of  the  second  signal; 
however,  an  attempt  was  not  made  to  calculate  the  vector  solution  for  the 
second  signal. 

In  addition,  based  on  working  with  this  problem  and  exposure 
to  the  numerical  results  obtained  so  far,  we  offer  the  following  conjecture 
concerning  the  performance  of  the  canonical  correlation  technique: 

Conjecture:  The  vector  solution  for  multiple  signals  will 
improve  for  more  than  two  sensors  per  sub-array,  since  angular  resolution 
is  in  general  improved  by  increasing  the  number  of  sensors  and  the  overall 
array  size. 

Support  for  Conjecture:  For  the  case  of  a  linear  array  with 
2m  sensors,  the  inner  product  between  the  direction  vectors  for  two  signals 
is  given  by  (2-45),  in  which  Dm(0i-82)  is  the  half-array  beam  pattern. 
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For  large  m  and  e1  =  02J  this  pattern  response  is  small;  let  us  suggest  the 
approximation  that  the  inner  product  is  negligible.  The  consequence  of  zero 
inner  product  is  seen  by  substituting  z2  =  0  into  (2-41)  and  (2-42),  giving 
the  matrix  in  (2-41)  the  form 


r_mp_ 
p  Ll+mp 


*  +  r 


m  r 
1+m  r 


—i—i 


=  a  I  +  bv_  v_*  +  cw  w;  * 
where,  again  it  is  assumed  that 


0. 


vl 


(3-12a) 
( 3- 12b ) 


Using  ( 3- 12b )  and  (2-25),  the  determinant  of  the  matrix  (3-12a)  whose  roots 


give  the  canonical  correlations  is  found  to  be 


det  [a  I  +  bV  _v  *  +  cw  w 

=  am  det  [I  +  |  VjVj*]  det  [I  +  |  (I  +  f 

=  am_1  (a  +  mb)  det  [I  +  j  (I  -  const,  WjWj*] 

=  am_1  (a  +  mb)  det  [I  +  *  w. w  *]  for  v  *w.  =  0 

=  am~2  (a  +  mb)  (a  +  me).  (3-13) 


This  result  indicates  that  the  two  nonzero  roots  for  v  are  identifiable 
with  the  two  signals: 


.  =  f —  =  mr 

x2  l+mr 

This  gives  confidence  to  expect  a  similarly  good  result  for  the  vector 
solutions  which  accomplish  the  correlations  \\  and  X2  . 


(3-14) 
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4.0  RECOMMENDATIONS  FOR  FURTHER  STUDY 

The  original  motivation  for  investigating  the  canonical 
correlation  technique  was  to  determine  whether  its  application  to  array 
processing  would  yield,  directly  and  automatically, simul taneous  detections 
and  beam  steering  solutions  for  multiple  target  sources.  "Directly"  is 
used  in  the  sense  of  not  requiring  a  physical  beam-steering  mechanism  to 
isolate  a  narrow  range  of  source  directions  for  detection  testing  of  the 
beam  output,  but  rather  a  numerical  procedure  on  the  covariance  matrix  of 
the  sensor  data.  "Automatically"  is  used  in  the  sense  of  the  algorithm's 
not  requiring  knowledge  of  the  sensor  positions. 

The  numerical  results  obtained  in  this  study  were  performed 
for  a  minimum  number  of  sensors  (four)  to  implement  the  concept,  in  order 
to  restrain  the  computational  aspects  of  the  problem.  This  choice  was 
sufficient  to  demonstrate  that  algorithm  works  in  detecting  and  beamforming 
on  a  single  source,  in  agreement  with  the  analysis.  Unfortunately,  this 
small  number  of  sensors  yielded  generally  unsatisfactory  performance  for 
two  sources.  In  Section  3.3,  strong  support  is  given  for  the  conjecture  that 
this  poor  performance  is  due  to  the  small  number  of  sensors  used  in  the 
calculations  and  that  with  a  larger  number  of  sensors,  the  algorithm  will 
successfully  isolate  multiple  sources. 

We  recommend  that  this  study  be  continued,  using  larger  numbers 
of  sensors  and  addressing  the  computational  aspects  of  the  problem.  The 
further  study  should  also  compare  in  some  reasonable  fashion  the  canonical 
correlation  algorithm  with  conventional  techniques  such  as  multiple-beam  or 
beam-scan  in  terms  of  equipment  complexity. 
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